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Microwaves—A Review of Progress in Great Britain 
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A. E. KARBOWIAKT{ 


Summary—tThe fields covered by the review are: Electromag- 
netic theory as applicable to wave propagation, theory of wave- 
guides and components, microwaves in fundamental measurements 
(such as time and atomic constants), designs of components, meas- 
urements, solid-state devices and a summary of utilization. Future 
trends are also indicated. The review (with some exceptions) does not 
cover the following fields: Antennas, propagation and microwave 
tubes. In reviewing technical publications, preference is given to such 
as ects of work which are novel, fundamental or are of controversial 
nature or are likely to influence future trends; where appropriate, 
criticism of the work is given. 


INTRODUCTION 


ICROWAVES is a term (now almost generally 
IV\ accepted) which stands for a vast field of ac- 

tivity. Indeed, so vast that many of us who 
have been active in it for many years have been sur- 
prised by its tremendous rate of growth. Not surpris- 
ingly, therefore, one finds that a substantial proportion 
of published scientific papers, public lectures on science 
and various aspects of engineering (not only electrical) 
as well as subjects on utilization, directly or indirectly 
concern microwaves: microwaves per se, and micro- 
waves for fundamental measurements on atomic and 
material constants. The most precise measurements of 
length and time involve microwaves and vice versa, 
for some microwave uses, the skill of the mechanical 
engineer is tested to the full. Microwave power from 
fractions of a microwatt to many megawatts and fre- 
quencies corresponding to wavelengths from decimeters 
to fractions of a millimeter is now a common practice. 

Industries of diverse interests call on microwaves as 
the new source of power since microwaves, when lib- 
erated in the form of heat, can be applied with precision 
in the quantity and place required. Greater and greater 
powers also become available for communication pur- 
poses and this brings an investigation of industrial, 
biological as well as medical aspects of microwave 
radiation into prominence. 

Parametric amplifiers and masers are the fruits of 
research carried out during the last decade and they are 
often applied in the fields of radio astronomy and com- 
munications. In fact it is said that there has never been 
so much interest in liquid helium as there has been since 
the advent of low noise receivers! 

For applications involving microwaves in communica- 
tion systems, as well as in fundamental research, there is 
a trend towards higher and higher frequencies, yet lack 
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of suitable sources of power in the millimetric and sub- 
millimetric bands besets any rapid progress. Microwave 
tubes are inefficient, have short life and are not a com- 
mercial proposition in the short millimetric band. But, 
for submillimetric frequencies harmonic generators are 
the only sources of monochromatic power. These are 
inefficient and progress is slow and doubtful. Not sur- 
prisingly, new solid-state devices of unspecified nature 
are the hope; the stakes are high but odds are equally 
high. This is a field for fundamental research and so is 
the whole field of generation of coherent radiation in 
the band from submillimetric waves through infrared 
to frequencies corresponding to the visible spectrum. 
Are we going to include this entire field under the aus- 
pices of microwaves, particularly since our many meas- 
urement methods of the higher microwave frequencies 
already resemble those of optics? 

It is fitting to observe that among the many conven- 
tions, conferences and international meetings which 
have taken place in Great Britain during 1960, three in 
particular were concerned, within their own spheres of 
interest to a larger or smaller extent, with microwaves. 
These were: 


1) The XIIIth General Assembly of URSI held at 
University College, London, Eng., September 5— 
15, 1960. 

2) A Conference on Solid State Microwave Ampli- 
fiers held at Nottingham University, Eng., April 
6-8, 1960. 

3) The 3rd International Conference on Medical 
Electronics held at Olympia, London, July 21-27, 
1960, under the auspices of the IEE and the In- 
ternational Federation for Medical Electronics. 


Work in the microwave field can conveniently be re- 
viewed under the following headings: 


1) Investigation of fundamentals of physics such as 
measurements of length and time and atomic and 
material constants. 

2) Investigation of fundamentals of electromagnetic 
waves. 

3) Development in microwave components and meas- 
urements. 

4) Generation of microwave radiation, in particular 
solid-state devices. 

5) Propagation and antennas. 

6) Utilization including industrial applications and 
communication systems. 

7) Future trends. 
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The above topics cover an immense territory of our 
knowledge of which only a part was chosen for this 
review. 
| Under item 1 a great deal of work done in physics, 
‘such as spectroscopy, could be listed but because of its 
very specialized nature it will be omitted from further 
discussion. On the other hand the measurements of 
length and time deserve a mention, because such meas- 
urements are fundamental. 

Study of microwaves and waveguides has always in- 
volved a great deal of electromagnetic theory and not 
infrequently progress in the field of microwaves was a 
direct result of deeper understanding of electromagnetic 
theory; a step forward in the electromagnetic theory is 
always a step forward in the field of microwaves. 

Progress in microwave components and measure- 
ments is usually of the nature of an improvement or ex- 
tension in the existing knowledge and consequently no 
detailed discussion of such work will be given here. 
Nevertheless, such contributions are building-blocks of 
microwave techniques and some references are there- 
fore included. 

Any progress in the generation of microwave power, 
particularly in the upper frequency band is of utmost 
significance in further progress in the microwave field; 
but progress in conventional microwave tubes has 
been extensively reviewed recently by Harvey [1]. 
This review will be chiefly concerned with generation of 
microwaves using solid-state devices, a fascinating field 
full of vague promises. 

Propagation and antennas is a subject intimately re- 
lated to the activities of URSI and progress in this field 
can be gathered from the many excellent reports of the 
XIIIth General Assembly, which took place in the 
autumn of 1960. Mention should be made however of 
the work by Belatini [2] which does not seem to have 
received any publicity. The work which concerns 
tropospheric scatter phenomena, deserves attention be- 
cause there exists a considerable disagreement over the 
exact process of the scattering phenomenon. Belatini 
observes that the field received beyond the horizon 
shows an almost regular spatial interference pattern, 
implying that incoherent scatter cannot be the mecha- 
nism responsible, as is so widely accepted, for signals 
received beyond the horizon. The author’s thesis is that 
the heat from the sun is responsible for vertical motion 
of air producing in effect torus-shaped “Bénard cells.” 
These cells act as weak biconcave lenses causing the 
radio waves to diverge in the vertical plane and in this 
way enable radio signals to reach beyond the horizon. 
The hypothesis is not only original, but the evidence 
against “tropospheric scatter” is strong and the evi- 
dence for Belatini’s theory good. 

Industrial applications, and communication systems 
in particular, are the important uses for microwave 
equipment. Much of the interest and financing of re- 
search in industries is taking place because of such 
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applications. The existence of research, therefore, hinges 
to a large extent on possible industrial and commercial 
applications although research per se is being accepted 
on a larger and larger scale. Microwave links are, to a 
large extent, the main commercial interest, but, in 
addition, numerous defence projects such as radar, etc., 
occupy the activity of many. But new communication 
systems such as long distance communication by wave- 
guide and satellites occupy the imagination of many. 
Some other industrial applications have been discussed 
adequately by Harvey [3]. 

Assessment of future trends, which is of great im- 
portance to industry and those responsible for plan- 
ning future activities, involves a specific form of crys- 
tal-gazing. But, there is one certain lesson that we can 
learn from the past: few things in life are less predict- 
able than the activity of a research worker; in fact the 
game is so speculative that staking a Derby winner is 
not an unjust analogy! 

This review is primarily concerned with work which 
in the opinion of the reviewer is either of a fundamental 
nature or is likely to have a significant impact on the 
progress in microwaves, has features of novelty or is 
simply fascinating. 


FUNDAMENTALS 


A concise review of progress in electrical units and 
standards has been prepared by Vigoureux [4]. It 
transpires that the fundamental units of electricity and 
magnetism can now be measured with a reproducibility 
of about 1 part in 10° using electrical measurements on 
constants such as gyromagnetic ratio of the proton. 

Length and time—the two fundamental units—can 
be measured and the units maintained and reproduced, 
to a much higher order of accuracy than was possible 
only a few years ago. The measurements are carried 
out in electromagnetic spectrum; but whereas the 
“standard of length” is conveniently found in the visible 
spectrum (for example, the wavelength of the orange 
line of krypton—86), the most accurate “clock” op- 
erates at a fundamental frequency falling in the micro- 
wave band. The caesium frequency standard which 
operates at about 9.2 Ge, is a result of three or four years 
of research carried out in various countries. 

In Great Britain there are now a number of caesium 
clocks in operation following the fundamental work 
carried out by Essen and Parry at the National Physical 
Laboratory [5]. Such clocks operate with an error not 
exceeding about 1 part in 101°. 

In March, 1959, the International Committee of 
Weights and Measures adopted the second of Ephemeris 
Time as the official standard. The astronomical stand- 
ard, however, cannot be reproduced to better than about 
1 part in 10°, even when averaged over a period of sev- 
eral years. The caesium beam standard, on the other 
hand, is readily available and there is a possibility that 
the second of time may be defined in terms of frequency 
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of an atomic frequency, and this lies in the microwave 
band. 

Quartz crystal oscillators are still a very convenient 
intermediary between the caesium clock and _ astro- 
nomical observations: But with the recent develop- 
ments in the atomic frequency standards and also 
MSF transmissions in the United Kingdom and Western 
Europe, the quartz crystal oscillators are likely to be- 
come superfluous. 

Present work on atomic frequency standards is con- 
cerned chiefly with design of more compact units (some 
so small that they can be fitted inside rockets or in- 
stalled inside an artificial satellite), study of different 
designs, sources of errors and intercomparison of meas- 
urements. There is at present an unresolved discrep- 
ancy of about 2 parts in 10'° between clocks of different 
type or build in different countries [6]. 

A characteristic feature of many meetings or lectures 
concerning microwaves is that the subsequent discus- 
sions are frequently diverted towards the fundamentals 
of electromagnetic theory. There again one is frequently 
involved in discussion on basic units [7], or even the 
choice of coordinate systems and nomenclature. The 
opinions are frequently diverse and sometimes no con- 
clusions are reached. This is all to the good, yet never- 
theless a little surprising, for it is an outward sign of an 
inward uneasy feeling about the fundamental concepts. 
Most of us, at one time or another, when engaged in an 
investigation of an advanced nature have experienced a 
feeling of complete standstill, only to find out that the 
root of the trouble was inadequate comprehension of 
fundamentals and not the apparent complexity of the 
particular problem on hand. We should, rightly so, 
have a look at our fundamental formulas and examine 
them from time to time in the light of gathered experi- 
ence. It is reassuring to find that this spirit is prevailing 
and our concern about the fundamental concepts is a 
healthy one. 

The Department of Scientific and Industrial Research 
has issued a booklet containing tables [8] of the refrac- 
tive index of air for a wide range of temperatures, pres- 
sures and water-vapor pressures based on formulas de- 
rived by Essen and Froome. 


THEORY OF WAVE PROPAGATION AND WAVEGUIDES 


The foundations of waveguide theory were laid down 
more than half a century ago by Lord Rayleigh, but it 
was the widespread adoption of the powerful analytical 
techniques from the field of quantum mechanics di- 
rectly to microwaves that led to the elaboration of 
waveguide theory as we know it today. In principle, 
therefore, any imaginable wave-guiding structure could 
be analyzed in detail, yet it would be a gross over- 


. The implied great faith in the stability of atomic vibrations 
which we seem to cherish, cannot be overemphasized, particularly 
since our knowledge in connection with atoms and molecules is 
still so new. 
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statement to say that there is fundamentally nothing 
new that could be discovered. For there is more to un- 
derstanding than mere writing down of general mathe- 
matical expressions. One must be able to comprehend 
what the mathematical expressions stand for and be 
able to formulate useful physical pictures before a clear 
physical interpretation can be given; and later on one 
must be able to apply it to the design of a particular 
structure, be it a single component or a complete com- 
munication system. It is clear that during the last few 
years we have progressed towards that goal, but we 
have not reached it entirely. Indeed, it is doubtful if 
such an ideal state of understanding can ever be 
reached. But, it is only through discussion, analysis, 
and a continual examination and re-examination of our 
accumulated facts that a progress towards unification 
of our knowledge can be made. At the same time we 
may be able to derive some enjoyment from it. 

There are two outstanding aspects of the theory of 
waveguides: 1) propagation in waveguides with walls 
other than perfectly conducting, and 2) propagation in 
nonuniform waveguides. 


Propagation in Homogeneous Uniform Waveguides 


The two above mentioned factors have one thing in 
common: we do not seem to be entirely sure what the 
best way is to describe an electromagnetic field inside a 
waveguide and how such fields should be classified. A 
particular case is that of surface waves.? 

Some contributors, in their own way, maintain that a 
rational unified and useful classification of waveguide 
modes is possible. Thus, Karbowiak [9] shows quite 
generally that the electromagnetic field can be con- 
veniently resolved into pole waves and_ branch-cut 
waves; but in addition he also introduces the concept of 
a quasi-mode. The findings can be briefly summarized 
as follows. 

To every cylindrical waveguide, there corresponds 
[10] a certain infinite but discrete distribution of poles 
in the complex plane of the propagation coefficient, y, 
and one or two branch points (usually only one) de- 
pending on the number of geometrically possible in- 
finite directions, in the transverse plane of the wave- 
guide. Integration around the branch-cut yields the 
branch-cut wave, and the poles represent the possible 
modal (or guided wave) propagation. This picture is 
agreed on by most workers in the field. Karbowiak 
however accepts for the definition of a pole wave the 
following function: ;=a;-f,(u, v) -A;-exp ( — iz), where 
a; is the amplitude coefficient, f,(w, v) a function 
describing the field pattern in the waveguide trans- 
verse plane and A; a multiplier, which is related to the 


* It was the function of a Study Group of Commission VI of URSI 
to try to discuss various types of surface waves with a view to classi- 
fying them. But, it transpired that there are nowadays so many 
variants of surface wave concepts that a useful definition could not 
be formulated. It was agreed that to pursue the matter of terminology 
applied to surface waves would be unprofitable. 
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complex Fresnel Integral. This definition differs from 
the more commonly accepted one by the multiplier A. 
If A=1, then the expression defines a quasi-mode. The 
point of such a definition will transpire from an ex- 
‘ample. 

Imagine a waveguide with a given geometry and ex- 
citation having arbitrary but physically realizable 
walls. Then in the plane of y there will be a distribution 
of an infinity (discrete) of poles and, say, one branch 
point. Of the poles, some will be placed on the right and 
some on the wrong Riemannian Plane [10]. The quan- 
tity A will be close to unity for all poles on the Sree? 
plane and it will be very nearly zero for poles lying on 
the “wrong” plane. But, A is a continuous function of 
frequency and parameters describing the physical 
medium and, therefore, if either of these quantities is 
~ gradually varied, then the poles will move from place to 
place and the value of the multiplier will vary. Some 
poles, however, will pass gradually through the branch 
cut onto the other Riemannian plane, but the total 
number of poles will remain unaltered, and, therefore, 
all pole waves exist provided the excitation is such that 
a; is finite. The constitution of the field changes con- 
tinuously, because depending on the value of the multi- 
plier A some pole waves become more significant than 
others. If we were to expand the field in terms of func- 
tions with A=1 then clearly the modes would either 
exist or not depending on whether the pole were just to 
the left or the right of the branch cut. Since all events 
in nature are continuous such a classification could not 
be recommended. Yet, functions with A =1 are useful in 
describing fields in finite regions and for such modes the 
term quasi-mode is reserved; modes of this kind need 
not, on their own, be a solution to the waveguide prob- 
lem, but in finite regions (of interest) the field may very 
closely approximate such modes and hence their signif- 
icance. A leaky wave is one example, a Zenneck wave 
another. 

Some investigators, however, do not share Karbo- 
wiak’s viewpoint. Barlow [11] for example, proposes to 
classify electromagnetic waves into three types: 1) 
TEM waves, 2) surface waves and 3) so-called wave 
guide modes. Such classification, while having the merit 
of apparent simplicity, would be rather restricted, as 
pointed out recently by Waldron [12]. Because, if we 
consider the pole representing (say) the TEM wave in 
a perfect waveguide, then as we change gradually the 
perfectly conducting medium and make it imperfect, 
the pole will move slowly away from the imaginary 
axis and the associated A will become a little less than 
unity. Eventually, if the properties of the media are 
sufficiently altered, then the mode will not in the least 
resemble the TEM mode. The wave due to that particu- 
lar pole may even be more like some other mode or even 
form a hybrid with one of its neighbors, etc. As a pole 
wave, however, it exists all the time, though for some 
properties of the media it may appear as one quasi-mode 
and for some as another. 
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Surface waves in various forms have always been a 
fruitful field for research—surface waves for transmis- 
sion lines or for antennas, surface waves for measure- 
ments on materials and surface waves for their own 
sake. 

Clarricoats has made an interesting and exhaustive 
study of wave propagation along dielectric rods [13]. In 
the past, several investigators have analyzed propaga- 
tion of selected modes on dielectric cylinders with and 
without a concentric metallic shield. In this study, 
however, the author was able, through a detailed anal- 
ysis, to show how the various modes (which sometimes 
are very complicated) are related to each other and 
how the mode spectrum is affected by the cross-sec- 
tional dimension of the rod and the proportions of the 
metallic shield. This is yet another example of a uni- 
fication of the concept of modes and a greater under- 
standing of the structure of electromagnetic field. A 
graphical method of evaluation of propagation co- 
efficients is given and application of the results to ferrite 
devices is also indicated. 

The particularly interesting feature of the analysis 
is the physical picture. The mode spectrum of the struc- 
ture can be followed step by step as the diameter of the 
dielectric rod and the diameter of the metallic shield are 
varied. In fact, the modes of an empty cylindrical 
waveguide are linked to those of a shielded dielectric rod 
and these in turn are related to those of an unbounded 
dielectric rod, the surface wave modes. It would be in- 
teresting to analyze the excitation problem in such 
waveguides. From such an analysis, we venture to 
guess, it would transpire that there is a branch-cut wave 
representing the radiation field, and a discrete infinity 
of pole waves. The pole waves would always be present - 
but their proportion in the total field and their charac- 
ter would vary as the proportions of the dielectric rod 
and the metallic shield were varied. In this way modes 
of one waveguide could be related to modes in another 
waveguide by a continuous process without any am- 
biguity. Furthermore, over certain finite regions some of 
the modes (quasi-modes) would be appreciable and 
some other modes would be negligible (the quasi-mode 
would have no physical significance). 

A number of other surface wave modes and support- 
ing structures have also been analyzed. The hybrid EHu 
mode of a dielectric rod has received detailed attention 
with a view to utilizing this mode to measurements of 
permittivity of samples in the form of circular rods [14]. 
Surface wave antennas have been analyzed on many 
occasions previously but Hirsch [15] has carried out an 
analysis of higher-order surface wave modes supported 
on a dielectric coated conductor. These modes are all 
hybrid modes and most of them are very complicated. 
But although such modes can have low attenuation, on 
suitably proportioned guides, the exploitation is dif- 
ficult because of analytical complexity and lengthy 
computations. 

Two different surface wave excitation problems have 
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been analyzed. One by Wait and Conda [16] concerns 
essentially launching efficiency of a surface wave from a 
magnetic line source placed above a reactive surface. 
The novel part of the investigation lies in the extension 
of the analysis to launching of a surface wave over a 
corrugated elliptic cylinder. It is shown that for certain 
proportions of the cylinder strong resonance charac- 
teristics are obtained. The results, no doubt, can be ex- 
tended to other structures such as radiation over the 
surface of a sphere. Another excitation problem was 
analyzed by Angulo and Chang [17]. This concerns 
excitation of surface waves along two infinite, identical, 
separate and parallel dielectric slabs by a magnetic 
line source. A modal analysis is carried out and an ex- 
pression for the far field is obtained. Launching ef- 
ficiency, radiation loss and radiation pattern are also 
obtained. It is shown that there exists an optimum 
thickness and spacing of the slabs for maximum ef- 
ficiency. A surface wave problem of a different kind 
was analyzed by Williams [18]. He investigates the 
excitation of surface waves by plane electromagnetic 
wave incident on. a wedge with surface impedancy 
boundary conditions. It is shown that surface waves 
will be excited even though the surface impedance 
values on the two sides of the wedge are unequal, 
provided that certain discrete values of the wedge angle 
are excluded. , 

There has been, and still is, a considerable amount of 
interest and discussions on how surface waves should be 
classified. But, in the opinion of the reviewer an im- 
portant aspect is still not receiving adequate attention, 
that is, we are calling by the same name “surface wave” 
physically different things.* On one side there are the 
modal waves (pole waves) of propagation and these 
could be classified on their own. On the other side there 
are various types of waves associated with the radio 
wave propagation over the ground. We do not think 
that it really matters whether we call the first group 
surface waves and call the other group by a different 
name or the other way round. Provided that a distinc- 
tion between the two groups (let it be three groups, 
if the need arises) is made, classification within the 
group is of secondary importance and will follow much 
more easily. After all, biologists, for example, do not call 
different species of the animal kingdom by the same 
name, but find it easier to differentiate between them, 
despite the fact that borderline cases clearly exist. 


Propagation in Nonuniform Waveguides 


There are many variants of nonuniform waveguides 
and these are discussed under three separate headings: 
1) Single mode propagation in nonuniform lines, 2) 
Propagation in slowly varying nonuniform multimode 
waveguides and 3) Propagation in waveguides with 
statistically distributed irregularities. 


* A geod summary of the situation is given in a review p. 
J. R. Wait presented at the XIIIth General Agombloet URSI. 2 
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The foundations of the theory of wave propagation in 
nonuniform waveguides were established by many 
workers in different countries (e.g., Stevenson, Schel- 
kunoff, Katsenelenbaum, Riter, etc.). The subject is 
very complicated, but by making various simplifying 
assumptions it is possible to solve a number of specific 
problems and to go a long way towards forming a co- 
herent physical picture. 

Problems involving only a single mode can be reduced 
to the problem of design of a taper between two trans- 
mission lines of unequal characteristic impedance. It is 
well known that Tchebycheff design leads to minimum 
length for a given reflection coefficient over an infinite 
bandwidth. But, for most applications only a finite 
bandwidth is required and then it is possible to designa 
shorter taper. Solymar shows that there is an infinity of 
possible solutions to tapers of as short a length as de- 
sired but that the shape of the taper would be extremely 
complicated (a violently oscillating function), and 
therefore the design would be impractical [19]. By in- 
troducing the concept of “function complexity” he was 
able to limit the class of functions and give design cri- 
teria of practical significance. The theoretical results are 
interesting in that there exists a Fourier transform rela- 
tion between the shape of the taper and the reflection 
coefficient. 

A problem somewhat related to tapers is that of 
bends. Here again over a given bandwidth a small re- 
flection coefficient is specified and the problem is to de- 
sign a bend as small as possible. The usual procedure is 
to design a binomial bend; but Wray and Hastie show 
that by going to the limit and taking infinitesimal sec- 
tions, a considerable improvement in the design is pos- 
sible [20]. The analogy with taper is striking, and one 
would guess that further improvement should be possi- 
ble using Solymar’s approach. Presumably many other 
gradually varying transitions could be treated in a 
similar manner, e.g., twisted waveguides [21]. 

An entirely different approach to tapers was pro- 
posed by Barlow [22]. The idea is to maintain a com- 
pletely undisturbed field pattern not only inside but 
also outside the guide, irrespective of its cross section. 
The method has been demonstrated by application to 
the H.» and E,,, modes but it is doubtful whether it will 
be possible to extend it to other modes because phys- 
ically it would be difficult to provide the required 
boundary conditions. We can examine the proposal by 
reference to Fig. 1. Here A and B are arbitrary planes 
extending right across the infinite pattern of the Ho, 
(or Eo,) mode. The planes A and B could then be re- 
placed by suitable impedance sheets and the field: out- 
side omitted. Clearly, A and B form, geometrically, a 
taper, but not from an accepted electrical point of view. 
The reason is that the wave is not really “guided” 
since the energy crossing planes A and B is of the same 
order of magnitude as the component of energy flow 
along the planes. The wave pattern does remain un- 
disturbed but the planes 4 and B are not surface im- 
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pedance boundaries and they do not “guide” the wave. 
No doubt, some applications of this proposal will be 
forthcoming and any experimental work in this field 
would be of interest. 


Design of mode transducers is perhaps one of the 

most difficult problems in the theory of nonuniform 
“waveguides and has eluded a satisfactory solution for 
a long time. Such designs as have been produced have 
been based on engineering intuition. Solymar and 
Eaglesfield have recently proposed a new approach to 
the design of mode transducers [23]. In essence the 
method of design is as follows. If a mode ¥4 in wave- 
guide A is to be transformed into a mode y? in wave- 
guide B then we form a wave function ~=gi(z)W4 
+ g2(z)¥2" where gi and ge are slowly varying functions 
of the axial coordinate z with terminal conditions 
gi(0) =go(L) =1 and g.(0)=gi(L) =0, where LZ is the 
length of the transducer. There is, of course, an in- 
finite number of possible solutions and the optimum 
solution is yet to be discovered. For simplicity the 
eigenvalues of A and B are made equal by preceding 
either one or the other waveguide with a simple taper. 
The shape of the transducer is determined by zeros of 
the W function; this part of the design for all but the 
simplest cases is best accomplished by a graphical pro- 
cedure. The authors have designed a number of dif- 
ferent transducers and the simplest example is that of 
the Ho: rectangular to Ho: rectangular [24]. There does 
not seem to be any limitation to the method, and trans- 
ducers so designed are all wide band. It is most in- 
structive to work out a complete design of a transducer 
because it is then that one gains a deeper and clearer 
understanding of modes in nonuniform waveguides. 

It remains to be proven,‘ but it would appear that 
the complete set of pole waves is undergoing a con- 
tinuous transformation as one continuously changes the 
shape of the waveguide. There is a one-to-one corre- 
spondence between modes in one section of the wave- 
guide and another. No new poles can be created by such 
a continuous process but some poles may merge into one 
or form hybrids, while others, which may have been 
presupposed to be simple poles, will split. It is not 
possible to subdivide modes in general waveguides 


4 The description given in this paragraph represents the reviewer’s 
opinion. No rigorous proof in support of it has yet been produced. 
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into E or H modes or in any other way that would be 
generally valid, but one can label the modes just the 
same. During such a process one should not be sur- 
prised to find that a pole representing an Ey; mode in a 
circular waveguide will become transformed by a one- 
to-one correspondence into (say) an He; mode of a 
circular waveguide, but, if the transformation of the 
waveguide is carried out in a different way, a one-to-one 
correspondence with another mode will be established. 
The type of one-to-one correspondence between the 
modes is purely a function of the shape of the trans- 
ducer. 

Transmission properties of waveguides with statis- 
tically distributed irregularities are of particular im- 
portance in connection with long distance communica- 
tions by waveguide. Various analyses have been carried 
out in Great Britain and the United States during the 
last few years, but recently Larsen (in Germany) has 
prepared yet another study [25]. 


WAVEGUIDE COMPONENTS, DEVICES 
AND MEASUREMENTS 


The field of interest here can be divided into theoret- 
ical analysis of components and various configurations 
of obstacles and apertures, properties of ferrites and 
measurements of dielectric properties and microwave 
power. 

Strip-lines continue to be a subject of interest [26], 
[27]. With most strip-line configurations, it is difficult 
to obtain a good correlation between theory and experi- 
ment. This is understandable because the structure, 
from an electrical point of view, is rather complicated. 
In addition, strip-line forms an open structure and, 
therefore, any discontinuities such as matching stub or 
corner cannot be calculated using conventional wave- 
guide methods because of the associated radiation field. 
Some progress, however, can be made using various 
simplifying assumptions, but the results are under- 
standably approximate. Assuming a substantially TEM 
transmission, Lewin was able to assess the properties 
of a post, a corner and a short circuit; he also ex- 
amined the radiation patterns associated with such dis- 
continuities [28]. He concluded that a short circuit is 
appreciably worse, from the point of view of radiation 
loss, than an open circuit. Only a few of his results, 
however, have received experimental confirmation to 
date. 

The theory of waveguide junctions is now well under- 
stood through Schwinger’s Variational Approach. But 
all waveguide discontinuities could be analyzed from 
first principles by matching the complete waveguide 
fields on both sides of the discontinuity. The matrix 
relating the modes is infinite, but for a single mode 
waveguide the elements of the matrix—which, with the 
exception of the dominant mode, represent only 
evanescent modes—is quickly convergent. Sinha used 
this fact and approximated by summing up just the 
first few terms and in this way was able to obtain the 
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equivalent circuit parameters for a waveguide junction 
[29]. His formulas agree very well with experimental 
results. 

A problem of great practical significance and yet 
difficult to analyze for the general case is coupling of 
waveguides through apertures. Some years ago Steven- 
son obtained a rigorous formulation of the problem, 
but his equations are difficult to apply in most cases. 
Bethe’s approach is relatively simple and has been used 
successfully for the design of many waveguide and 
cavity configurations. Bethe’s method, however, is re- 
stricted to small apertures only, and a new approach, to 
handle the: numerous practical cases where large aper- 
tures are used, is clearly needed. Lewin’s approach 
forms a useful extension to Bethe’s work [30]. The 
idea is to use the results obtained in the antenna 
theory (quasi-static) and apply them wa Babinet’s 
principle to apertures in a waveguide wall and obtain 
corrections to Bethe’s aperture polarizabilities. It would 
be difficult to deal, in this way, with arbitrary apertures 
but the author manages to analyze long narrow slots 
and also a case of two waveguides coupled by crossed 
slots. 

The phenomenon of traveling-wave resonance has a 
number of specialized applications. The theory and 
some of its properties were investigated by Twisleton 
isi 

Theory of microwave propagation in ferrites is of 
great importance in the design of numerous microwave 
components such as isolators, gyrators, phase shifters, 
etc. But, because of the tensorial character of ferrite 
permeability, the general theory is crowded with dif- 
ficulties. Although in a number of specific cases, useful 
solutions can be obtained [32]. Sometimes, even what 
could appear, on superficial examination, to be quite a 
manageable problem, seems to present difficulties. 
Thus, for example, the determination of the input 
impedance of a ferrite-filled waveguide becomes a siz- 
able analytical problem to which only partial answers 
have been obtained. For certain range of parameters 
Lewin was able to obtain a correct solution, but outside 
that region the formulas lead to a physically unac- 
ceptable solution (complex impedance) and have to be 
rejected [33]. Recently however, by taking into ac- 
count Bresler’s surface waves, which propagate in the 
gap between the ferrite and the wall, the author was 
able to justify the validity of his solutions. 

Various measurements on ferrites and their inter- 
pretation are yet another prolific source of mathe- 
matical problems. One of the most favored measure- 
ment methods is one in which the sample of ferrite to be 
measured is prepared in the form of a small sphere and 
is inserted inside a cavity [34]. From the measurement 
of the change in the resonance conditions of the cavity 
the ferrite permeabilities are determined. In the past a 
number of perturbational techniques were developed, 
and there have been many discussions concerning the 
relative merits of various methods. Underlying all 
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theoretical formulas there are various assumptions and 
these have been reviewed and classified by Waldron 
isis 

Various measurements on waveguide components, 
power measurements using resistive films and Q-factor 
measurements have also been described [36]-[40]. 
These incorporate various extensions and modifications 
to the existing methods. A somewhat different method 
of measuring permittivity has been proposed by Sinha 
and Brown [41]. Here the sample, which is in the form 
of a thin rod, is inserted through a hole in the end 
wall of a circular cavitiy operated in the Ho, mode. 
The other end wall is formed by a movable plunger. At 
the desired frequency the loss tangent and permittivity 
are deduced, from the cavity Q factor and the change 
in the plunger position, as a function of the sample 
insertion. 

In the millimetric band, yet another method of meas- 
uring dielectric properties becomes feasible; Caicoya 
refers to the method as the “High-Order Mode Inter- 
ferometer” [42]. In essence, the method consists of 
transmitting a parallel beam of millimetric radiation at 
an angle to a parallel slab of material (under test) 
placed between metallic plates. From the transmission 
loss, as a function of angle of incidence and the geome- 
try of the sample, the permittivity and loss tangent are 
deduced. 

This is an illustration of quasi-optical methods in the 
microwave band. With the growing interest in the upper 
millimetric band, such methods are likely to be applied 
to the design of microwave components and, perhaps, of 
complete microwave systems [43]. General mechanical 
design and manufacturing techniques have been re- 
viewed by Harvey [44], [45]. 


SOLID-STATE DEVICES 


The fields of interest here are transistors, diodes in 
general and tunnel diodes in particular, parametric 
amplifiers and masers. 

Transistors are really beginning to reach microwave 
frequencies [46] but it is doubtful whether transistors 
of the types available to date will be a commercial 
proposition. A real success of transistors at microwave 
frequencies hinges on the development of new types. 

It is probably true to say that none of the semicon- 
ductor devices has undergone such a rapid development 
and wide acceptance as has the tunnel diode [47]; yet, 
for microwave applications the prospects for the tunnel 
diode are not too bright. Currently available tunnel 
diodes are capable of delivering only tens of micro- 
watts of power at the lower microwave frequencies, and 
although research is taking place in a number of es- 
tablishments, there does not seem to be much chance of 
inducing a tunnel diode to deliver more than about one 
milliwatt of power at the lower microwave frequencies 
or tens of microwatts at the upper microwave fre- 
quencies. Clearly, junctions with much higher current 
densities are needed. This suggests the use of high 
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mobility compounds such as gallium arsenide, indium 
antimonide, etc. Using wafer construction, a number of 
gallium arsenide diodes have been produced [48]. Al- 
though the technology of such diodes is more involved 
' than that of the germanium diodes, these diodes are 
preferred for microwave applications. 

Microwave parametric amplifiers and masers form a 
very new field. Parametric amplifiers entered the micro- 
wave field in 1957 (in the United States) through the 
work of Suhl and Hines but its rapid progress was made 
possible as a result of Uhlir’s work on parametric 
diodes. Parallel with the diode work, ferromagnetic 
amplifiers and parametric beam tubes were developed. 
The age of maser devices is about the same as that of 
parametric amplifiers, but the chief impetus to the 
maser work was received from Bloembergen’s (1958, 
United States) proposal of a three-level solid-state 
maser. 

The interest in parametric amplifiers and masers 
quickly spread to other countries, and in Great Britain 
there are a number of establishments working in this 
field. Parametric amplifiers were discussed in a number 
of review articles (e.g., [49], [50]). Clarricoats analyzed 
the traveling-wave ferromagnetic amplifier [51] and 
Thompson reported on an investigation of spin wave 
excitation in YIG [52]. Various other aspects of para- 
metric amplification were also studied, such as gain, 
bandwidth and noise figure [53], circuit aspects [54], 
electron beam amplifiers [55] and measurements on 
broad-band amplifiers [56]. Cullen was one of the 
first to analyze the traveling-wave parametric amplifier 
[57]. He was able to take into account a number of 
factors influencing the performance of amplifiers, such 
as ohmic losses and saturation phenomenon [58], and 
deduced the consequential maximum gain. 

The commercial success of parametric devices at 
upper microwave frequencies—apart from, possibly, 
parametric beam devices—depends largely on success- 
ful development of parametric diodes with cutoff fre- 
quencies in the region of 500 Gc. Thus, for example, 
efficient harmonic generation in the millimetric band 
is not possible with presently available diodes because 
of the proximity of cutoff. On the other hand, it is in 
the upper millimetric band that microwave tubes are 
also inefficient and, therefore, diodes as harmonic 
sources might yet succeed. 

A different approach to harmonic generation has been 
reported by Froome. He obtained good harmonic yield 
using the nonlinear voltage current and relationship of a 
mercury arc [59]. Subsequently a somewhat different 

structure was investigated. Using tungsten anode and 
calcium cathode in argon atmosphere, Froome was suc- 
cessful in obtaining frequencies in the region of 300 Gc 
(most recent reports state that nearly 600 Gc were 
produced) from a source of 35 Gc [60]. The efficiency 
at present is low (it still compares favorably with any 
other source available at present), but further work is 


proceeding. 
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One of the many reasons why interest in parametric 
amplifiers is so great and is sustained is because of the 
very low noise figure realizable with such devices. In 
this respect, however, masers (despite their greater 
complexity) have an overwhelming advantage; yet the 
choice depends on the particular application envisaged 
[61]. 

General review articles concerning masers have ap- 
peared (e.g., [62], [63]). The activity can be divided 
into 1) investigation of materials and modes of opera- 
tion, 2) two level devices, 3) three level devices and 
applications, 4) ammonia maser oscillators, 5) traveling 
wave masers, 6) optical and infrared masers. 

Although from a scientific point of view masers are a 
great success, any great commercial success is intimately 
related to the development of more suitable materials— 
materials with smaller losses and more favorable re- 
laxation times, materials for operation at elevated tem- 
peratures (say 60°K), materials suitable for operation 
in millimetric band, etc. 

Solid-state masers using materials with suitably di- 
luted paramagnetic ions are of chief interest. Of these 
materials, those containing Cr ions dissolved in 
Al,O3, KsFe(CN)s, TiOs, etc., are preferred. A number 
of investigators have carried out studies of relaxation 
phenomena (e.g. [64], [65]), while rutile was investi- 
gated by Thorp for operation in the millimetric band 
[66]. Two-level masers are also of interest in the study 
of properties of paramagnetic ions in a host crystal 
lattice [67]. 

A number of engineered versions of masers, mainly 
for operation at lower microwave frequencies, have been 
described [68]-|70]. Ditchfield reported on an X-band 
maser operating at 77°K; this is a significant step for- 
ward, but, of course, at the expense of efficiency [71]. 
He also reported one study of variation of relaxation 
times with temperature. 

Cavity masers are substantially (but not necessarily 
so) narrow-band amplifiers of negative resistance type. 
For many applications this is a limitation. For this 
reason there is a wide interest in traveling-wave masers 
because such amplifiers are reputed to be free from such 
disadvantages [62], [72]. The preferred construction 
employs a comb slow-wave structure with ruby as the 
active material. 

The power output of ammonia masers is rather small 
and the bandwidth very narrow. For this reason such 
devices are not likely to be used as amplifiers. As an 
oscillator, the ammonia maser is characterized by a 
high-frequency stability. The prospect of ammonic 
maser oscillator as a frequency standard is, therefore, 
good and at least in one research establishment work is 
proceeding in that direction [73]. The ultimate fre- 
quency stability may probably be of the order of a few 
parts in 10". 

There are many more fields of interest to those en- 
gaged on microwaves: experiments with plasma, cyclo- 
tron resonance, etc. The published information is 
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scarce however and largely confined to review articles. 
Infrared and optical masers, however, ought to be men- 
tioned [74] because of the large impact such devices 
are likely to have on our future developments. 


UTILIZATION 


Apart from a number of industrial applications for 
a variety of uses [3] and applications in defense proj- 
ects, microwave links remain the chief consumer of 
microwave equipment. 

There is, however, a widespread interest:in communi- 
cation with and wa satellites; this is evidenced by the 
large number of lectures and public talks on the subject. 
The chief effort in this direction seems to be made in 
the United States; understandably so, because means 
for putting satellites in space are not, at present, at the 
disposal of European countries, except, of course, 
Russia. 

In Britain the activity in satellite communication is 
largely confined to considerations of electronic equip- 
ment, investigations of upper atmopshere and theoret- 
ical assessments of future prospects; the published in- 
formation, however, is very scarce. 

There is active interest in long-distance communica- 
tion by waveguide. In fact it transpires that Great 
Britain, Japan, Russia and the United States are the 
four countries where research in this field is taking place 
on a large scale. But, there is interest in a number of 
other countries, notably France and Germany. 

Progress in the field of long-distance communication 
by waveguide can be gathered from an article by 
Karbowiak and Solymar [75]. It appears that for ef- 
ficient operation a waveguide communication system 
must utilize frequencies corresponding to the upper 
millimetric band. The authors also point out that the 
use of special waveguides—either of the dielectric 
coated variety or a helical construction—is, from a 
practical point of view, a necessity. It is also stressed 
that modulation methods such as PCM, where signals 
can be regenerated at repeaters spaced about 20 miles 
is also a necessity, but that PCM is not the only 
possible choice and, furthermore, for short-haul ap- 
plications a variety of modulation methods could be 
used with success. 

An incidental application of helical waveguides to 
the construction of high-Q wavemeters has been sug- 
gested by Barlow, et al. [76]. In this way the authors 
propose to overcome harmful interference from other 
modes and thereby improve the wavemeter perform- 
ance. 


FUTURE TRENDS AND CONCLUSIONS 


The field of microwaves is now so immense that 
one approaches the task of a critical review with a 
certain amount of diffidence. It is difficult, indeed im- 
possible, to read and listen to everything that concerns 
microwaves and yet give a critical account of the work 
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done even though the work concerns a period of time 
as short as one year. A great deal of preliminary selec- 
tion has been carried out and it is therefore possible 
that someone’s important contribution has been 
overlooked. The material has been selected, and the 
space devoted to the review adjusted, in the following 
order: new or controversial matters, work of funda- 
mental character or one that is likely to have a signif- 
icant impact on future developments, miscellaneous 
fascinating discoveries and modifications. Fields which 
have not been included in the review (with some excep- 
tions) are: antennas, propagation and microwave tubes. 
Whereas most of the space was devoted to the theory of 
microwaves and waveguides, various component meas- 
urements and solid-state devices occupy the remainder 
of the review. Utilization has also been briefly sum- 
marized. 

In assessing the future trends, it is important to 
observe that most of the progress achieved in the micro- 
wave field is a direct result of research work. Thus 
future activity and trends will be, in a way, a reflection 
of the attitude towards research prevailing throughout 
the country; the more we are research-minded the 
greater are likely to be the fruits. It is the attitude of 
mind that determines the future. Admittedly, research 
is an expensive business, but some years ago the atti- 
tude towards research was such as to justify the saying 
that research is yet another form of hobby. 

While not disparaging the spirit behind such an atti- 
tude, the reviewer welcomes the change which has taken 
place during recent years. For it is only by recognizing 
that research is an essential part of industrial, and in- 
deed of national activity, that a rapid progress in new 
fields can possibly be made. 

Which fields should be exploited first in order that 
the harvest be worthwhile? This is the big question 
asked over and over again by many. The simple answer 
is that in the majority of cases there can be no answer 
on the available evidence. In fact, it is wrong to ask 
such a question; for, what is worthwhile? 

It is evident from the large volume of published ar- 
ticles and lectures on research management that there 
is a growing concern (see, for example, [77]) about 
most aspects of research. It would appear, however, 
that the surest way to success in research, if that is 
the aim, is faith. Faith of one member of a research 
team in another, faith in our profession, and faith in our 
aims. There is, really, no other recipe. While nowadays, 
a lot of fine work is still being done by individual people, 
research is yet another form of social intercourse, but of 
a very specialized nature. 

It is then the research scientists and engineers them- 
selves that have the right qualifications, training and 
education to guide and influence the future of research. 


° In such cases the reviewer would be grateful if his attention 
were drawn to the fact, so that the contribution could be taken into 
account at a later date. 
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Yet, so few of them are in sufficiently senior positions 
to undertake this duty. In this connection Lord Wol- 
ton’s words spoken in the House of Lords debate on 
November 21, 1956, may be recalled without comment: 


§ 


If industry wants to get first-class scientists and 
engineers into its service, it has to make perfectly 
clear to the scientist what can happen to him in the 
future. The scientist does not want to stay in the 
laboratory all his days, and just to give advice. He 
wants to get into a position in which, equally with 
the finance director or sales director, he will be able 
in the board room to influence the policy of the 
business. 


Having now established that any future forecast of 
research progress is next to impossible, we can now 


_summarize the future trends in the microwave field as 


follows. 
There is an indication that in the next few years the 


fundamental standards of physics will be reviewed; 


microwaves will play an important part in it. 

There is already a greater understanding of theory of 
electromagnetic fields, and an even deeper understand- 
ing of guided waves in waveguides of arbitrary struc- 
ture should follow. More and more powerful mathe- 
matical techniques (such as the Wiener-Hopf technique) 
are being applied to the solution of waveguide prob- 
lems. We should, therefore, witness a solution to many 
outstanding problems. 

Microwave techniques are also developing, but there 
is a general tendency towards higher and higher fre- 
quencies. New developments in the upper millimetric 
band are anticipated and some of the engineering prob- 
lems connected with the exploitation of the submilli- 
metric spectrum should be overcome. 

Efficient generation of coherent radiation in the milli- 
metric and submillimetric spectra is an outstanding 
problem whose solution will probably lie in the invention 
of suitable solid-state devices. We may even witness the 
closing of the present gap between submillimetric and 
infrared spectra. 

In utilization, the application of microwaves to com- 
munication systems will dominate our activity (except- 
ing defense projects). In this connection, satellite com- 
munication seems to have captured the imagination of 
many, chiefly because satellites will be launched and 
rocket research will proceed irrespectively of any other 


- uses, and also because the problems are fascinating and 


there is a great deal to learn. The commerical success, 
however, is at present doubtful, but in the infrared or 
the visible spectrum the chances are perhaps better. 

Long distance communication by waveguide will be- 
come a reality in a not-too-distant future, particularly 
since the engineering in the upper millimetric band has 
reached, now, a high degree of development. 

While on the subject of communication systems, it 
should be observed that although telephone traffic con- 
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tinues to increase, other uses (such as television, data 
transmission, etc.) which the communication systems 
serve, increase even faster. It may be that in the years 
to come telephone traffic will form but an insignificantly 
small proportion of the activity in communication sys- 
tems and we will have established even better ways, not 
yet conceived, of exchanging our ideas. 
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A New Type of Circular Polarizer Using 
Crossed Dipoles* 


M. F. BOLSTER, MEMBER, IRE 


Summary—A method of obtaining a circularly-polarized wave by 
use of two orthogonal dipoles driven in parallel by a common trans- 
mission line is shown. The lengths of the dipoles are so chosen that 
the real part of their input admittances are equal and the angle of the 
input admittances differ by 90°. When these two conditions are met 
the resulting radiated wave in a normal direction will be circularly 
polarized. 

The method is applicable both to a circularly-polarized radiating 
antenna and to the problem of producing a circularly-polarized wave 
of the TE;, mode in a round waveguide. For the first case, an analysis 
and a method of design are shown, and for the second case an ex- 
perimentally developed example is given. The second case employs 
monopoles rather than dipoles for convenience in energizing from a 
coaxial line. 


INTRODUCTION 


polarization using crossed dipoles.! Previous meth- 

ods have used identical resonant (approximately a 
half wavelength long) dipoles in each plane, and the 
proper power and phase relationships in the dipoles 
were obtained by use of an input phasing and matching 
network so as to drive the two dipoles from a common 
generator. 

The method presented here uses two unequal length 
dipoles connected in parallel to a common generator and 
requires no matching network. The required power and 
phase relationships are obtained by proper choice of 
the two dipole lengths. For the first case (dipoles in 
free space), the lengths are calculable; for the second 
case (monopoles in waveguide), the lengths must be 
determined experimentally. 


| T has been common practice to produce circular 


CONDITION FOR CIRCULAR POLARIZATION 


In Fig. 1 are shown two orthogonal cylindrical 
dipoles fed in parallel and a representative equivalent 
circuit. The mutual interaction between the two di- 
poles should be negligible. as long as their lengths are 
much greater than their radii and each dipole is placed 
in the neutral (or zero) potential plane of the other di- 
pole (as is shown on Fig. 1). 

The calculated input admittance of a single dipole 
as a function of its length is shown in Fig. 2 for lengths 
in the neighborhood of a half wavelength. These values 
of admittance are second-order calculated values as 


* Received by the PGMTT, April 16, 1961. | 
+ General Engineering Lab., General Electric Co., Schenectady, 


i GFH. Brown, “The turnstile antenna,” Electronics, vol. 9, 
pp. 14-17; April, 1936. 


SUSCEPTANCE (B) MHOS (I0)~> 


8 10 12 |4 16 


4 
CONDUCTANCE (G) MHOS (10)7> 


Fig. 2—Dipole input admittance. 


given by King? for 


he ft 
ae = 75.2. 
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In general, the two dipoles will produce an elliptically- 
polarized wave. In order to have a circularly-polarized 
wave, there must be equal power input to each dipole 
and 90° phase difference between the input currents. 
This is equivalent to having G,=G, and arg Y,=arg 
Y,+90° (where Y.=G,+JB,Y,=G,+JB,). This con- 
dition can be determined from an accurate plot of 
dipole admittance by a graphical solution as shown in 
ea as 


2R. W. P. King, “The Theory of Linear Antennas,” The Harvard 
University Press, Cambridge, Mass., p. 172; 1956. 
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Fig. 3—Real and imaginary parts of dipole admittance. 


If the admittance locus were extended to show longer 
dipole lengths, it would be evident that there are many 
other combinations of (longer) lengths that would also 
produce circular polarization. 


VARIATIONS WITH FREQUENCY 


When the frequency of the driving signal is changed 
the radiated wave will become elliptically rather than 
circularly polarized. To determine how the input ad- 
mittance and ellipticity varies as a function of fre- 
quency, it is convenient to replot the real and imaginary 
parts of the admittance in Fig. 2 as functions of 
Bh=(27r/d)h as shown in Fig. 3. 

From Fig. 3 the two dipole admittances V,=G,+jB; 
and Y,=G,+jB, can be obtained graphically as func- 
tions of frequency by simply reading values of 
Bhz=Bohz(f/fo) and Bhy=Boh,(f/fo), where fo is fre- 
quency at which circular polarization occurs. Perform- 
ing this operation over a frequency range of f=fyp+15 
per cent the values of G and B shown in Fig. 4 are ob- 
tained. 

The total input admittance Y;=G,;+jB;=Y,+Y, is 
shown as a function of frequency in Fig. 5. 


POLARIZATION ELLIPSE OF RADIATED WAVE 


From the data of Fig. 4 the polarization ellipse of the 
radiated wave can be calculated. The power in the x or 
y component waves will be proportional to G, or G,, and 


the amplitude of the component electric fields E, or Ey 
will be proportional to WG, or /G,. Other nomencla- 
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Fig. 5—Input admittance of dipoles connected in parallel. 


ture concerning the polarization ellipse shown in Fig. 6 
is defined as follows: 


E,=electric field component in x (dipole) direction 
E,=electric field component in y (dipole) direction 
E,’ = electric field component in x’ direction 
E,'=electric field component in y’ direction 
7 =polar angle between xy and x’y’ frames 
T)=polar angle between dipoles and principal axis of 
polarization ellipse 
A,=</G,=amplitude of electric field component in 
x (dipole) direction 
A,=</G,=amplitude of electric field component in y 
(dipole) direction 
R=A,/A, 


A,sin6,cos7 + A, sin 6, sinr 
62 =stan 2 : = : ee (5) 
A,cos6,cost + A, cos 6, sint 


The angle tT) between the x dipole and the principal 
axis of the polarization ellipse is obtained by dif- 
ferentiating A,’ with respect to x, setting the result equal 
to zero, and solving for 7 =79. Carrying out these opera- 
tions and simplifying the result yields 


To = 1/2 tan“! sao ) (6) 
1 — R? 
where 
R=A,/Az and 6= 4, — 6. (7) 


This is the same formula which Born*® obtained by a 
different method. 
The axial ratio of the polarization ellipse is defined as 


A,’ 
Bee a : 
oe SEN ae (8) 


where A,’ is given by (4), and A,’ is obtained in similar 
manner from 
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6,=phase of electric field component in x (dipole) 
direction 
6,=phase of electric field component in y (dipole) 
direction. 
The field in the x’ direction will be 
E,’ = E, cost + E, sin r, e5) 
where 
E, = A,cos (wi + 6,) and E, = A, cos (wt + 6,). (2) 
This can be written 
E,’ = A,’ cos (wt + 4,’), (3) 
where Fig. 6—Polarization ellipse. 
A, = (A; cos 6, cos r + A, cos 6, sin)? + (A, sin 6, cosr + A, Sin dy sin 7)?, (4) 
and The axial ratio is obtained from this equation when 


2R cos 6 


— 

1/2 tan rage 

is substituted for 79. This calculation was performed 

using the data of Fig. 4 and gave the axial ratio (in db) 

and orientation of the polarization ellipse shown in Fig. 
7 over the 30 per cent band. 

If the signs of cos 6 and 1—? in (6) are considered, 
then the angle 7) would have to be in the second quad- 
rant for values of f>jfo; however, due to the symmetry 
of the polarization ellipse, plotting 7) in the fourth 
quadrant as shown will cause no error in describing the 
ellipse. 

The shape of the curve for tT) shown in Fig. 7 in the 
region of f=f) would probably not be duplicated in 
measurements on an actual antenna. It is unlikely that 
there would be any cusp at this point. The cusp prob- 
ably appears in the calculated values because small 
deviations in the input values of RK and 6 in the neighbor- 
hood of f=fo result in relatively large changes in the 
calculated values of to due to the indeterminate form 
of (6). 


E,’ = — E,sint + E, cost = Ay’ cos (wt + 4,’) (9) 
A,! = V/(— A; cos 6, sint + A, cos 6, cos 7)? + (— Azsin6é,sinz + A,sin 6, cos 7)?. (10) 
; ¢ ‘ . cos 6 0 
After some simplification, one obtains lim | | eG) (12) 
1+ R Ro Ll — Kk? 0 
a (11) 590° 


| Ry? =—— to . 5 
cos? 79+ R? sin? ro +2R cos 79 Sin To CoS 6 


8 Max Born, “Optik,” Verlag Julius Springer, Berlin, Germany, 
p. 23; 1933. 


The order of the zero of cos 6 is greater than that of 
1—R? when R and 6 are considered as functions of fre- 
quency. 
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Fig. 7—Axial ratio and tilt angle of polarization ellipse. 


APPLICATION TO WAVEGUIDE 


The same technique can also be used to produce cir- 
cular polarization in round waveguide, however, the 
characteristics are not as readily calculable. Consider a 
pair of crossed monopoles mounted at 90° to each other 
and connected together at one end. The two could then 
be fed in parallel from the center conductor of a co- 
axial line. If the assembly is next placed inside a 
round waveguide it is possible to obtain circular polari- 
zation in the waveguide by experimental adjustment of 
the monopole parameters. 

An experimentally-adjusted circular polarizer of this 
type using a dielectric-filled waveguide is shown in Fig. 
8. This polarizer was designed to operate over a 10 per 
cent frequency band from 3000 to 3300 Mc. The mono- 
pole elements are arranged to rotate in a circular slot 
in the dielectric (e=2.54), and the waveguide diameter 
(1.666 inch) was chosen so as to propagate only the 
TE mode over the desired band. 

For the polarizer and coaxial matching transformer 
shown in Fig. 8, the measured input VSWR and axial 
ratio of the wave propagated in the waveguide are 
shown in Fig. 9 for a +10 per cent frequency band 
centered at fo =3.15 Ge. 

It may be noted that there is approximate agreement 
in the rate of change of ellipticity between the simple 
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Fig. 8—Circular polarizer in dielectric-filled round waveguide 
with input matching transformer. 
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Fig. 9—VSWR and axial ratio of polarizer in dielectric-filled 
round waveguide with input matching transformer. 


dipoles and the more complex waveguide polarizer of 
Fig. 8. For a bandwidth of fo+10 per cent, the dipole 
ellipticity is about 9 db (Fig. 7), while it is about 5 or 6 
db for the waveguide polarizer (Fig. 9). 


CONCLUSIONS 


It is possible to obtain circular polarization by proper 
choice of the dimensions of the radiating elements them- 
selves, rather than by use of power splitting and phasing 
networks. The principle is applicable to any type of ele- 
ments which have input admittance characteristics 
similar to the dipole of Fig. 2, and it may be accom- 
plished either in free space or inside waveguide. 

The free-space type of polarizer would be useful as a 
circularly-polarized antenna. The fact that no input- 
phasing network is necessary is a definite advantage. 
The waveguide type could be used as a waveguide feed 
for a horn antenna or as a phase shifter by mechanically 
rotating the monopoles about ‘the axis of the coaxial 
line. The small moment of inertia would allow rapid 
rotation with better angular stability and smaller 
driving torque than previous methods (rotating helix or 
rotating dielectric plate). 
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Rounded Corners in Microwave High-Power Filters 
and Other Components’ 


SEYMOUR B. COHN}, FELLOow, IRE 


Summary—Microwave high-power filters must be operated with 
internal air pressures of at least one atmosphere, or with a good 
vacuum. Pressures between these extremes result in reduced power- 
handling ability. The breakdown processes for both high air pressure 
and vacuum are discussed, and it is made clear that any sharp corner 
on which the electric field would concentrate must be rounded if 
high-power operation is to be achieved. For good results in vacuum 
operation, the surfaces must be especially smooth and free of con- 
tamination, while in high-pressure operation, minor irregularities are 


_ less important. 


Various high-power filter configurations of importance are de- 
scribed, and the structural corners at which electric-field concentra- 
tions occur are pointed out. A number of simplified geometries are 
then shown that can represent the essential portions of the practical 
structures with sufficient accuracy for ordinary purposes. Formulas 
and graphs for these simplified geometries are presented that give 
the ratio of the maximum electric field strength on the boundary to a 
uniform reference field strength at a point sufficiently removed from 
the corner. In some cases, the boundary curve is an approximation to 
a circular arc, while in other cases a boundary shape is derived such 
that the electric field strength along the curve is constant. These 
constant-field-strength boundaries are optimum shapes from the 
standpoint of power-handling ability. 


I. INTRODUCTION 


HE maximum allowable power flow in high-power 
[ites or other components is limited by ionization 

breakdown in regions of high electric-field con- 
centration. The critical value of electric field above 
which breakdown occurs depends upon a number of 
interdependent factors: 1) the composition and pressure 
of the air or other gas filling the device; 2) the signal 
frequency; 3) the size and shape of the region over which 
the electric field approaches its maximum concentra- 
tion; 4) the presence of nearby conducting surfaces, 
their shapes, and their spacings; and 5) the pulse length, 
shape, and repetition frequency. 

The breakdown process as a function of the above 
parameters has been analyzed by Gould.! His work 
shows that for usual sizes of uniform waveguides and 
coaxial lines operated near standard-atmosphere air 
pressure, the ratio of the peak RF electric field strength 
at breakdown in volts per centimeter to the pressure 
in atmospheres is approximately independent of fac- 
tors 2)—5), and is near to the previously accepted value 


* Received by the PGMTT, March, 7, 1961; revised manuscript 
received, May 25, 1961. This work was performed at Stanford 
Res. Inst., Menlo Park, Calif., under the support of the Rome Air 
Dev. Ctr., under Contract No. AF 30(602)-1998. 

+ Rantec Corp., Calabasas, Calif. ; ‘ 

11. Gould, “Handbook on Breakdown of Air in Waveguide 
Systems,” Microwave Associates, Inc., Burlington, Mass. Rept. on 
Contract Nobsr 63295; April, 1956. F ‘ 

Also, L. Gould and L. W. Roberts, “Breakdown of air at micro- 
wave frequencies,” J. Appl. Phys., vol. 27, pp. 1162-1170; April, 1956. 


of 29,000 volts/em/unit pressure in atmospheres.? 
Gould has not considered more complex geometries, but 
his data for coaxial lines indicate that this figure would 
hold approximately for rounded-corner radii as small 
as 1 mm at one atmosphere pressure, with the break- 
down field strength being higher for smaller radii or 
smaller pressure. (For constant error, radius times pres- 
sure is a Constant.) As the pressure is decreased below 
one atmosphere, factors 2)—5) begin to have much 
stronger effects. However, since a high-power filter 
would usually be operated with at least one atmosphere 
of pressure, it appears a reasonably good approximation 
to treat the value of 29,000 volts/cm/atmosphere as a 
constant. Fortunately, this is a conservative approxi- 
mation, since the effect of factors 2)—5) is to raise rather 
than lower the breakdown field strength. 

At a pressure of the order of one millimeter of mer- 
cury, the breakdown field strength as a function of 
pressure passes through a minimum. At lower pressures, 
the breakdown field strength increases rapidly. At pres- 
sures lower than about 10-* mm Hg, corresponding to a 
good vacuum, ionization of the remaining air molecules 
is no longer important and another mechanism of break- 
down becomes limiting. Much less is known about RF 
breakdown in vacuum than at high pressures. However, 
current evidence indicates that the RF breakdown 
mechanism in vacuum is field emission of electrons 
from conducting surfaces. Field emission requires 
field strengths of the order of megavolts per cm.* Such 
field strengths may develop at minute irregularities, 
thus initiating breakdown, even when the macroscopic 
field strength is much lower. Because of this it is ex- 
tremely important that the surface be smocth and free 
of contamination. 

It has thus been established that the attainment of 
prespecified values of electric field strength will result in 
breakdown of high-power microwave filters at high 
pressures. In vacuo, it is not certain whether a definite 
breakdown field strength can be prespecified, but never- 
theless the breakdown power can certainly be in- 
creased by minimizing electric-field concentrations. 
Thus either in high pressure or in vacuo, any edges at 
which the electric field may concentrate should be 
rounded. As an aid in determining the degree of round- 


2H. A. Wheeler, “Nomogram for Some Limitations on High- 
Frequency Voltage Breakdown in Air,” Wheeler Labs., Inc., Mono- 
graph No. 17; May, 1953. 

3 W. S. Boyle, P. Kisliuk, and L. H. Germer, “Electrical break- 
down in high vacuum,” J. Appl. Phys., vol. 26, pp. 720-725; June, 
1955. 
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ing necessary, a number of rounded-corner geometries 
important in microwave-filter structures have been 
analyzed by conformal-transformation methods, and the 
results are given in this paper. For each geometry con- 
sidered, the maximum field is related to a reference field 
at a sufficient distance from the rounded edge that the 
reference field will be essentially uniform. This informa- 
tion will enable filter designers to relate the maximum 
field strength occurring in a given filter to the power 
passing through it. The breakdown value of power may 
then be determined in the case of high air pressure by 
setting this maximum field strength equal to approxi- 
mately 29,000 volts/em multiplied by the pressure in 
atmospheres. A safety factor should of course be ap- 
plied to the calculated breakdown power in order to 
arrive at a safe power rating for the filter. 

Another phenomenon which can, under special cir- 
cumstances, reduce the power transmitted through an 
evacuated high-power filter is called multipactor.4 Multi- 
pactor is a resonant secondary-emission process in an 
evacuated region that occurs when an electron under 
the action of an RF electric field has a time of transit 
between opposite walls of the region that equals ap- 
proximately one-half the period of an RF cycle. Limited 
evidence available to this author indicates that in high- 
power filters utilizing low-Q elements, multipactor 
would ordinarily absorb no more than a few watts of 
power, and hence would not greatly affect performance. 
If high-Q elements are used however, the loss may be 
large. 


Il. TypicaL HicgH-PowEr FILTER CONFIGURATIONS 


Fig. 1 shows simplified sketches of a number of basic 
high-power filter structures in which electric-field con- 
centration occurs at corners. At all outside corners in 
the £ plane, such as those marked A, the field strength 
would approach infinity if the radius of the corner 
were made to approach zero. At all inside sharp corners, 
such as those marked B, the field strength is zero and 
is hence no problem. Fig. 1(a) is typical of the corru- 
gated, or varying-impedance filter.6* Fig. 1(b) is a 
modification of Fig. 1(a), in which the angle of the 
corner is changed to decrease the field concentration.? 
Figs. 1(c) and 1(d) are leaky wall filters, in which slots 
in the broad walls of the main waveguide open into 


*W. G. Abraham, “Interactions of Electrons and Fields of Cavity 
Resonators,” Ph.D. dissertation, Dept. of Elec. Engrg., Stanford 
Univ., Stanford, Calif.; 1950. 

®S. B. Cohn, “Design relations for the wide-band waveguide 
filter,” Proc, IRE, vol. 38, pp. 799-803; July, 1950. 

_°B. M. Schiffman and S. B. Cohn, “Wide-Stop-Band Waveguide 
Filters,” presented at the 1959 PGMTT Natl. Symp., Harvard Univ., 
Cambridge, Mass.; June 1, 1959. 

Also see S. B. Cohn, et al., “Design Criteria for Microwave Filters 
and Coupling Structures,” Stanford Res. Inst., Menlo Park, Callif., 
Tech. Rept. No. 2, Contract No. DA36-039-SC-74862; June, 1958. 

7 J. H. Vogelman, “High-power microwave filters,” IRE TRANS. 
ON Microwave THEORY AND TECHNIQES, vol. MTT-6, pp. 429-439; 
October, 1958. 
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SECTION C-—C 


Fig. 1—Typical rounded-corner configurations 
in high-power filters. 


secondary waveguides of higher cutoff frequency.*? 
Fig. 1(e) is a coupled-cavity filter with rounded edges at 
the apertures.§® 

Rigorous solutions of the structures in Fig. 1 would 
be virtually impossible to achieve. However, the solu- 
tions for a number of simpler cases shown in Fig. 2 have 
been obtained that may be applied to the actual filter 
structures with accuracy sufficient for practical pur- 
poses. The two most important assumptions are: 1) 
the rounded corner geometry may be considered to be 
composed of infinite cylindrical surfaces, so that the 
solution for a single two-dimensional cross section is 
sufficient, and 2) the essential portions of the cross 
section are small compared to a wavelength, so that the 
field distribution in those regions may be considered 
to be very nearly a static field distribution. The first 
assumption applies exactly to cases like Fig. 1(a) and 
1(b), and with fair accuracy to cases like Fig. 1(c)—1(e). 
The second assumption will usually apply with good 
accuracy in the pass band of these filters, even though 


8S. B. Cohn, “Design considerations for high-power microwave 
filters,” IRE TRANS. ON MICROWAVE THEORY AND TECHNIQUES, vol. 
MTT-7, pp. 149-153; January, 1959. 

*V. Met, “Absorptive filters for microwave harmonic power,” 
Proc. IRE, vol. 47, pp. 1762-1769; October, 1959. 

0 'V. G, Price, R. H. Stone, and V. Met, “Harmonic suppression 
by leaky-wall waveguide filter,” 1959 IRE WESCON ConvENTION 
REcoRD, pt. 1, pp. 112-118. 
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(a) 


E, 


Fig. 2—Two-dimensional geometries considered in this paper. 


the over-all waveguide or cavity dimensions are large 
parts of a wavelength. 

The two-dimensional geometries considered in this 
paper are shown in Fig. 2. The angular designation ap- 
plied to each corner in Fig. 2 indicates the change in di- 
rection of a tangent line as it is moved around the 
corner. The array of 180-degree rounded corners of 
Fig. 2(a) has most direct application to the leaky-wall 
filters of Fig. 1(c) and 1(d). The individual conductors 
of the array are assumed to be at the same potential. 
Their center-to-center spacings are assumed to be sulf- 
ficiently smaller than the spacing to any other conduct- 
ing surface that the field well below the array of corners 
may be considered uniform, and equal to a reference 
value Ey. The 45-degree rounded corner of Fig. 2(b) 
applies to the varying-impedance filter of Fig. 1(b). By 
image arguments, the unsymmetrical and symmetrical 
forms of the 45-degree corner may, of course, be shown 
to have identical solutions. In this case, the reference 
field is the uniform field E, in the parallel-plane region 
at a point sufficiently removed from the corner. The 
90-degree corner near an electric wall, Fig. 2(c), applies 
to varying-impedance filter structures when the spacing 
between the broad walls of the main waveguide is small 
compared to the width of the slots. The 90-degree 
rounded corner near a magnetic wall, Fig. 2(d), may be 
applied to slots in the broad wall of a waveguide when 
the RF voltage across the slot is relatively small. This 
would usually be true in the pass band for the leaky- 
wall filters of Figs. 1(c) and 1(d), but not usually in the 
pass band of the varying-impedance filter of Fig. 1(a). 
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In that case, when the dimensions are about as shown 
in the figure, the field in the vicinity of the slot opening 
may be obtained to a good approximation by super- 
imposing the field distributions of Figs. 2(c) and 2(d). 
The coupled-cavity filter of Fig. 1(e) is an example of 
a more complex structure to which the rounded-corner 
data of this paper may be applied. The maximum elec- 
tric field strength will occur in the central longitudinal 
E plane of the filter. In the absence of the coupling 
apertures, the field would be greatest at the centers of 
the cavities, and zero along their vertical sides. In the 
presence of the coupling apertures, however, there will 
be electric field concentrations along the top and bot- 
tom edges of the apertures that could exceed the field 
strength at the cavity centers unless the edges are suf- 
ficiently rounded. The maximum field strength at the 
input and output couplings of Fig. 1(e) may be deter- 
mined from the solution for Fig. 2(c), and that at the 


central coupling from the solution for Fig. 2(a), with 
t/1&1. 


III. ARRAY OF 180-DEGREE ROUNDED CORNERS 


In the case of an array of 180-degree rounded corners, 
a method of analysis due to Wheeler! has made possible 
the derivation of the curved-boundary shape upon 
which the £ field is constant (see Part A of Appendix). 
Figs. 3 and 4 show to scale the resulting boundary 
shapes for different values of plate thickness ¢ to center- 
to-center spacing /. It is seen that the radius of curva- 
ture varies greatly over each curved portion of the 
boundary, being relatively large at the center point and 
decreasing toward the ends. This variation in curya- 
ture is necessary in order to maintain constant field 
strength over the curved boundary. On the straight, 
vertical portions of the boundary, the field drops off 
rapidly from the corner. If the curved boundary had any 
other shape, the field would be nonuniform, attaining in 
some regions values higher than the constant value oc- 
curring on the boundaries of Figs. 3 and 4. Conse- 
quently, the boundary shapes in these figures are par- 
ticularly significant for the design of high-power filters 
requiring this basic configuration. 

The constant field strength Emax turns out to be very 
simply related to the uniform field Eo well below the 


array, as follows: 
J Drees Els 
= / ae (1) 
Eo t 


The square-root relationship makes Emax/Ho increase 
slowly as t/l is decreased. For example, if ¢// equals 
0.5, Emax/E£o will equal 1.414, and if ¢// equals 0.1, 
Emax/E9 will equal 3.16. Points x, y on the curved 


1H. A. Wheeler, “Conformal Mapping of Rounded Polygons by 
a Wave-Filter Analog,” Wheeler Labs., Inc., Great Neck, N. Y., 
Rept. No. 667; August 8, 1955. 
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boundary may be computed by means of the following 


parametric equations: 
L—% 
«= —— [p tan (p tan¢) — 4] (2) 
T 


Se [1 + (p2 — 1) sin’ 4], (3) 


y =S 
Tv 
where p=(1+#/l)/(1—t/l), and where ¢ is an inde- 
pendent variable. The co-ordinate system is shown in 
Fig. 5. As @ is increaséd from —7/2 to 7/2, the point 
x, y moves from the left end of the curve to the right 
end, resulting in a curve that is a symmetrical function 
of @. The coordinates of the end points are 
x1 1 M1 


11-—é¢/ 1+42/] 
Ss, —) = n ° 
t ae. eg aR 1—<2/1 


(4) 


The value of y;/t decreases from 0.318 for t//=0 to 0 for 
t/l=1, which results in the curve shapes being more 
flattened for the larger t// values, as may be seen in 
Fig. 4. 

There are two special cases of constant-field-strength 
boundaries that may be derived from (2)—(4). The 


+: 0.25 


Fig. 3—Arrays of rounded 180-degree corners shaped 
for constant electric field. 


0.4 “05 


Fig. 4—Shapes of rounded 180-degree corners for 
various values of ¢/I. 
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first is that of an isolated plate with rounded edge, 
obtained when //t-«. The curve for this case, which 
was computed previously by Wheeler," is plotted in 
Fig. 6. The end points are x.= +t/2 and y,=0.3182. 
This curve may be used in other situations than that of 
Fig. 6 as an approximation to the ideal curve, if the 
edge of a plate is removed from any other boundary sur- 
face by a distance equal to at least several plate thick- 
nesses. 

The other special case is that of an isolated slot with 
rounded edges where the two sides of the slot are at the 
same potential. The boundary shape, obtained when 
t/l—1, is shown in Fig. 7. It is seen that the curve falls 
monotonically to the left, but approaches zero slope as 
x—»— «, Eq. (1) shows that the field strength on the 
curved boundary is equal to the uniform field strength 
far below it. However, because of its infinite extent, the 
boundary must in practice be altered at some point, in 
which case a somewhat greater maximum field strength 
will occur. 


a 


~My xy 
> x 


Fig. S—Coordinate system for array of 180-degree corners. 


O.318t 


Fig. 6—Isolated rounded 180-degree corner plotted to scale. 
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WALL 


Fig. 7—Limiting case of rounded 180-degree corner for t/l1. 
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IV. ROUNDED 90-DEGREE CORNER NEAR 
ELEctTRIc WALL 


Fig. 2(c) shows a rounded 90-degree corner near an 
electric wall, and its symmetrical equivalent. Two 
curved-boundary shapes are of interest: 1) a circular 


arc, and 2) a shape resulting in constant field strength 
on the rounded surface. 


A. Approximate Circular Boundary 


It is not feasible to solve exactly the case of a 90- 
degree circular-arc corner, but a solution has been ob- 
tained for a shape that is a good approximation to a cir- 
cular arc.’-“ The actual curved boundary obtained in 
this case is shown in Fig. 8 for two values of 7/b, where 
r and 0 are defined in the figure. These boundaries were 
plotted with the aid of formulas given by Weber.” 
(Note, however, that Weber’s plot of one of the bound- 


--aries'® is grossly inaccurate.) The formula for the ratio 


of Emax/Eo as a function of r/b is also given by Weber, 
and is plotted in Fig. 9. For example, Emax/Ey)=1.40 at 
r/b=0.6. Thus, in a rectangular waveguide containing 
such a corner, breakdown will occur at about half the 
breakdown power of the uniform waveguide itself. 

An abrupt step in height of a waveguide, or in di- 
ameter of a coaxial line, has an equivalent circuit con- 
sisting simply of a shunt capacitive susceptance at ref- 
erence planes coinciding with the step itself. Graphical 
data have been published giving the value of this sus- 
ceptance for sharp-corner rectangular-waveguide and 
coaxial-line configurations.'!7!8 If the corner is rounded, 
the discontinuity susceptance will be less than that of a 
sharp corner. A formula for the change in susceptance, 
AB, was previously derived by this author’? for the ap- 
proximate circular-arc case. The resulting plot of AB 
versus r/b is given in Fig. 10 for a step in height in 
rectangular waveguide. The total shunt susceptance of 
the rounded step is 


B=B8B 


r=0 25 AB, (5) 


where B\ ,-o 1s the susceptance of a sharp step.'® Note 
that AB is a negative quantity, and hence B is smaller 
than leer Rounding of the corner also results in an 


2 L, Dreyfus, “Zur Berechnung von Durchslags und Uberschlags- 
spannung,” Arch. Electrotech., vol. 13, p. 131; 1924. | ; 

18 J. D. Cockcroft, “The effect of curved boundaries on the dis- 
tribution of electrical stress round corners,” J. TEE, vol. 66, pp. 
385-409; 1928. : 

4 E. Weber, “Electromagnetic Fields,” John Wiley and Sons, Inc., 
New York, N. Y., vol. I, pp. 373-377; 1950. 

16 Tbid., p. a 

16 Tbid., p. 376. : 

17 N. Miscavibes “Waveguide Handbook,” McGraw-Hill Book 
Co., Inc., New York, N. Y., Rad. Lab. Series, vol. 10, pp. 307-312; 
1951. : 4 oe ; 

18 J. R. Whinnery, H. W. Jamieson, and T. E. Robbins, “Coaxial- 
line discontinuities,” Proc. IRE, vol. 32, pp. 695-709; November, 
1944, 

19 S. B. Cohn, unpublished note; March, 1958. 
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Fig. 8—Actual shape of approximately circular, 90-degree 
rounded corner near an electric wall. 
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Fig. 9—Plot of Emax/E) for rounded 90-degree corner near an 
electric wall; also, 72/7 for a uniform-field-strength boundary. 
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Fig. 10—Capacitive susceptance correction and equivalent 
circuits for rounded step in waveguide. 
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increase in stored magnetic energy. This may be taken 
account of by adding a series inductive element in the 
equivalent circuit at the reference plane of the step, or 
by a shift of one of the reference planes and modifica- 
tion of B. ‘These alternative equivalent circuits are 
shown in Fig. 10. Values of AB taken from the curve in 
Fig. 10 will be accurate as long as 7 is sufficiently small 
compared to the dimension (6:—0:)/2. Good results 
should be expected for 7 as large as one half that dimen- 
sion. It is also necessary that AX/Z¢. be small; e.g., less 
than 0.3 should suffice. 


B. Constant-Field-Strength Boundary 


A solution in the case of a rounded 90-degree corner 
near an electric wall has been obtained yielding the 
boundary shape having uniform field strength on the 
curved arc.!23 Cockroft gives formulas from which this 
boundary shape may be plotted, and shows in his Fig. 
19 one such boundary drawn to scale. The shape of the 
boundary is very much like that of Fig. 8, but com- 
pressed in height. The plot of Emax/Ho versus 7:/b for 
this case is included in Fig. 9. Also plotted is the ratio of 
r, and 72, the dimensions that define the horizontal and 
vertical extent of the curved boundary. As should be 
expected, the Emax/EH) curve for the uniform-field- 
strength case falls below that for the approximate circu- 
lar-arc case. The two cases coalesce as 1/b approaches 
zero. 


V. ROUNDED 90-DEGREE CORNER NEAR 
MAGNETIC WALL 


A rounded 90-degree corner near a magnetic wall is 
shown in Fig. 2(d), along with its symmetrical equiv- 
alent. As in the case of an adjacent electric wall, 
Section IV, two curved-boundary shapes are of particu- 
lar interest: 1) a circular arc, and 2) the shape resulting 
in constant field strength on the rounded surface. 


A. Approximate Circular Boundary 


The approximate-circular-arc solution discussed in 
Section IV-A has been modified to apply to an adjacent 
magnetic wall (see Part B of Appendix). The boundary 
shape is exactly the same as in Section IV-A, and hence 
Fig. 8 applies to the magnetic-wall case as well as to the 
electric-wall case if r/d=r/b. For the adjacent magnetic 
wall, the ratio of the maximum field strength on the 
boundary to the distant uniform field strength was 
found to be 


tan 


Se btN) \/ 1oErp, (6) 
0 

where A and # are parameters plotted by Weber" in his 

Fig. 28.8. Eq. (6) is plotted in Fig. 11. The maximum- 

field-strength point occurs at the end of the curved 

boundary farthest from the magnetic wall. 
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Fig. 11—Plot of Emax/E» for rounded 90-degree corner near a mag- 
netic wall; also, 72/m for a uniform-field-strength boundary. 


B. Constant-Field-Strength Boundary 


The rounded 90-degree corner shape having constant 
field strength Emax on the curved portion of the bound- 
ary has been solved for the case of an adjacent mag- 
netic wall [see Part B-2) of Appendix]. The curved- 
boundary end-point dimensions 7 and 7 are given by 


777s eat) a 
WD 


1 
Fae [tanh a/1 = 4? = 4/13], (8) 


Tv 


where J is a parameter with value between 0 and 1. In 

terms of A, the ratio of constant field strength Enax on 

the curved boundary to the distant uniform field Ep is 

ae 1 + r 

Ea Vee ©) 
0 —iX 


The resulting curves of Emax/Eo and 72/7; are shown in 
Fig. 11. Note that the Enax/E» curve falls below the 
corresponding curve for the approximate circular- 
boundary case, as would be expected. 

Detailed boundary curves for this case may be 
plotted from the following equations: 


a 1 
— =— {tan /u—1—JS/u— 1} 


d 1 


(10) 
aes! aes [a 
Fo Sl tanh VA Poy ae 


where d is the same parameter as in (7), (8), and (9), u 
is a variable in the interval 1<u<1/)X, and the coor- 
dinate system is as shown in Fig. 11. The boundary 
shapes yielded by (10) and (11) resemble those in Fig. 8 
very closely, if the horizontal and vertical scale factors 
are altered to make the end points correspond to the 
appropriate values of 7:/d and r2/d. 
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VI. RouUNDED 45-DEGREE CORNER 


The solution for a rounded 45-degree corner adjacent 
to an electric wall is outlined in Part C of the Appendix. 
This configuration and its symmetrical equivalent are 
‘shown in Fig. 2(b). The resulting curved boundary is 
neither a circular arc, nor a shape having constant field 
strength. For r/b<1, the curve has the appearance 
shown in Fig. 12. The formula for this curve is 


r 
ieee Ue 14], 


a a) 
where =x-+7y is a point on the curve with coordinate 
system as shown in Fig. 12, and tis a real variable in the 
interval —1<¢<1. The field strength along the curved 
portion of the boundary is given, in the case r/b<1, by 
F)1/2 
E(2)= ; Sonia (13) 
[Gp P+ A) P+ V2 DEM? 

in which &, is the field strength at the center point of 
the curve. Because of the assumption that r/b<1, 
E(z) is a symmetrical function of ¢ near the corner. Eq. 
(13) is plotted in Fig. 13, where it is seen that the field 
strength on the boundary is almost constant except very 
near the junction points between curved and straight 
parts of the boundary. For r/b<1, the ratio of field 
strength at the center of the curve to the uniform field 
strength in the parallel-plane region is 


E. 2 2\115/ \ Us 
By) (+ 4/2)? (=) (-) 

b 1/5 
wey 

Y 


Clearly, if the transitions from the curved-to-straight 
portions of the boundary were made slightly more 
gradual, the sharp rise at the junction points could be 
greatly reduced, with the maximum field strength on 
the boundary perhaps 1.1 times the value at the center 
of the curve. 

Because of the one-fifth-root dependence in (14), the 
increase in E,/Ep is very slow relative to 6/r. For ex- 
ample, E,/E)=1.35 when r/b=0.1, while for r/o=0.01, 
E,/E,=2.14. These are much smaller field-strength 
values than occur with 90-degree rounded corners. 


(14) 


VII. CONCLUSIONS 


In high-power filters, whether operated with high 
internal pressure or vacuum, the power-handling ability 
is limited by electric-field concentrations at corners. By 
increasing the radius of curvature of these corners, the 
power rating of the filter may be increased. The formu- 
las and graphs contained in this paper give quantitative 
data for various practical rounded-corner configura- 
tions. With these data, a filter designer can predict the 
power rating of a given structure or, if a particular 
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Fig. 12—Shape of rounded 45-degree corner drawn 
to scale for r/b<1. 
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Fig. 13—Electric field strength on boundary of the rounded 45- 
degree corner shown in Fig. 12 (¢ is a “length” parameter meas- 
ured from the center of the corner). 


rating is specified, he can select the necessary rounded- 
corner dimensions to meet that requirement. 

Some of the rounded corners considered in this paper 
are approximations to circular arcs, while others are 
shaped so that the electric field strength is uniform on 
the rounded portion of the boundary. The uniform- 
field-strength corners may be considered to be optimum, 
since any other boundary curve starting and ending at 
the same points would have a greater field strength 
somewhere on the curve. 

Machining with high precision the exact theoretical 
shapes treated in this paper would not be a straight- 
forward operation. Fortunately, minor deviations in 
shape will result in only small changes in the value of 
maximum field strength, as long as the boundary curve 
is smooth. 


APPENDIX 
A. Analysis of Array of 180-Degree Rounded Corners 


Following a suggestion of Wheeler,’ one may show 
that the periodic z-plane boundary in Fig. 14(a) is 
transformed into the real axis of the w=w-+jv plane in 
Fig. 14(b) by the following relation: 

dz 


— = /1—m’' tanw + jm tanw, 


(15) 
dw 


where m is a real constant between 0 and . The single 
period of the z-plane boundary between x=0 and x=T 
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maps into the w axis, also between 0 and T. Inspec- 
tion of (15) shows that dz/dw is constant 1n magni- 
tude but variable in angle for w real and in the range 
| 2| <tan-'(1/m). This range of w corresponds to the 
curved part of the z-plane boundary. The electric field 
strength in the z plane is proportional to | dw/dz| , and 
hence the field strength is constant along the curve. 
Further examination of (15) shows that this constant 
value of field strength is the maximum occurring in the 
cross section. Integration of (15) and substitution of the 
appropriate boundary conditions yields (1)—(4). 
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Fig. 14—Transformation for array of 180-degree 
rounded corners. 
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Fig. 15—Transformations for rounded 90-degree 
corner near magnetic wall. 


B. Analysis of Rounded 90-Degree Corner Near Magnetic 
Wall 


1) Approximate Circular Boundary: Weber™ gives the 
following transformation (with notation changed) re- 
lating the g- and ¢t-plane boundaries of Fig. 15: 


oer 


p+aAvi — q}, (16) 


{V/ 


dz Cy 
dt t 


where Ci, ~, g, and 2 are real, positive constants, with 
q>p. The formula for z is obtained by integration. 
Relations between #, g, and ) are established such that 
r1=1r2=r. These constants may be computed from his 
formulas as functions of 7/b, or read from his graph.2° 
Weber applied this transformation to the case of an 
electric-wall boundary adjacent to the rounded corner 
by assuming the real axis of the ¢ plane to be an elec- 
tric wall with a discontinuity in potential at point 
2. However, the same transformation may be used 
for an adjacent magnetic wall by transforming the ¢ 
plane into the w plane as shown in Fig. 15. The vertical 
magnetic wall in the w plane maps into the adjacent 
magnetic wall of the z plane. Since the electric field in 


20 Weber, op. cit., p. 375. 
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the upper-right-hand corner of the w plane is uniform, 
the electric field strength in the z plane is proportional 
to |dw/dz|, and was found to be 


We 2 


From this, (6) for Emax/Eo was derived. 

2) Constant-Field-Strength Boundary: Inspection of 
(17) shows that if \ is set equal to p/q, E(z) will be 
constant for ¢ real and in the interval p<t<g. When 
\=+/p/g is substituted in (17) and in Weber’s formulas, 
(7)—(11) result. 


(17) 


E(z) « 


C. Analysis of Rounded 45-Degree Corner 


By application of the Schwartz-Christoffel method, 
the z-plane boundary in Fig. 16(a) was transformed into 
the real axis of the ¢ plane by the following relation: 


dz CiI4 


— = 18 
dt t—p Oy 


After integration and substitution of boundary condi- 
tions, the following formula for z as a function of t was 
obtained: 


(19) 


where p is a positive constant. This transformation 
yields an abrupt 45-degree corner. To round the corner 
as in Fig. 16(b), (19) was modified into the following 
form: 


z= Az(’) + (1 — A)z(t”), (20) 


where 
eK — 09 
f= (’ aT — 6; and’ 2?” = (¢ \i+ da. (21) 
p p 


The changes in variable of (21) are such that t/=¢!’ =f 
at the critical points t=1 and », but not at ¢=0. At 
t'=0, t= pb1/(p+4:), while at i” =0, t= —pd:/(p—6)). 
As a result of this modification, a gradual change of 
slope from 0 to 45 degrees occurs with (20) for ¢ in the 
range —6:/(1—8:) <#<6,/(1+46)). 
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Fig. 16—Transformations for 45-degree corner. (a) Abrupt 45- 
degree corner. (b) Rounded 45-degree corner. (c) Transforma- 
tion between ¢ and w planes. 


The procedure to follow in determining the con- 
tants A, 61, and 62 for a given value of r/d is to let (20) 
take on the following values: 


6 bo 
Z2=-—r when # =0 and pas (22) 
p01 
1 } 6 é: 
= ae, when f/f = ees and 4° =0, (23) 
/2 p — b2 


Since these are only two conditions and there are three 
constants to be determined, one additional condition 
may be specified. This could be, for example, that the 
curve should be most nearly symmetrical about its 
center point. 

The case r/b<1 is of particular practical interest. In 
the limit r/b-30, the constants have values A =1/2 and 
6, = 62 =6, so that 6 becomes the only unknown. Further- 
more 6p so that (20) and (21) become 


z= 42(t — 6) + $2(¢ + 8). (24) 
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In the limit 0, (19) and (24) reduce to 


g> — 


{ (4 — 6)5/4-+ (¢ + 8) 5/4}, (25) 


omp*!* 
Eq. (25) may be used to calculate the shape of the z- 
plane boundary in the vicinity of the rounded corner, 
but it does not, of course, yield the adjacent electric 
wall. 
A discontinuity in potential occurs at point 1 in the 
z and ¢ planes, so that the electric field in the ¢ plane is 
nonuniform. The real axis of the t plane may be trans- 
formed into the w-plane boundary of Fig. 16(c) con- 
taining a uniform electric field by means of this dif- 
ferential relation 


dw 1 
ie eae 


(26) 


The electric field within the z-plane boundary near the 
corner may then be obtained by differentiation of (23) 
and combination with (26), with | | <p, 


dw Qe p'!4/b 
dz | (¢ — 6)/4 + (t+ 84] 


(27) 


E(z) « | 


The uniform field Hy) far from the corner is obtained 
from (18) and (26) with i>) and C= —b/rp"4 


dw 


dz 


a (28) 


Ey « 


Combination of (27) and (28) gives the following for- 
mula for E(g)/E» valid in the vicinity of the rounded 
corner for the case 7/b<1: 


Ee) _ apis 
Fo Ge 2 


(29) 


The maximum-field-strength points occur at = +6. At 
t=6, (25) yields 


x (26)5/4 d 1/4 (= ) ea (30) 
= an SS Z ° 
f Sr pelt Y Sar 


Eqs. (12)—(14) now follow readily from (25), (29), and 
(30), when 6 is set equal to unity. 
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A Periodic Structure of Cylindrical Posts in a 
Rectangular Waveguide* 


EDWARD E. ALTSHULER?, MEMBER, IRE 


Summary—The propagation characteristics of a rectangular 
waveguide loaded with uniformly spaced cylindrical posts (periodic 
structure) are investigated at a frequency of 2840 Mc. A qualitative 
discussion on the expected behavior of the effective guide wavelength 
of this type of periodic structure is presented, and it is shown that the 
presence of the posts reduces the guide wavelength of the waveguide. 
The guide wavelength is then measured as a function of post diame- 
ter, post depth, and post spacing; and curves enabling one to design 
periodic structures which have guide wavelengths in the region of the 
free space wavelength are presented. 


INTRODUCTION 


HE propagation characteristics of a rectangular 
‘| eee are modified when the guide is loaded 

with uniformly spaced posts. The purpose of this 
investigation is to determine how the effective guide 
wavelength of the periodic structure shown in Fig. 1 
varies as a function of post geometry and its depth of 
penetration. 

To understand the behavior of this periodic structure 
it is first necessary to examine the equivalent circuit of 
a single element of the structure. A cylindrical post of 
variable height and diameter in a rectangular wave- 
guide can be represented by the four-terminal network 
shown in Fig. 2. An approximate theoretical solution 
for this post has been obtained by Suzuki;! however, the 
final expressions for the reactances are very difficult to 
evaluate. Experimental results have been obtained by 
Marcuvitz? and although they cannot be scaled di- 
rectly to the waveguide and post dimensions being in- 
vestigated, they can be used to describe the qualitative 
behavior of the post well enough to analyze the periodic 
structure. 

In Fig. 2 7x, is a relatively small series reactance that 
increases slightly with post depth and diameter. The 
main effect of the post is to introduce a shunt reac- 
tance jx, across the waveguide. In Fig. 3 jx, is plotted 
as a function of post depth h/b for various post di- 
ameters d. As the post penetration is increased, the 
capacitive reactance of the post approaches zero (when 


* Received by the PGMTT, April 10, 1961; revised manuscript 
received, May 31, 1961. : 

{ Electromagnetic Rad. Lab., AF Research Labs., Office of Aero- 
space Research, USAF, Bedford, Mass. 

? M. Suzuki, “Circuit Parameters of a Tuning Post in a Rectangu- 
lar Waveguide and its Applications,” AF Cambridge Research Cen- 
ter, Bedford, Mass., AFCRC-TN-57-764 [Rept. R-591-57, PIB, 519 
AFCRC Contract No. AF 19(604)-2031]; July, 1957. chet 

2, N. Marcuvitz, “Waveguide Handbook,” M.I.T. Rad. Lab. Ser. 
McGraw-Hill Book Co., Inc., New York, N. Y., vol. 10, p. 271; 1947. 


Fig. 2—Equivalent circuit of post. 
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r 


Fig. 3—Shunt reactance of post. (Experimental curves 
from Marcuvitz.?) 
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it behaves as a short circuit) and then becomes induc- 
tive. It can be shown from the theorem of images (see 
Fig. 3) that the shunt reactance of the double post being 
investigated is equivalent to that of a single post in a 
, waveguide of half the height. 

The effective guide wavelength \,’ and characteristic 
impedance Z’ of a periodic structure are related to the 
guide wavelength \, and characteristic impedance Z of 
the unloaded waveguide as follows: 


2a Qa 1 Qa 
cos — L = cos - sin — L, (1) 
9g Ng Xa NG 
tan — L 
Zz! ws. , 
Z 1 (2) 
tan we Tb 


9 


where L is the distance between posts, and x, is the 
‘normalized shunt reactance of the posts. (It is as- 
sumed that x =0.) 

Eqs. (1) and (2) are valid only when the posts are far 
enough apart so that the coupling between them is 
negligible. At a fixed frequency (A,=constant), A,’ is a 
function of x, and ZL. For a fixed post separation of 
L/Xg, aS X%_ is made less negative, cos (27rL/\,’) ap- 
proaches —1. When cos (27L/),’) is greater than —1, 
it becomes indeterminant and this corresponds to a cut- 
off region in the waveguide. For a fixed L/),, therefore, 
A,’ can be decreased by only a limited amount regard- 
less how small x, is made. Upon examining the limiting 
condition where 

2rL 
cos = —l, 
AG 


however, it can be seen that \,’=2Z and so ),’ can be 
made smaller by decreasing the distance between posts. 
When this distance becomes too small, the coupling be- 
tween posts cannot be neglected, (1) is no longer valid, 
and ),’ cannot be further decreased. 

When x, is made inductive (positive), cos (27L/),’) 
approaches +1 and ),’ approaches infinity. If the posts 
become too inductive 


[cos (2rL/Xy’) > 1], 


the TEi) mode does not propagate. 

The behavior of the periodic structure of cylindrical 
posts can be summarized as follows: When the post 
penetration is small the post can be represented by a 
shunt capacitance, and the guide wavelength },’ of the 
periodic structure is less than the guide wavelength 
X, of the unloaded waveguide. As the post penetration 
is further increased, ),’ is further decreased until a cut- 


3 J. Brown, “The design of metallic delay dielectrics,” Proc. IEE, 
vol. 97, pt. 3, p. 45; January, 1950. 
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off region is reached and propagation of the TE) mode 
ceases. With a further increase in post penetration, the 
shunt reactance passes through a resonance (x,=0) and 
then becomes inductive. When the inductive reactance 
increases to a large enough value such that cos 
(2rL/X,’)<+1, propagation resumes and ),’ decreases 
from infinity down to a value larger than \, as the post 
is extended all the way across the waveguide. 


EXPERIMENTAL PROCEDURE 


The guide wavelength ),’ of a periodic structure of 
known length can be determined by measuring the phase 
of a traveling wave at the beginning and end of the 
structure if the guide wavelength A, of the unloaded 
waveguide is known. The phase measurement is made 
using the experimental setup shown in Fig. 4. The input 
signal is divided into two parts; one part travels through 
the periodic structure under test, and the other part 
travels through a slotted line that is used as a reference 
phase line. To obtain precise data it is necessary that 
both signals be traveling waves. Each waveguide is 
therefore terminated in a matched load. The input 
signals to arms 1 and 2 of the magic T are adjusted so 
that they are approximately equal in amplitude. A 
sharp null will appear at output arm 3 when the input 
signals are 180° out of phase and arm 4 is terminated in 
a matched load. The phase difference between two sig- 
nals is determined from the distance that the probe on 
the phase line (slotted line No. 1) has to be moved to 
produce a null at output arm 3. 


MATCHED 
LOAD 
SLOTTED 
LINE #1 
E—- PLANE 
TEE 
VARIABLE 
ATTENUATOR 
SLOTTED 
LINE#2 


PROBE 


DETECTOR 
# 2 


KLYSTRON POWER 
2 840 MC SUPPLY 


FERRITE 
ISOLATOR 


3 


DETECTOR 


Ya- WAVE PERIODIC 20-DB DIR 
TRANS. STRUCTURE COUPLER 


Fig. 4—Block diagram of equipment. 


The periodic structure consists of a uniform section 
of waveguide with adjustable posts and behaves as a 
uniform waveguide with a propagation constant and 
characteristic impedance different from that of the un- 
loaded waveguide. The phase of the signal from the di- 
rectional coupler at the end of the periodic structure is 
first measured with the posts at zero depth. The post 
penetration is then slightly increased and the meas- 
urement is repeated. The periodic structure has an 
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Fig. 5—Guide wavelength vs post depth for constant post spacing. 
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Fig. 6—Guide wavelength vs post depth for constant post diameter. 
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electrical length of 630° (1.75\,) with zero post pene- 
tration. As the post depth is increased, the shift in 
phase increases. Since this corresponds to an increase 
in the electrical length of the structure, the guide wave- 
length is decreased. The ratio of freespace wavelength 
Xo to the guide wavelength A, of the unloaded structure 
is 0.875. The ratio of Ao to A,’ is 


% - Ao 
—=().8/5|/1-+ ‘ 
Xa 630° 


where A@ is the change in phase in degrees due to the 
presence of the posts. 

Measurement of the input VSWR of the periodic 
structure checks the relationship. 


Tv 
tan a |b 
“\Qq 


This is possible because the input VSWR is composed 
of reflections from the interface from Z to Z’ and the 
interface from Z’ to Z. Since the characteristic im- 
pedances are assumed real, the VSWR that is due to 
each interface is R/R’ for R>R’. The resultant input 
VSWR can vary from a minimum of 1.0 to a maximum 
of (R/R’)?, depending on the electrical distance be- 
tween interfaces. It is 1.0 when the interfaces are an 
even number of quarter-wavelengths apart and (R/R’)? 
when they are an odd number of quarter-wavelengths 
apart. Due to the reflections at the interfaces, there is a 
slight difference between the phase angle of the trans- 
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mission coefficient that is measured and the electrical 
length of the periodic section. Since the reflection 
coefficients are small, the error introduced is negligible. 


DISCUSSION OF RESULTS 


The guide wavelength \,’ of the periodic structure 
shown in Fig. 1 was measured as a function of post di- 
ameter d, penetration p, and spacing L, at a frequency of 
2840 Mc (Ao =4.159 inches). In Figs. 5 and 6, Ao/A,’ is 
plotted as a function of p/A». Each set of curves in Fig. 5 
corresponds to a fixed post spacing LZ. Each set of curves 
in Fig. 6 corresponds to a fixed post diameter d. It can 
be seen that Ao/A,’ increases slowly with an increase in 
post diameter and depth and a decrease in post spacing. 
These results agree very well with those expected from 
the qualitative discussion that was presented earlier. 
Since this investigation was conducted primarily in the 
region where \,’=Ao, the actual cutoff conditions of the 
periodic structure were not established; but the fact 
that the slopes of the curves become very steep for large 
post diameter and penetration indicates that a cutoff 
condition is rapidly being approached. 

The results of the input VSWR measurements were 
found to be consistent with those expected from (2). 
The assumption that this periodic structure can be 
treated as a uniform waveguide with a real charac- 
teristic impedance is therefore valid. 
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Modes in Rectangular Guides Loaded with a 
Transversely Magnetized Slab of Ferrite 
away from the Side Walls* 


G. BARZILAT}, sENIOR MEMBER, IRE, AND G. GEROSA, ASSOCIATE, IRE 


Summary—The characteristic equation describing the general 
modal spectrum for a rectangular guide partially filled with a slab 
of ferrite transversely magnetized and situated away from the side 
walls is derived. This equation is numerically solved for particular 
cases and for modes of zero and first order with respect to the depend- 
ence along the direction of the dc magnetic field. Some experiments 
to verify the theoretical results are presented and show good agree- 
. ment with the theory. 


I. INTRODUCTION 


HE authors in a previous paper! have discussed 
dike modal spectrum in rectangular guides com- 
pletely filled with transversely magnetized ferrite; 
in a successive paper? they have discussed the modal 
spectrum in rectangular guides loaded with a slab of 
transversely magnetized ferrite against one side wall. 
The purpose of this work is to report on further 
theoretical study on a more general structure of a 
rectangular guide loaded with a slab of transversely 
magnetized ferrite away from the side walls and to de- 
scribe some experiments carried out to verify the theory. 
The only analyses available up to date for this struc- 
ture are, to our knowledge, relative to modes with no 
dependence along the direction of the dc magnetic 
field.*~* 
The method of solution we have used is similar to the 
one already discussed,! but the derivation of the char- 


* Manuscript received by the PGMTT, February 20, 1961; re- 
vised manuscript received, June 15, 1961. This paper was presented 
at the 13th General Assembly of URSI, London, Eng., September 
5-15, 1960. The theoretical investigation reported herein was spon- 
sored by the U. S. Air Force under Contract No. AF-61 (052)-101. 
The experimental part was made possible by the support of the 
Consiglio Nazionale delle Ricerche of Italy. 

+ Istituto di Elettronica dell Universita di Roma, Rome Italy. 

1G, Barzilai and G. Gerosa, “Modes in rectangular guides filled 
with magnetized ferrite,” L’Onde électrique, 38° Année, N° 376ter, 
Suppl. Spécial—Congrés International Circuits et Antennes Hyper- 
fréquences, Paris, France; pp. 612-617; October 21-26, 1957. 

2G. Barzilai and G. Gerosa, “Modes in rectangular guides par- 
tially filled with transversely magnetized ferrite,” IRE TRANs. oN 
ANTENNAS AND PROPAGATION, vol. AP-7, pp. S471-S 474; December, 
1959. Also see Istituto Electrotecnico dell’Universita di Roma, Tech. 
Note No. 1, Contract No. AF 61(052)-101; June 3, 1959. 

3M. L. Kales, H. N. Chait and N. G. Sakiotis, “A nonreciprocal 
microwave component,” J. Appl. Phys., vol. 24, pp. 816-817; June, 
1953. Erratum: p. 1528; December, 1953. . 

4B. Lax, K. J. Button and L. M. Roth, “Ferrite phase shifters in 
rectangular waveguide,” J. Appl. Phys., vol. 25, pp. 1413-1421; 
November, 1954. : ; 

5B. Lax and K. J. Button, “Theory of ferrites in rectangular 
waveguide,” IRE Trans. ON ANTENNAS AND PROPAGATION, vol. 4, 
pp. 531-537; July, 1956. ; : 

6 A. L. Mikaelyan, “Magneto-optic phenomena in a rectangular 
waveguide containing a ferrite plate,” Izvest. Akad. Nauk. USSR, 
Otdel. Tekh. Nauk., No. 3, pp. 139-144; 1955. 


acteristic equation is very involved and the relative 
numerical computations, even by using a medium speed 
electronic computer, require considerable length of 
time. We have therefore limited our numerical in- 
vestigation to some cases which we thought to be of 
interest. 


Il. THE CHARACTERISTIC EQUATION 


Let us refer to the parallel plate guide represented in 
Fig. 1. We shall assume the walls to be perfectly con- 
ducting. 


NS 
? 
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Fig. 1—Geometry of the parallel plate guide partially 
filled with magnetized ferrite. 


Let us assume the dc magnetic field, of sufficient in- 
tensity Hy to saturate the ferrite, to be directed along 
the z axis. We shall assume for the ferrite region a 
scalar dielectric constant €9 € and a magnetic tensor 
permeability u given by the following expression (time 
dependence exp [jwt]): 


Mi jue O 
u— Ho jue Mi O ) 
0 ie 
where 
p Tp Mo w 
re pees ae rpm he terest 


and Mo is the intensity of the saturation magnetization, 
® and w)»=—vyH) are the applied and the resonant 
circular frequencies, y is the gyromagnetic ratio for the 
electron, fo and €) are the permeability and the dielectric 
constant of the vacuum. 

As discussed in detail in our previous work,! to con- 
struct our modal solution we shall consider for the three 
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regions of width 6, by and 49’ fields having spatial de- 


pendence of the form: 
exp [j(Aaw + ayy t bas) |. (1) 


We shall add to the propagation constants ky, ky, &: a 
subscript f or 0 in order to refer to the ferrite or to the 
vacuum regions, 

The propagation constants will be measured in units 
of © \/uoep ANd consequently lengths will be measured 
assuming as unit 1/@ Vv Moo: 

In order to satisfy Maxwell’s equations, we must 


have in the ferrite region 
pyit + [Cui + 1h? — e(ur? — we? + ar) Je + k.;' 
— Qurek.p? + A(ur® — ps?) = 0, (2) 


where 
&= ky + ky’, (3) 
and in the vacuum regions 

Ray? _ Ryo + | ts an (4) 


Each modal solution of our problem is labelled by a 
pair of values &., 8. In order to satisfy the boundary 
conditions at the interfaces we must have: k.;=k.» =k: 
and k.; =k. = kz. 

For the ferrite region (2) yields, for a given k,, two 
values of #*:42 and &*. From (3) we obtain therefore two 
values of &,7:%,.°7 and &,9*, for a given &,. Consequently 
we express the field in the ferrite region as a superposi- 
tion of four fields of the form (1) having the same &, 
and &, and arbitrary amplitudes. 

In the vacuum regions for a given pair k., k, we have 
only one value of &,o*, but, since a general field can be 
expressed as a superposition of TE and TM waves, for 
each vacuum region we shall have again four arbitrary 
amplitudes. 

By imposing the tangential component of the electric 
field to be zero at the metallic boundaries and the tan- 
gential component of both the electric and magnetic 
fields to be continuous at the interfaces between vac- 
uum and ferrite, we obtain a system of twelve linear al- 
gebraic homogeneous equations in the twelve unknown 
amplitudes. By setting equal to zero the determinant of 
the coefficients, we obtain the following characteristic 
equation: 


De 21MM & ama (bp) 
(i,k) (lym) ss 3 - = 
-. i,x(b0) Xim(— bo’) —— 0, (S) 
where the index pairs 7, k and /, m can assume the values 
$2241 331 A 23° 24-34 and (7, Rk). represents the 
pair which together with 7, k complete the set of first 


four integral numbers; and similarly for (1, m).. 


The functions X and Y have the following expres- 
sions: 


Xix(bo) = Crees — Cie, 
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where 

Cn = Re(1 — Re )Ryo Sin Ryodo Ci. = O 

¢y) = - RRA Ryo sin Ryobo 62> Ryo sin Ryobo 

¢33 = 0) 633: = (1 = k,”) cos Ryobo 
é) = RRyo® cos Ryobo C42 = kzk: cos Ryobo 


and 

— Ryrkyo(Cie,11Cim,.22 + Cik,22Cim,11) 

+ sin kyib; sin Ryob;(A ix,12A im, 12 

+ hy Ryo? Bix12Bim12 + RyV?Cx,21Cm,21 
+ Ryo? Ck 12C im, 12) 

+ RyRy: COS Ry1b; CoS Ry2b;( Bix ,12A im, 12 
+ Aix r2Bima2 + Cix,21Cimr2 + Cie,12Cim,21) 
+ ky cos ky1b; sin Ryobs[Cix,21 im, 12 

— A,12Cim21 + Ry2?(Cie,12Bim,12 

— Bik,12Cim,12) | 

+ Ry2 sin Ryiby cos Ryaby* [A sk,12Cim.12 

age Cx 2A Im,12 ={- ky? (Bik s2Cim,21 

— Cix,21Bim,12) ], 


Y 6.8) tam) (by) a 


where 
Aas = @ndye — Ap Ai2 
Byes = biubes — berdi2 
Cre = Giiber — Aidan 
Cre .22 = Giadkg — Axgdie 
Cxx21 = Arby — axabin 


Craze = Ginbkg — axrdie 
and similarly for 7, m, and: 


ait bok. (e€ = $s") = a2 


bu = — (ui — Ikke = dis 

dzy = — pole — ty*)Re G22 = — pole — fe”) Rez 

ba = kk. — mile — th’) bos = kz* — pi(e — Ze?) 
as, = ke[k2 — (€ — t1)] as = ky?[k.2 — (e — ty?) ] 


— e[k.? — wile — t1)] 
bsi =Q= bso 
G4 = — Ri(uie—k2— 12k, ag = — ke(uie — kz? — te?) he 


bar eS Rig oes bao. 


— ¢[k.? — wi(e — f2)] 


It can be seen that (5) contains only even powers of R,. 
It can be verified that the solution possesses reflec- 
tion symmetry along the s direction. Therefore a solu- 


tion for a rectangular guide of height a can be obtained 
by choosing 


i= 0 Tope (6) 
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For m=0, (5) breaks down into two equations cor- 
responding to TE and TM modes. TM modes cannot 
exist in the rectangular guide obtained by closing the 
structure of Fig. 1 with two metallic walls normal to the 
‘z-axis. TE modes can exist and the relative charac- 
teristic equation is the same as the one given by others. 
From m0, (5) does not in general break down into two 
equations and the resulting electromagnetic field can- 
not be resolved into TE and TM modes. 

The solution does not possess reflection symmetry 
along the x axis. In addition, (5) contains in general 
odd powers of kz, and the solution therefore is not 
reciprocal, 7.e., a solution +, does not necessarily 
imply the solution —k,. 

It should be noted that obvious solutions of (5) are 
kyi=0; ky2=0; kyo=0. These solutions however corre- 
spond to zero amplitude fields. It can be verified that 


hk. — wile — t1,2”) 
= € 
ky? — (€ — ty,2?) 


are always solutions of (5), but these solutions also cor- 
respond to zero amplitude fields. 

We note that by interchanging by with bo’ and k, 
with —k,, (5) remains unchanged as it follows from 
symmetry considerations. 

To discuss (5), it is convenient to find the asymptotic 
behavior of some of the solutions when bo’—0. More 
exactly we shall assume 0b; and bo finite and by letting 
by’—0 we shall look for solutions with real k, such that 
| Re| — oo. From (5) we obtain 


ke 


a eS 
a Mi(ui + 1) lee [2 


tanh | k.| bo’ = 


(7) 


z= 


oS ire 
Mit Take 


From (7) it is apparent that asymptotic solutions 
are possible only when the second member of the equa- 
tion is positive. In such a case the asymptotic solutions 
are represented by hyperbolas. It should also be noted 
that (7) is independent of k, and therefore the asymp- 
totic solutions considered are the same for any kz, 1.e., 
in the case of the rectangular guide for any m. 


III. NUMERICAL ANALYSIS 


There are several parameters which determine the 
solutions of (5), namely: quantities characteristic of the 
ferrite medium, 7.e., Mo and e; quantities describing the 
structure, 7.¢., a, b, by and bo’ (bb =b —by— bo’) ; impressed 
quantities, z.e., w and Ho. The three quantities Mo, 
w and H, enter our problem through the adimensional 
parameters p and 7, so that the actual parameters to be 
considered are p, €, T, a, 0, by and by’. 

When a set of these parameters has been chosen, (5) 
in virtue of (2), (3) and (4) becomes a relation between 
k, and k2. We shall fix k. by means of (6), which allows 
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us to label our modal solutions according to the integral 
values given to m. We shall call, therefore, modes of zero 
order those corresponding to k,=0, and modes of first 
order those corresponding to k,=7/a, etc. 

For a given mode order, (5) determines the relative 
propagation constants k,. In what follows we shall only 
look for real values of k,, 7.e., for unattenuated pro- 
pagating modes. 

With respect to the field configuration in the cross 
section of the guide the mode order determines the z de- 
pendence. The y dependence is determined by the 
values of Ry, Ry. and kyo, which can assume real or purely 
imaginary values. In fact, by assuming k, and k, real, 
(4) states that kyo? is real and (2), once solved with re- 
spect to ¢”, shows that ¢,;? and t.? are both real. 

The solutions of (5) can be classified by dividing the 
field of variability of 7 into six regions, exactly in the 
same manner as we have done in our previous work? rela- 
tive to the slab of ferrite against one side wall. 

We have numerically solved (5) for three cases which 
we thought to be of interest. Two of such cases are 
those corresponding to Figs. 5 and 6 of the quoted 
paper.” 

The results of the numerical analysis are recorded in 
the diagram of Figs. 2-4, which show k, vs by’ for the 
values of the parameters indicated in the figures and for 
modes of zero and first order. 

In the diagrams of Figs. 2-4 the various zones de- 
limited by the straight lines k,:=0, ky2=0 and ky» =0 
have been shaded in different ways in order to recognize 
at a glance when k,1, Ry», and kyo are real or imaginary. 
It is understood that when one of the three typical 
shadings indicated is present the relative k, is real and 
when it is not it is imaginary. For instance, when no 
shading exists the three k,’s are all imaginary; when all 
three shadings are present the three k,’s are all real, and 
so on. It should be noted however that for modes of zero 
order k,1 is associated with TM modes, which have zero 
amplitude in the rectangular guide. The shading rela- 
tive to k,, has been therefore omitted in the diagrams 
relative to zero order modes. 

The most interesting feature appearing from the 
diagrams is the fact that for suitable dimensions of the 
guide and of the ferrite slab unidirectional propagation 
exists for 6)’ =0 and, as soon as the slab of ferrite is 
moved sufficiently away from the wall, for zero and first 
order modes no propagation can occur in either direc- 
tion. It is reasonable to assume that higher order modes 
exhibit the same behavior. Experimental evidences de- 
scribed in the next section seem to agree with this 
theoretical result. 

In Fig. 4 we have indicated with bo.’ the cutoff dis- 
tance for modes of zero order. The cutoff distance for 
modes of first order is practically the same. 

With reference to the zero order modes of Fig. 4 we 
have investigated the sign of the group velocity for the 
region where there are two propagating modes. We note 
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Fig. 2—The mapping of the solutions of the characteristic equation 


for the indicated values of p, ¢, 7, @, 0 and by. These values may 
be taken to correspond to 
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Fig. 5—(a) T-shaped waveguide built to verify the unidirectional 
character of a guide of suitable dimensions loaded with magne- 
tized ferrite. (b) Experimental curves of the attenuation below the 
input level vs the external dc magnetic field for the structure and 
the frequencies indicated. All the dimensions are in millimeters. 


that, because of the normalization assumed, the group 
velocity u, is given by 
1 


Bea (Tes 
Mo€o ae 


Ug = 


By using the asymptotic expression (7) it is easily 
seen that by increasing the frequency, 1.e., 7, TR, de- 
creases and vice versa. On the other hand for by)’ =0 we 
have computed rk, for two values of the frequency 
slightly above and below the frequency relative to Fig. 4 
and we have found that the group velocity is positive. 
By using these results we can conclude that the two 
modes corresponding to the same value of bo’ have group 
velocity of opposite sign. It seems reasonable to extend 
the validity of this result to modes of the first order. 


IV. EXPERIMENTAL 


The experiments carried out had two different aims: 
1) to prove the unidirectionality of a structure corre- 
sponding to the case of Fig. 4 for bo’ =0; 2) to verify the 
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Fig. 6—(a) Waveguide section built to investigate the behavior of a 
guide loaded with a slab of magnetized ferrite. (b) Experimental 
curves of the attenuation below the input level vs bo’ for the values 
of the frequency and the external dc magnetic field indicated. All 
the dimensions are in millimeters. 


general behavior of a structure corresponding to the 
case of Fig. 4 when 6,’ is varied. 

For the experiment 1) we have built a T-shaped 
waveguide [shown in Fig. 5(a)] whose dimensions are 
given in Fig. 5(b). By sending the RF energy in the 
input arm we have recorded in Fig. 5(b) the atten- 
uation for both the output arms below the input level 
vs the external dc magnetic field.? This has been 
done for the three different frequencies indicated in 
the figure. From the experimental results of Fig. 5(b) 
the unidirectional character of the structure considered 
is apparent for a wide band of frequencies. 

For the experiment 2) we have built a guide section 
[shown in Fig. 6(a)] whose dimensions are given in 
Fig. 6(b). By means of insulating rods and screws the 
slab of ferrite could be moved parallel to itself inside 
the guide. By sending energy in the input arm we have 
recorded in Fig. 6(b) the attenuation below the input 
level at the output arm vs the distance of the slab 
of ferrite from one of the walls, for the two opposite 
values of the dc magnetic field, which correspond to 


7 Note that the de magnetic field assumed for the theoretical cal- 
culations is the internal one, which in the experiments is not uniform 
and its average is smaller than the external. 
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forward and reverse propagation. In the figure it is 
shown the attenuation for the empty guide, 7.e., for the 
same guide without the slab of ferrite. 

Let’s follow the forward propagation experimental 
curve from left to right. When the slab of ferrite is 
against the left side wall the attenuation is at a mini- 
mum. By moving the slab of ferrite away from the wall, 
after a region of low attenuation, the signal goes very 
rapidly below the level of the attenuation of the empty 
guide. This is justified by the theoretical results of Fig. 
4 and by the considerations about the group velocity. 
In fact, beyond the distance bo.’ no propagation can exist 
and below such a distance always exists a propagating 
mode with group velocity in the forward direction. 

Moving the slab further away from the left side wall, 
the signal remains at a level below the empty guide 
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attenuation level until we reach a distance approxi- 
mately equal to b—b;—bo-’ when the attenuation begins 
to decrease and then, after it has reached a minimum, 
increases again above the empty guide attenuation. 
This last behavior is easily explained by the theoretical 
results, since beyond the distance )—};—bo.’, energy 
begins to pass in the forward direction through modes 
having positive group velocity. However, since these 
modes have propagation constants going to — © as the 
slab of ferrite approaches the right-hand side wall and 
cannot therefore allow propagation in the limit, the 
minimum of the attenuation experimentally found is 
explained. 

From the preceding discussion we can conclude that 
there is good agreement between theory and experi- 
ments. 


Higher-Order Evaluation of Electromagnetic 
Diffraction by Circular Disks* 


W. H. EGGIMANN} 


Summary—tThe problem of the diffraction of an arbitrary electro- 
magnetic field by a circular perfectly-conducting disk has been 
solved by using a series representation in powers of kK=27/) and the 
rectangular disk coordinates. The surface current density is given in 
terms of the field and its derivatives at the center of the disk. General 
expressions for the electric- and magnetic-dipole moments, the far- 
field and the scattering coefficient for the case of a plane wave at 
arbitrary incidence are presented. The calculations agree with results 
published by other authors. A bibliography of the most recent publi- 
cations on this problem is included. 


I. INTRODUCTION 


HE problem of the diffraction by a circular con- 

ducting disk (or the complementary problem for 

a circular aperture in an infinite plane conducting 
screen) has occupied many workers in the field of dif- 
fraction theory. The problem can be formulated as 
follows: 


1) the electromagnetic field has to obey Maxwell’s 
equations, 

2) the boundary conditions on the surface of the disk 
have to be fufilled, e.g., for a perfectly-conducting 
disk the tangential electric field must vanish, 

3) the edge conditions [48] at the rim of the disk have 
to be obeyed; they require that the field energy 
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received, June 8, 1961. This work has been sponsored by the Elec- 
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remains finite, or that the energy density has to 
be integrable over any finite space. This leads to 
the requirement that the normal component of 
the electric field increases not faster than (1/r)!/? 
where ¢ is the distance from the edge, 

4) Sommerfeld’s radiation conditions [47] have to be 
fulfilled. 


In this paper a power-series solution in (ka) valid for 
the small disk problem (a <)/2z7, where a =disk radius, 
\=free-space wavelength) and an arbitrary incident 
field is given. It is essentially an extension of a proce- 
dure described by Bouwkamp [45]. The surface cur- 
rent density on the disk up to the third-order approxi- 
mation in (ka) is calculated in terms of the electro- 
magnetic field and its derivatives at the center of the 
disk. From these results expressions for the induced 
electric and magnetic dipole moments and the far-zone 
fields are derived. The scattering coefficient for a plane 
wave at arbitrary incidence has been calculated in agree- 
ment with formulas given by Lur’e [19] and Kuritsyn 
[20]. The essential advantage of the expressions ob- 
tained in this paper is that they can be used for any 
primary field. This is important in the case where in- 
teraction between several disks is considered. If the 
spacing between the disks is not large compared with 
the wavelength, the interaction fields cannot be ap- 
proximated by a plane wave and the interaction be- 
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between the induced dipole moments and higher-order 
multipole moments have to be taken into account. 
Some of this work will be reported at a later date. 


II. PowEr-SERIES SOLUTION 


General Formulation 


We shall consider a perfectly-conducting circular 
disk of radius a with its axis along the g direction as in 
Fig. 1. The magnetic vector potential A(R) for the dif- 
fracted field is given by 


—jkr 


Mo ee 
= ef se : 


where J(p’) is the electric surface current density and 
the integration is carried out over the disk D. Rand p’ 
are the coordinates of the field point P(x, y, z) and the 
surface element dS’ respectively, while r is the radius 
vector from dS’ to P(x, y, 2). The time dependence 
et is omitted throughout. The scattered field is found 
from the following relations: 


A(R) 


dS’ (1) 


1 
H=—VXA (2) 
Ko 
VV-A 
E = — jwA +. . (3) 
J@HOED 


The boundary conditions are satisfied if 


E,(R) = — E,*(R) 
E,(R) = — £,'(R) on the disk (4) 
HAR) = — H.(R) 


where the superscript 7 indicates the incident field. We 
now express (2) and (3) in rectangular coordinates and 
combine them with (4) 
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Laplace operator. 
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In (6) we used (5) and Maxwell’s equations for E,’. A 
similar expression as (6) is obtained for ,. These equa- 
tions can now be written in the form 


aH, 
VaAs aii RA, = 180 
OZ 

OH,’ 
VayrAy + RA, = — wo on the disk 

aA, dA, ” S 

a; ed 2” / 
ax oy : 


where k=w({o€o)!” is the free-space wave number. 
Eqs. (7) are the basic differential equations which we 
shall solve for the vector potential A(x, y) on the disk. 
Knowing A(x, y) we then obtain the surface-current 
distribution J(x’, y’) from the integral equation (1). 
In order to obtain a unique solution the conditions 
1)—4) in Section I have to be fulfilled. Conditions 1) and 
2) have been taken care of in (2)—(4). The edge condi- 
tion 3) will be fulfilled by assuming a suitable current 
distribution on the disk, and the radiation conditions 
4) are secured by (1) for any finite current distribution. 
In the following we attempt to find a power-series 
solution for the current J in terms of (ka), which is ex- 
pected to converge well for small disks where (ka) 
=2r(a/d) <1 or a<)d/2ar (A=free-space wavelength). 
First, we consider an expansion of A, J and the right- 
hand side of (7) in powers of (k). We obtain the follow- 
ing expressions: 


JH PRP ERP TBP HRTF oo (8) 
Jer = P+ RT! = jr) + RP = jr = 1/2) 
+ RJ? — jr? — 1/2P° J + 1/6jr8]°) 
+ AT! = jr — 1/2 P+ 1/657 P 
+ 1/24rJ) +++ (9) 
A = A°+ kA! + kA? + RAP 4+ BATH (10) 


where 
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For simplification we adopt the following notation: 


Ao = A? 
Ai = All — jA% 

A? = A® — jA? — 1/2A” 

A? = A®® — 7A — 1/2A¥ + 1/6jA% 


At = At — jA™ — 1/2A%4 + 1/6jA™% + 1/24A% (12) 


where the partial vector potentials A” are defined by 
(11) and (12). 


oH, ' 
Lo - = $9 + RS'+ k2S? + RS? --- 
z 
aH! 
— 0 ‘ =T+tRT!4+ PT? 4+ RTI+ .--- 
Zz 


—pwoH, = U°+ kU + RUPE + UE. (3) 


The coefficients S”, 7”, and U” in (13) are, of course, 
functions of the field coordinates x and y and can, there- 
fore, be written as 


S” = Sot + Syme + Syry + Syta® + Saray 


Oey ics ea 

[Pores GSS IE ae Jie a ARN aes TE nese 
Sela = 

CU Ua U sty Uist Use, 
ay nye tl ee (14) 


The coetiicientsS,,”, 1,,° and U,,” are constants and 
can be calculated from the primary field. For the term 
k"x"y’ we obtain 


i Gre ge OH 
Sn” = Mo ee 
mir!s! Ax™dy® Ok” Oz k=0 
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Where all derivatives are taken for x=y=z=k=0. 
The subscript m can be found by inspection of (14) or 


from the relation 
r+s 
m = ( SS ) ase 


t=1 


(15a) 


We now substitute (10) and (13) in the basic equa- 
tion (7). Equating equal powers of k we obtain for the 
nth-order approximation, 7.e., for k", the following rela- 
tions: 
VWAEe + A,” = sn 
VA = A eT 
OA," GAS 
Ox oy 


(16) 
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It is interesting to note that A” does depend on the 
value of A”. Consider now (16) in terms of the partial 
vector potentials as defined in (11) and (12). It is read- 
ily seen that the second term in these expressions, 
AD”, is independent of the field coordinates x and y, 
so that all its derivatives vanish identically. Summariz- 
ing we find that in the calculation for A” the expressions 
for A”! and A” or, in terms of the surface currents, 
the term J”~, do not appear. From this follows that the 
zero- and first-order approximation for J are not re- 
lated to each other and can be evaluated from the ex- 
pansion coefficients of the primary field only. For higher- 
order approximations, however, we have to use the re- 
sults of the previous calculations. This can be clearly 
seen, if we rewrite (16) and keep only the term A” on 
the left-hand side. We then obtain 


V?A,™ = yn (17a) 
V2A, = Ww (17b) 

0A,” GA, 
= = Xn (17c) 


Ox Oy 


where V”, W” and X™*” are functions of the field co- 
ordinates x and y and can be written as follows: 


Ve = Vet aa fee Very 4 ee ey 
+ Vey? t.-- = Vp erys 

We = Wor + Wirty + Worx + Wty? + Waryx 
Wet = >a Wnryrxs 

AX” = Xo Xie Noty + Xe a” ay 
+ X5"y" ae —— Ss age se 


ca 


(18) 


Again m is found from (15a). 
Vin®, Wr and Xm" can be calculated from the co- 
efficients S,”", T,” and U,,." for the primary field and 
from the results of the (#—2) and lower-order approxi- 
mations as will be shown explicitly in (29). The 4,”" 
and A,”" can then be obtained from (17) by straight- 
forward integration in form of power series in x and y. 
The main problem is now to find a solution for the 
integral equation for the current J”. This equation is 
obtained from (11) and (12): 
A(x, y) 
/ / 
2 nig es 
4dr J p [eo = ey a 


(19) 


Formal Solution for the Surface Current Distribution 


The kind of integral which appears in (19) has been 
investigated by Bouwkamp [45], [49], [50]. He found 
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that if we write the current distribution in the follow- 
ing form: 


Ho a(x’, 9’) 
a ee aan 

7 Tacs =p) 422 
M ful#’, 9’) 


0 
BPG | a, 1 ee 
As le’, 9’) (a2 — p!2)3/2 


(Peo y)2 (20) 


where f, and fy are polynomials in x’ and y’, the integral 
in (19) is also a polynomial in x and y, and is of the same 
order as f(x’, y’). Table I in the Appendix gives a solu- 
tion for the following integral: 


G(x, y) 


. 1 f oe ao) 
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The integration is over the surface D of the disk; 
f(x’, »’) is a polynomial in x’ and y’. Only terms up to 
the fourth power in x’ and y’ have been calculated. 
Computations of higher-order terms are possible yet 
they become, in principle, extremely tedious. It is note- 
worthy that G(x, y) is also a polynomial of the same 
order as f(x’, y’). 

In order to obtain a unique solution, we have to as- 
certain that the edge conditions are fulfilled. It can be 
shown that they are equivalent to the requirement that 
the current component normal to the edge vanishes at 
the rim of the disk. The radial current is 


J,=J,cosg’ + J,sing’ = (1/p’) («Jae + y'Jy) 
and thus we must have 


cae «fe + y'fy 
—— tn ——— ———— = ) 
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(22) 


We further require that the total charge on the disk 
remains finite. The charge density o is given by 


j jp 4 f 
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The last term in (23) is only integrable over the disk if 
Gis 19 fx) = Dix’, ye — vp”) (24) 


where D(x’, y’) is also a polynomial in «’ and y’ as will 
be shown in the following calculations. This condition, 
however, is sufficient to satisfy (22) and is, therefore, 
equivalent to the edge condition. 

The current-distribution function is now expanded in 
terms of & as follows: 


f=P+R + PP ++: :. (25) 
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We then write for the polynomial 

In(x', y/) = fara’, y')as + fur(x', yay 

a(x", y') = ao" + ara! + agry! + agra!® + ayra’y! + asry” 
H aghae!® + ayra!2y! + a gta!y!® + agry’® + arora’! 
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The edge conditions (24) become: 


(x'fao” + yf") 
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Calculations of the Current Coefficients a and b 


Eqs. (17)-(20), (26), and (27) determine now a 
total of 40 coefficients a” and 6”, if in f(«’y’) terms up to 
the fourth power in x’ and y’ are considered. Detailed 
calculations are presented in a separate report [52]. 
Here a short outline of the procedure is given: the 
coefficients Gm” and bm” are found in terms of V»”, 
W,,.”, and X,” defined in (17) and (18) 


In” = in Vig Woe’, DG) 


bn” = Dak m!”y Wa Gri) (28) 


The prime on the subscripts m in the brackets indi- 
cate that terms with a different subscript than m are 
involved. 

Our next problem consists in expressing Vn”, Wm", 
and X,,” in terms of the field quantities S,,.”, Tm", and 
Un". If we use the partial vector potentials given in 
(11), (12), and (16) and keep only the terms with A” on 
the left-hand side, we obtain together with (17), (20) 
and (26) 
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where only current coefficients of (7—2) and lower order 
are involved. Putting these expressions back in (28) and 
replacing all the current coefficients of lower order by 
an iterative process we finally obtain an expression 
which involves only the field coefficients Sn”, Tm”, and 
Un” 
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Following (25) we define the total current coe ffi- 

cients by 

Ge hag Ak On atk Gate oe 

Bplit= ROn Nie keD. ak Ory cal mene 


am = 


om = (31) 


and similarly for the expansion coefficients of the inci- 
dent field 
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The right-hand side is obtained by using (15) and 
Taylor’s expansion theorem. 

It turns out that all current coefficients belonging to 
different powers of k” but having the same subscript m 
have the same functional dependence on Sm”, Tm”, and 
U,”. This allows us to write the total current coeffi- 
cients in terms of the primary field directly by using 
(32). After some calculations we finally obtain 
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The superscript 2 for the incident field has been 
omitted. All magnetic field components and_ their 
derivatives should, of course, be evaluated at the center 
of the disk. 

Eqs. (33) and (34) represent the solution of our prob- 
lem from which we shall derive the results in the fol- 
lowing sections. 
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Ill. InpucED ELectRIc-DIPOLE MoMENT 


The induced electric-dipole moment is defined by 


P = { po(p)ds (35) 


where 


¢= carey (36) 
® 

is the electric-surface-charge density and p the radius 

vector on the disk. 

Obviously only odd-power terms in x’ and y’ for 
a(x’, y’) (or even-power terms for J(x’, y’)) contribute 
Lou. 

Using the results in Section II and Maxwell’s equa- 
tion we obtain 


16 ier (ka)? (132 : 3 Olas 
P,= ate £3 30 x“ Re Be 
hy Oleh 87 
ey ) 2) (ka)*E ] 
We Or 0 
16 Cae SP Geary 
Ie = a «| Bat + (1s igor ey Be az? 
De Olek 87 : 
wanes ) ape (ha) "Fy | eee 
weg OX Or 0 


The subscript 0 at the right-hand side of the main 
bracket in (37) indicates that all terms have to be 
evaluated for x=y=z=0. 

The first terms in (37) represent the well-known first- 
order approximations calculated by Bethe [21]. 

Eq. (37) shows also that in the first-order approxima- 
tion P is in the direction of the tangential electric field 
at the center of the disk. The higher-order terms, how- 
ever, lead to cross polarization, due to the terms 
0H,'/dy and 0H,'/dx. 

For the case of an incident plane wave the electric- 
dipole moment depends also on the angle of incidence 
6; and is given by (see Fig. 1) 


IP sate[ 1 +( = sin 6) (ka) 
Zea? i 7a i) (Ra) 
a) .; 
(he) | Ex'(0, 0, 0) 
Or 


16 Sonata 
Py = Foe E + (= — % sin? ) (ka)? 


8 
— j= (a) | E,(0,0,0). (38) 
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Fig. 1. 


The case for normal incidence has been treated by 
Bouwkamp [23] up to the fifth-order approximation. 
From his current distribution, one obtains 


P= Sateo[ 1+ = (ha)! — j= (bo)* + (bas 
= — — (ka)? — 7 — — 
ries ie fog seu Osa 


SLO el coro 39 
=i (eo) | (0,0,0). (39) 


IV. InpUcCED MaGnetic-D1PoLE MOMENT 


The magnetic-dipole moment is defined by 


1 
M = —| » X j(p)dS. (40) 
2/p 
In the case of a plane disk M has only a component 
M, along the disk axis. It is easily seen that only odd 
power terms in J(x’, y’) contribute. 
Again we can express M, in terms of the primary 


field 
8a’ 1 1 0°H,' 
M, = — —4d.t — — (ka)?| 30 + — 
YE = ag a) [sent SE] 
. 4 . 
+5 (aya . (41) 
Or 0 


The bracketed term is evaluated at x=y=z=0. 


For a plane wave incident at an angle 6; (see Fig. 1), 
(41) becomes 


8 1 
eae seh — 75 (2 + sin? 0) (ka)? 


4 
els. (ka) ‘ H,*(0;0,.0).. (42) 


For normal incidence the induced magnetic-dipole 
moment vanishes. 
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V. Far-Zone FIELps 


Given the current distribution we can find the scat- 
tered field in a straightforward manner from the vector 
potential as given in (1). Unfortunately a general solu- 
tion of this integral is not available. It is, however, 
possible to evaluate it for field points which are at dis- 


tances large compared with the disk radius. In this case 
we can write 


Ho al iL 


as’ > 
R—- 0 4r R 


A= 


At D Yr 


. il J exp [jkp’ sind cos (¢ — ¢’) |dS’ (43) 
D 


where R, #, ¢ are the spherical coordinates of the field 
point, p’ and ¢’ are the integration variables as shown 
in Fig. 1. Expanding the exponential in powers of k up 
to the term k* and using (8) one obtains an expression 
for A correct to the third order in k. We express the 
field in spherical coordinates 


A, = A,cosgsin’d + A,sing sin? 
As = A,cosgcosd + A,sin g cos? 
Ag — 


— A,sing+ A,cos¢. (44) 


It turns out that for large distances the field com- 


ponents become 


Hs = 


€o\ 1/2 jk 
He = |— Es = — — Az. 
Ho Ho 


After performing the straightforward but tedious 
computations one finally obtains in terms of the inci- 
dent field 


e-ikR 
ka)?a 
rR uo} \ 


(45) 
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1 4 
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2 €0 aster 2 ~~) 
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— (-s tt ts fe Sin g 
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x 


is 
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wur0 


Eggimann: Higher-Order Evaluation of Electromagnetic Diffraction by Circular Disks 


415 
eo ikk 
Hy, = (ka)?a — EA Cos P,sin 

7R ures) a8 [ ¢ + P,sin ¢] 

(ka)?— 7 OF." One 
= | sino] 6 = 8=——— cose 

45 Lwpo Oz oy 

Of | OE OE Nae 
= 8 sin? g + 8{- + ) sin ¢ cos 6 | 

Ox oy Ox 


/ 


Ces 
= 64/—sin'o[—E,s COS ~ 
Ko 


— ho sin el} cos J. 
0 


Again all the fields are evaluated for x =y=z=0. 
For a plane wave at oblique incidence we obtain 


(46b) 


e ikk 


D 
Hei = ka)*?aH) <—sin ¢; sin 6; sin? 
3 R (ka) “4 3 oy) 


4 
+ = cos ¢; cos 6; Sin g — oy sin $; cos @ 


+ (ka)? - 2 — 6sin? @,) cos ¢; cos 6; sin g 
— (32 — 10sin?@,) sin ¢; cos g 

+ sin (s cos ¢; cos 9; sin 6; sin g cos g 
+ — (2 — 3 sin? 6; — 8 cos? g) sin ¢; sin ) 
+ sin? d (- 10 cos ¢; cos 6; sin g 

+ (10+ 4 sin? 6;) sin ¢; cos °) 


— 3sin*? sin ¢; sin | 


Ca gkR 


(47a) 


4 
Hy, = (ka)*aH cos 0 iz (cos $; Cos 8; Cos ¢ 


TT 
+ sin ¢; sing) 
1 

+ (ka)? ie | @ — 6sin? 6;) cos ¢; cos @; cos g 

+ (32 — 10 sin? 9;) sin ¢; sin ¢ 

+ sin 9( (6 — 8 sin? ¢) cos ¢; cos 8; sin 9; 

+ 8sin ¢;sin 0; cos ¢ sin °) 

— 6sin20(cos ¢; cos 6; cos ¢ + sin ¢; sin | . (A47b) 
Here ¢; and 6; are the angle of incidence as shown in 
Fig. 1. These results do not seem to agree with results 
obtained by Stevenson [24]. They lead, however, to the 


same expression for the scattering coefficient as ob- 
tained by Lur’e [19] and Kuritsyn [20]. 
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VI. SCATTERING COEFFICIENT FOR INCIDENT 
PLANE WAVE 
The scattering coefficient 7 is defined as the ratio of 
the energy scattered from the disk to the energy inci- 
dent on the disk. The scattered energy is found by inte- 
grating the Poynting vector of the scattered field over 


the sphere at infinity. 
Using (47) one finds for an incident plane wave 


[of Rec we 


ra*H,? cos 6; 


24 | Hd |) sinddddy 


T= 


128 : cane 
= — (fa)? 1-4 sin 9; | — sin’ ¢; — 1 
271? cos 0; 4 
bay? 
ae ue) [(22 — 5 sin? 4,) cos? 6; cos? ; 


+ 4(88 — 54sin?6; — 5 sin‘ 6,) sin® oi) : (48) 


For normal incidence Bouwkamp [23] obtained 


6 [145 (ka)? 
2 — (ka 
Qin? x D5 


eRe eae | (49) 
18375 : eal 


A general expression for 7 for arbitrary primary fields 
would involve the separation of the real and imaginary 
parts of the incident fields and their derivatives at the 
center of the disk. The resulting expression would be 
very cumbersome and of little practical value. 


VII. DIFFRACTION BY A CIRCULAR APERTURE 


Using the results of the generalized Babinet’s princi- 
ple [51] it can be shown that the disk problem and the 
aperture problem are equivalent if the following sub- 
stitutions are made 


substitute Hell Amontune for eo E Disk 
SUDstiiute=—9 en Ea nerture for woH'pisx 
substitute Holi ectane for Pipe 


substitute Pa perture for woMpisk. (50) 
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It is customary to designate as primary field 2°, H® 
the field which would exist if the aperture is replaced by 
a solid screen. We obtain thus 

Ey oe pee 


Atan +s stelios (51) 


This leads to the following expressions for the in- 
duced electric and magnetic dipole moments: 


* eofns — = (bo| ane += | 
PP, == US = (er 3 — 
oe a4 10 Raz? 
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‘ ¥ 30 vm peas 
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| one (ka)*H9} ~ 163) 
wig Ox Or 0 

The bracketed terms are evaluated at the center of 
the aperture. 

The scattering coefficient ¢ of an aperture is usually 
defined as the ratio between the energy incident on the 
aperture and the energy transmitted through the aper- 
ture. We obtain thus 


i= 


bole 


is (54) 
where 7 is the scattering coefficient of a disk given by 
(48). 
APPENDIX 
EVALUATION OF INTEGRALS 


In the foregoing discussion we need to evaluate some 
integrals of the form 


G(x, 9) 


os fh fe’, ¥) dS’ (55 
Disk T(a?— x!? —y!?) 1/2((¢— 321)? (y—y/)?) 1/2 
where f(x’, y’) is a polynomial in x’ and y’. Bouwkamp 


[50] has given a solution for this kind of integral in the 
form: 


2m , 
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TABLE I TABLE II 
G(x, y) = IL 1 F(x’, y')p'dp'de’ = 
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where 
p M2 a yl Zt vy!2 <q? 
pP=x*?+y"* <a? 


0<o<2r 

P3 =associated Legendre functions 

n,m 7 and u=integers subject to the condition u>0, 
O0<m<n. 


The coefficients A, are: 


A,(n, m, p) 
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Low-Noise Properties of Microwave Backward Diodes’ 


SVERRE T. ENG{, MEMBER, IRE 


Summary—This paper describes, for what is believed to be the 
first time, the low-noise properties of backward tunnel diodes in 
microwave applications. The physics of the diodes are reviewed to- 
gether with some of the characteristics and equivalent circuit param- 
eters. The diodes are then considered as mixer diodes with IF in 
the audio range and also the standard 30-Mc IF. Another promising 


application considered is the use of the backward diodes in low-level 
detection. 


The results show that the noise figure at 13.5 kMc with a 1-kc IF 
is around 15 db better than any commercially available mixer 
diodes. Using 30-Mc IF, the noise figure of backward diode mixers 
is without special optimum design, comparable to the best mixer 
diodes on the market. Of great importance, especially in micro- 
miniaturization, is the fact that these diodes may be used with a 
very low local oscillator power (50 uw or less). The high nonlinearity 
of the I-V characteristic at the origin and the low 1/f noise properties 
of these diodes are also of benefit in crystal video receivers and other 
low-level detector applications. 


I. INTRODUCTION 
S | XN HERE are three basic types of noise which occur 


in semiconductor devices. Thermal noise is caused 

by random motion of electrons. When a drift 
velocity is superimposed by means of an electric field, 
we get shot noise. Then we have 1/f noise which is de- 
tected over and above thermal and shot noise, and which 
is distinguished by its spectral intensity. 

The investigation reported in this paper is con- 
cerned with low 1/f noise (excess noise) and thermal 
noise in experimental samples of microwave backward 
diodes. The problem of noise reductions in mixer diodes 
is discussed in Section IJ, and the ideas which led to the 
fabrication of the backward diodes are also described. 

Backward diodes made by diffusion or alloying tech- 
niques have not yet been accessible for microwave 
applications because of the large capacitances associ- 
ated with the junctions. However, the pulse-bonding 
technique used in our experiments proves to be well 
adapted for microwave versions of these diodes. The 
physics of the diodes are also reviewed together with 
some of the characteristics and equivalent circuit pa- 
rameters. These diodes have a nearly conventional 
diode forward I-V characteristic and the back char- 
acteristic of a tunnel diode. 

A reduction of 1/f noise is of interest, for example, in 
Doppler radar systems where the IF is in the audio 
range. Since light weight, small size and simplicity are 
very important in this type of antenna-navigations 
radar, it does not seem desirable in some cases to 1n- 
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crease the complexity of the system and go to higher 
IF (30 Mc) and thus possibly eliminate the 1/f noise 
problems. In some mixer diodes there is indication of a 
slight contribution of excess noise at 30 Mc. However, 
the main contribution to the noise at this frequency is 
shot and thermal noise. Because of the low spreading 
resistance in the semiconductor wafer, the thermal 
noise will be reduced somewhat compared with that of 
ordinary mixer diodes. Thus, an improvement of the 
sensitivity of systems using 30-Mc IF may also be ex- 
pected. Present video crystal receivers and other low- 
level detector applications may also benefit in sensi- 
tivity from an 1/f noise reduction in diodes. All these 
application aspects of the backward diodes are de- 
scribed in the last part of this paper. 


Il. REpucTION oF DIopDE NOISE 


Although the origin of the 1/f noise in semiconductors 
is not precisely known, a variety of experiments have 
indicated that the surface plays a significant part in 
connection with this noise. It has been observed that 
1/f noise increases with higher density of slow surface 
states, and, for a given density of slow states, the 1/f 
noise will be higher for an inversion layer than for an 
accumulation layer [1]. 

One approach for decreasing the 1/f noise is a sta- 
bilization of the surface by etching, ambient gas con- 
trol, or thermally grown oxide layers. These techniques 
for stabilization and evaluation have been successfully 
applied to conventional mixer diodes in our laboratory. 

However, the author is suggesting another approach 
to the problem of 1/f noise reduction in mixer diodes. 
For nearly degenerate material the relative change in 
surface potential is much less pronounced than in the 
nearly intrinsic case. Thus, if nonlinearities in the I-V 
characteristics of diodes made of much lower resistivity 
material than conventional mixer diodes were utilized 
for mixing purposes, the surface contribution to the 1/f 
noise should be reduced. Also, since the resistivity of 
the semiconductor wafer is reduced, the thermal noise 
contribution from the spreading resistance may be 
smaller. 

The backward diodes used in our experiment have 
doping levels one and two orders of magnitude higher 
than conventional mixer diodes. Another point of in- 
terest is that the 1/f noise decreases when the rectified 
current decreases [2]. Since the conversion loss can be 
preserved at very low local oscillator power when using 
the backward diodes as mixers, improvements in the 
over-all noise figure of systems using IF frequencies in 
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the audio range should be obtained. The reduction of 
the spreading resistance may also be of value in de- 
creasing the noise figure of receiver systems using 30- 
Mc IF. 


III. Tue Backwarp Diope [3], [4] 


Since electrons can be considered to have wave prop- 
erties, they have the ability to penetrate potential bar- 
riers which would be impossible for particles according 
to classic theory. This effect is known as quantum- 
mechanical tunneling. A backward diode is a special 
case of a diode exhibiting quantum-mechanical tunnel- 
ing phenomena. The most familiar example is the tunnel 
diode, where the impurity concentrations on both sides 
of the p-” junction are such that an increase in tunneling 
current is followed by a decrease, producing a negative 
resistance. The expression for the tunneling current is 


(1) 


| 
F 5) 


J exp | - 


where 


A =numerical constant 
m* =reduced mass of holes and electrons 
Eg=band gap of semiconductor 
F=average electric field across the space-charge 
region. 


F is inversely proportional to the width of the space- 
charge region of the p-m junction. A wide space-charge 
region means low tunneling current; a narrow one 
means high tunneling current. For a backward diode, 
it is necessary to reduce the tunneling current from the 
high levels of a tunnel diode to such levels that a nega- 
tive resistance is either nonexistent or is so high that it 
can be disregarded. The expression for the width of the 
space-charge region of an abrupt junction is 


2K Ko(N N 
w =| (Np + Na) 


eNpN a 


iL 
Wo+Va| (2) 


where 


Np =donor concentration 
Na =acceptor concentration 
Vp=built-in “diffusion” voltage 
V4 =applied bias voltage 
K Ko=dielectric constant of semiconductor 
e=electronic charge. 


If, in fabrication, one doping level is kept constant at 
soine high level, the space-charge region can be widened 
by simply decreasing the impurity level of the semi- 
conductor wafer. This causes a decrease of the tunneling 
current, and enables the device designer to adjust it to 
the proper level. 

It is theoretically possible to make backward diode 
junctions by diffusion techniques. However, it is very 
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impractical to control a diffusion of only several hun- 
dred atomic layers into the wafer. The technological 
difficulties of making electrical connections to such a 
shallow diffused layer are also severe. In addition, dif- 
fused junctions are undesirable because the tunneling 
probability decreases exponentially with the width of 
the space-charge region. 

Tunnel diodes and also backward diodes are gen- 
erally made by alloying techniques. An appropriate al- 
loy pellet is placed in contact with a semiconductor 
wafer and heated to such a temperature that wetting 
and dissolutions of some of the semiconductor into the 
alloy occurs. 

The alloy contains the impurity that will give opposite 
conductivity type to that in the wafer. On cooling the 
dissolved semiconductor will epitaxially grow on the 
base semiconductor. This regrowth will now be doped 
to the opposite conductivity type. There will be an 
abrupt transition between the doping levels on the two 
sides of the junction. The capacitance of this type of 
junction may be too large for practical use in the micro- 
wave region. However, by subsequent etching, the junc- 
tion diameter can be reduced. The main advantage of 
this technique is that fairly good control of the current- 
voltage characteristic can be obtained. The disadvan- 
tages are that capacitances in the 0.2 wuf range are diffi- 
cult to reproduce and that the structure is mechanically 
weak. 

The pulse-bonding technique seemed to be the best 
way of fabricating our diodes. A gallium-plated gold 
wire was brought into contact with an n-type germa- 
nium semiconductor wafer with a resistivity in the range 
of 0.001 to 0.004 ohm/cm, and the junction was formed 
by passing a pulse through it. The pulsing equipment 
has certain parameters to vary, such as the heights and 
lengths of the pulse. The I-V characteristic of the 
diode can be displayed on a scope so that the junction 
can be pulsed to the desirable current-voltage char- 
acteristic. The capacitances obtained are well below one 
picofarad. The capacitance-voltage relationship is that 
of a step junction. 

Unlike tunnel diodes, the package is not too critical. 
A standard computer glass package supplemented with 
gold-plated brass adaptors and with over-all dimensions 
confirming with ordinary mixer diode packages was 
found to be satisfactory at 13.5 kMc. 

The equivalent circuit of the diode used in the experi- 
ment is shown in Fig. 1. C; is the space-charge capaci- 
tance, which was around 0.2-0.5 uwuf at zero bias. The 
series wafer resistance Rs was found to be approx- 
imately 10 ohms when estimated from the shape of 
the I-V curve in the forward direction. The lead in- 
ductance Ls is muh when measured at 1 kMc, and C 
is the stray capacitance of the package. The variable 
resistance Ry varies between 300 ohms and 80 ohms for 
local oscillator powers in the range from 10 pw to 1 mw. 
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The nonlinear barrier resistance is believed to be inde- 
pendent of frequency from dc to well beyond the 
microwave region and is also considerably less tempera- 
ture dependent than that of a conventional mixer 
diode. 

The backward diode I-V characteristics are shown 
in Fig. 2. As the reverse bias is increased, the supply 
of electrons which are able to tunnel increases without 
limit. At zero bias no current flows through the junc- 
tion, At forward bias only a very small current is able to 
flow by the tunneling process. The excess current is, in 
some cases, fairly constant until carrier injection takes 
place giving rise to the normal forward characteristic of 
a p-n junction. 


(pramps) 
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DIODE CURRENT 
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Fig. 2—The backward diodes’ I-V characteristics. 


IV. APPLICATION AND MEASUREMENT 
A. Mixer With IF in the Audio Range 


Radars for navigation and moving target indications 
use the Doppler effect in one form or another. Since 
light weight, small size and simplicity are very im- 
portant, it does not seem desirable in some cases to in- 
crease the complexity of the system by using 30-Mc 
IF to eliminate the 1/f noise problems. One simple ar- 
rangement is shown in Fig. 3. Power from the trans- 
mitter is mixed with the CW echo from a target. The 
beats between the two frequencies f and f’ can be heard 
in phones or indicated on other types of indicators. The 
echo frequency is given by the well-known Doppler 


formula, 


(3) 
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The beat frequency is then 


(4) 


where v is the target velocity and is the wavelength of 
the transmitted signal. The beat frequency is usually 
in the ke range. 

The operating noise figure of our Doppler system can 
be written as 


> Bi 
Lon = ha a 1) B. =F IEE, ar Phy = 1); (5) 


U 


where 


fag =scuice temperature normalized with respect to 
room temperature 
B,=useful channel bandwidth 
B,=total channel bandwidth 
Ly» =mixer conversion loss 
F\y=IF amplifier noise figure. 


The properties of the crystal diode itself are involved 
explicitly or implicitly in this equation. In our system 
the total bandwidth is approximately twice the useful 
bandwidth (B;/B,2), since the signal is entering only 
one of the channels, while the antenna noise is appearing 
in the signal channel as well as the image channel. 
Since 13.5 kMc is used in several types of navigation 
Doppler radar systems, the noise figures of the diodes 
were tested at this frequency in the experimental system 
shown in Fig. 4. In this way the noise figure can easily 
be measured with the so-called “signal generator 
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Fig. 3—Schematic diagram showing a simple Doppler radar system. 
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Fig. 4—The Doppler radar noise-figure measurement system. 
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method.” The results obtained agree, within measure- 
ment errors, with calculated values of the noise figure 
obtained from separate measurement of the conversion 
loss and noise temperature. Microwave power at 13.5 
kMc travels from the signal generator to the magic tee. 
One part of the power is delivered to the test diode and 
may be thought of as the local oscillator power of the 
mixer. Another part of the power travels through the 
precision attenuator to the “Blitz Wheel” mechanical 
modulator. The modulator reflects power at 13.5 kMc 
+990 cps, and part of this signal is available for deliv- 
ery to the diode. This power may be considered the in- 
put signal, and it beats with the local oscillator sig- 
nal in the test diode to produce the difference frequency 
of 990 cps, which is then amplified together with the 
audio noise generated by the diode. The calibration of 
the system is made by careful measurement of the in- 
dividual insertion losses of the various microwave com- 
ponents and the response characteristic of the output 
meter. The signal amplitude is controlled by-the preci- 
sion attenuator in front of the modulator, and this in- 
put signal is adjusted so that the noise plus signal out- 
put power is double the output power obtained when the 
input signal is zero. Then the noise figure is obtained 
from the ratio of the input signal power used to the 
available thermal noise power from the source. 

In order to properly evaluate the capabilities of the 
diodes, it was also found necessary to measure the con- 
version loss, noise temperature and the noise tempera- 
ture spectrum of the experimental crystals [6]. 

The conversion loss L,, is defined as the ratio of the 
available RF input signal power to the measured IF 
output power at the mixer. In the “amplitude modula- 
tion method” used, the output of an oscillator at 13.5 
kMc was amplitude modulated at 990 cps. The modula- 
tion envelope, after crystal detection, was developed 
across an output load. The conversion loss was then 
found from this voltage, since the percentage modula- 
tion and the power of the oscillator were known from 
careful calibration. 

The mixer noise temperature ¢,, is defined as the ratio 
of available noise power output of the crystal to that of 


a resistor at room temperature. The “ Y-factor method” 
[6] in which 


dre eee i eal eed (6) 


was used in our case. Y is the ratio of system noise out- 
put with the test diode in the crystal holder to the 
noise output when a resistor equal to the IF dynamic 
resistance of the crystal was substituted for the crystal 
in the holder. The input unmodulated signal was at 13.5 
kMc, and the crystal was followed by a 990-cps VSWR 
(bandwidth 40 cps, /ir~6 db) and a transistor post 
amplifier. 

The noise temperature spectrum was measured under 
the same conditions as the noise temperature, except 
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that the crystal was followed by a wide-band transistor 
amplifier and a spectrum analyzer. 

In Fig. 5 the performance characteristics of a back- 
ward diode mixer with 990-cps IF are plotted. The 
noise figure of a commercially available low 1/f noise 
diode (1N1838) is also shown for comparison. First of 
all, it is indicated that the backward diode has a consid- 
erably lower noise figure which is caused by a reduction 
of 1/f noise (see Fig. 6). Another important considera- 
tion is that the diode is capable of operating with very 
low local oscillator powers without any dc bias. The rea- 
son is that the nonlinear region of the I-V characteristic 
is in the vicinity of the origin, while in ordinary mixer 
diodes the nonlinearity occurs around the contact volt- 
age, which may be 0.3 to 0.6 volt in the forward direc- 
tion. Thus, the backward diode may be used in a low- 
noise mixer using a tunnel diode or variable-capacitance 
diode harmonic generator as the local oscillator. 

A more direct comparison of the two different diodes 
is made in Table I. Since the RsC; product is higher in 
the BD No. 5 diode than in the 1N1838 diode, the con- 
version loss would be expected to be higher since the 
conversion loss is an increasing function of this product 
[6]. However, it should be pointed out that no real at- 
tempt has been made to optimize the fabrication tech- 
nique. 

Fig. 7 shows the noise temperature spectrums for the 
two diodes. It is interesting to note that the spectrums 
follow approximately the 1/f law. The turnover point 
from 1/f noise to shot and thermal noise is around 20 
ke for the backward diode, and the noise temperature 
seems to be roughtly proportional to the square of the 
rectified current in the 1/f noise range. This current 
dependence of the noise temperature spectrum corre- 
sponds with results obtained for conventional point- 
contact mixer diodes [2]. 

It can be concluded that the real improvement in the 
noise figure obtained by the backward diode compared 
with that of a conventional mixer diode is caused by 
reduction in noise temperature. There is reason to be- 
lieve that 


1) 1/f noise is less pronounced in the backward di- 
odes since they are made of lower resistivity ma- 
terial than ordinary mixer diodes. 


2) 1/f noise is lower at the optimum noise figure be- 
cause less local oscillator power (and thus recti- 
fied diode current) is needed for satisfactory 
operating performance. 


B. Mixer with 30-Mc IF 


In optimum receiver design it has been found that 30 
Mc is a good choice for the intermediate frequency us- 
ing present commercially available mixer diodes. This 
optimum IF is mainly determined by the noise figure of 
the IF amplifier and the mixer noise temperature. The 
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990-cps IF. The noise figure of a 1N1838 diode (low 1/f noise) 
is also shown for comparison. 
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Fig. 6—The noise temperature spectrums for a backward 
diode and a 1N1838 low 1/f noise Doppler mixer diode. 
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reduction of 1/f noise in mixer crystals may therefore 
influence receiver design. 

At 30 Mc the main contribution to the noise tem- 
perature is shot and thermal noise. However, in some 
cases a very small excess noise contribution may be 
present. Nevertheless, since the backward diodes are 
made of lower resistivity material than conventional 
mixer diodes, a small improvement may be expected in 
the noise temperature since the thermal noise caused by 
the series bulk resistance will be reduced. 

The experimental results are shown in Fig. 7. A 
comparison is made with the 1N23WE mixer diode. The 
noise figure was calculated according to (5) when f,=1, 
and the conversion loss was measured as outlined in 
Section IV-A. The noise temperature was measured 
with a 30-Mc IF amplifier (7tp=1.7 db) following the 
crystal. 

In Table II the comparison is made at 300-uw local 
oscillator power, which seems to give the best noise 
figure for both diode types, and 30-ww, which may be 
of interest in lightweight systems since this small 
amount of power may be readily available from tunnel- 
diode oscillators or harmonic generators using variable- 
capacitance diodes. 
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the conversion loss. This was expected, since the RsC; 
product is highest for the backward diode. However, the 
1N23WE seems to have a lower nonlinearity coefficient 
for the barrier resistance than the backward diode (this 
can be verified by inspection of the I-V curves for the 
respective diodes). Also, since the backward-diode noise 
temperature is, indeed, lower, there may be reason to 
believe that, by optimizing the fabrication technique so 
that a lower value of capacitance can be obtained, the 
backward diode may be superior to the 1N23WE in this 
application with 300-uw local oscillator power. 

At lower local oscillator powers the backward diode 
has much better noise figure than the 1N23WE diode. 
The main reason for this is that the 1N23WE diode does 
not have the same ability to rectify small powers as the 
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backward diodes. Thus, the conversion loss will be sub- 
stantially higher in the 1N23WE diode at low ‘local 


oscillator powers. 


C. Low-Level Detection 


Another major application of the backward diode is 
in video detection. In a crystal video receiver, for ex- 
ample, the incoming modulated RF signal is detected 
immediately at the input from the antenna, and the re- 
sulting video signal is amplified by a high-gain video 
amplifier. The advantages of this type of receiver are 
simplicity, small size, low cost and broad RF band- 
width. The price paid for these attractive features is a 
large loss of sensitivity compared with that obtainable 
with a superheterodyne receiver. However, there are 
applications where the lower sensitivity is acceptable 
and where the small size or broad band is required 
(beacon radar receivers, square-law detectors in power 
measurement instruments, etc.). 

The SNR of a detector crystal is proportional to a 
figure of merit [6] M, which combines the rectification 
properties and the noise generation of the crystal and 
the amplifier. A higher figure of merit indicates a 
greater sensitivity and thus a lower minimum detect- 
able signal. 


BR, 


ES — 7 
V Rs + Ra 2 


where 


6 =short-circuit current sensitivity of the diode 
R,=video impedance of the diode 
KR, =equivalent noise resistance of the video amplifier. 


The current sensitivity @ was measured at 13.5 and 
6 kMc by observing the short-circuit rectified cur- 
rent 2 for a given input power P. Then 8=1/P. The 
video resistance R, is the dynamic dec resistance of the 
crystal at the operating point, and it was measured 
with a bridge at 100 ke. 

The short-circuit current is plotted in Fig. 8 as a func- 
tion of the absorbed input power at 13.5 and 6 kMc for 
a backward diode and a detector diode MA408B. The 
backward diode has a higher current sensitivity and a 
square-law response of a larger range than the conven- 
tional detector diode. A more direct comparison is made 
in Table III. 

At 13.5 kMc the backward diodes seem more ad- 
vantageous according to the figure of merit. The rea- 
son for the higher current sensitivity is partly due to 
the lower video impedance R, compared with that of the 
MA408B diode. However, the product RsC; is lower 
and the I-V nonlinearity is also higher for the backward 
diode. Both factors will result in higher current sensi- 
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Fig. 8—The short-circuit current as a function of the absorbed in- 


put power for a backward diode and a detector diode MA408B. 


TABLE III 
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BD NO. 5 10 (ee) 730 Dik 45 iii 180 
BD NO. 7 8 0.5 880 2.4 47 (pPcse |) eS 


tivity. At 6 kMc the MA408B shows a higher figure of 
merit which may be contributed by the higher video 
resistance. There are good reasons for believing that the 
video resistance of the backward diode may be increased 
by suitable optimization of the shape of the I-V curve. 

Other possible advantages of the backward detector 
diodes compared with conventional detector diodes 
may be the expected smaller variations in the per- 
formance characteristics with variations in tempera- 
ture and nuclear radiation. 


V. CONCLUSION 


Experimental microwave backward diodes have been 
made, and the low-noise properties of these diodes in- 
vestigated and compared with the best commercially 
available crystals. The performance characteristics 
show that backward diode mixers are of considerable 
advantage in receiver systems where the IF is in the 
audio range. The improvement is mainly caused by a 
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reduction in 1/f noise. It should be noted that a back- 
ward diode mixer can be operated with an order of mag- 
nitude lower local oscillator power than mixer diodes 
used today. Very satisfactory performance was also 
obtained using 30-Mc IF and in low-level detection, in 
spite of the fact that no real attempt was made to 
optimize the diode fabrication technique. In fact, 
there are good theoretical reasons for believing that 
other wafer materials may well prove to be superior to 
germanium for the applications considered in the 
paper. Since the backward diodes are virtually inde- 
pendent of the lifetime of minority carriers or of the sur- 
face treatment, they can tolerate larger doses of nuclear 
radiation than conventional mixer diodes. The tunnel- 
ing portion of the I-V curve is also expected to be sub- 
stantially independent of temperature. 
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Design Theory of Up-Converters for Use as 
Electronically- Tunable Filters* 


GEORGE L. MATTHAEI}, MEMBER, IRE 


Summary—The up-converters discussed use a single diode, a 
wide-band impedance matching filter at their signal input, a moder- 
ately wide-band impedance matching filter at their pump input, and 
a narrow-band filter at their sideband output. With a narrow-band 
filter at the sideband output, the frequency which will be accepted by 
the amplifier can be controlled by varying the pump frequency. Anal- 
ysis of the impedance matching problem involved shows that tuning 
ranges of the order of a half-octave to an octave are possible. Theory 
is presented for both the lower-sideband and upper-sideband types 
of tunable up-converters and for the design of the required imped- 
ance-matching networks. It is shown that, because of the pump input 
bandwidth required, it will generally be necessary to accept some 
mismatch at the pump input. But, by use of a properly designed im- 
pedance-matching filter, the reflection loss can be kept nearly con- 
stant across the pump band, and the incident pump power required 
is not unreasonable. It is seen that properly designed devices of this 
type using voltage-tunable pump oscillators should have wide tuning 
range, fast tuning capability, a useful amount of gain, no image re- 
sponse, and a low noise figure. 


* Received by the PGMTT, April 17, 1961; revised manuscript 
received, June 19, 1961. The research reported here was sponsored by 
the U. S. Army Signal Res. and Dev. Lab., Ft. Monmouth, N. J., 
under Contract DA 36-039 SC-74862. 

+ Stanford Res. Inst., Menlo Park, Calif. 


I. INTRODUCTION 
A. Description of the Proposed Devices 


REVIOUS work dealt with the application of 
filter theory to the design of wide-band para- 
metric amplifiers and up-converters.' The present 
discussion applies a similar theoretical approach to a 
different but closely related problem. The objective will 
be to obtain an electronically controlled wide tuning 
range using up-converters having a wide-band input- 
impedance-matching filter, a narrow-band output- 
impedance-matching filter, and a voltage-tunable pump 
oscillator such as a backward-wave oscillator. 
Defining f as the input frequency, f’ as the sideband 
output frequency, and f? as the pump frequency, for a 
lower-sideband up-converter the output is at the lower- 


1G. L. Matthaei, “A study of the optimum design of wide-band 
parametric amplifiers and up-converters,” IRE TRANs. on MIcrRo- 
WAVE THEORY AND TECHNIQUES, vol. MTT-9, pp. 23-28; January, 
1961. 
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sideband frequency 
f= ff. (1) 


For an upper-sideband up-converter, the output is at 
the upper-sideband frequency 


ame (2) 


For either type of up-converter, tuning action can be 
achieved if a narrow-band filter is used at the output so 
that only frequencies equal or very nearly equal to a 
specific frequency fo’ can be passed. If the pump fre- 
quency f? is varied, then the input frequency which will 
be accepted by the amplifier will be given by 


Pe Sh (3) 
for the case of lower-sideband up-converters, and 
J ged (4) 


for the case of upper-sideband up-converters. In both 
cases, the amplifiers will yield gain. The lower-sideband 
type introduces some negative-resistance amplifica- 
tion in addition to the up-converter amplification,?? 
however, and will therefore generally give more gain. 
For gain to be achieved, the variable-capacitance diode 
must see proper terminations at both frequencies, f and 
fo’. Since fo’ is a fixed frequency, it should be relatively 
easy to maintain proper termination at that frequency. 
The input frequency f varies, however, and the tuning 
range of the amplifier will be determined largely by the 
range of f for which proper terminating conditions can 
be maintained. 


B. Factors Permitting Large Tuning Range 


For convenience, let us first consider the case of an 
upper-sideband up-converter with a fixed pump fre- 
quency. We analyze a wide-band, upper-sideband up- 
converter using a fixed pump frequency f?, with an in- 
put-impedance-matching filter whose pass band is cen- 
tered at fo and an output-impedance-matching filter 
whose pass band is centered at fo’. In order to obtain 
maximum operating bandwidth, it is necessary that the 
condition 

Ww a 

neh is (5) 

WwW fo 
be satisfied, where w is the fractional bandwidth Af/f, 
of the input filter and w’ is the fractional bandwidth 
Af’/fo’ of the output filter.! Eq. (5) says, in effect, that 
the input and output bands Af and Af’, respectively, 
must be equal. From the theory of wide-band up-con- 


* J. M. Manley and H. E. Rowe, “Some general properties of 
nonlinear elements—part 1. General energy relations,” Proc. IRE 
vol. 44, pp. 904-913; July, 1956. 

° H. E. Rowe, “Some general properties of nonlinear elements. IT. 
Small signal theory,” Proc. IRE, vol. 46, pp. 850-860; May, 1958. 
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verters,! it is seen that if, for example, Af were made 
larger than Af’, the output filter bandwidth Af’ would 
have to be shrunk to be smaller than it could be if Af 
and Af’ were made to be equal. Since the smaller band- 
width will be the one that limits the over-all transmis- 
sion, the operating bandwidth necessarily becomes 
smaller if Af and Af’ are made unequal. Thus it is seen 
that for the usual case where the pump frequency is 
fixed, wide-band impedance matching must be ac- 
complished both at the signal input and also at the 
upper-sideband output band, and best results are ob- 
tained if the bandwidth for impedance match is equal 
for both of these bands. 

For the tunable up-converters considered herein, a 
quite different situation exists. In this case, narrow 
rather than wide bandwidth is desired. A wide-band 
impedance match is desired at the input channel, how- 
ever, so that the input acceptance band can be moved 
about by varying the pump frequency. Under these 
conditions, we may expand the input bandwidth Af and 
be quite happy with the required shrinkage of the side- 
band output bandwidth Af’. With such an amplifier, the 
bandwidth Af’ of the upper-sideband filter will be the 
operating bandwidth of the amplifier, while Af will be- | 
come the tuning range. If there were no practical con- 
siderations, such as the effects of losses and the need 
for practical impedance levels, the tuning range Af 
could be made to approach infinity while the output 
bandwidth Af’ would approach zero; practical consid- 
erations, however, do limit the obtainable tuning range. 
Nevertheless, it is clear that the tuning range obtainable 
with such a device is considerably greater than the 
bandwidth obtainable in an up-converter using a fixed 
pump frequency. As is seen in Matthaei,! bandwidths 
of the order of 16 per cent appear to be feasible for 
upper-sideband up-converters having fixed pump fre- 
quencies. The tuning range which is feasible for up-con- 
verters with varying pump frequencies is thus likely to 
be several times this value. 

An additional problem, not yet mentioned, arises in this 
case of an up-converter with varying pump frequency. 
Additional impedance matching for the pump input 
channel must be provided in order to pump the diode 
satisfactorily at frequencies throughout the required 
pump frequency range. It can be shown that it will gen- 
erally be impossible to obtain a good match between 
the pump signal oscillator and the small resistance in 
the diode while tuning across the required pump fre- 
quency band. However, it is possible to design impe- 
dance-matching filters which will give a nearly uniform, 
minimum reflection level throughout the required band. 
Then it is merely necessary to have a pump oscillator 
with sufficient available power to compensate for the 
reflection loss. 

The discussion above has been phrased in terms of 
upper-sideband up-converters. The same principles, 
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however, apply to lower-sideband up-converters. The 
main difference between the two devices is that the 
lower-sideband type tends to reflect a negative resist- 
ance through the coupling action of the varying ca- 
_pacitance. This negative resistance has the advantage 
that it contributes to the gain of the amplifier, but it 
also has the disadvantage that it could lead to oscilla- 
tion if the amplifier were improperly terminated. Calcu- 
lations indicate that good results should be obtainable 
by keeping the negative resistance small, so that stabil- 
ity will not be a problem. 


II. DEFINITION oF. DIODE PARAMETERS 


The specification of the diode parameters will be 
the same as that used previously,!? but will be sum- 
marized here for easy reference. The diode capacitance 
parameters Cy and C; are defined by 


C(t) = Co + 2Ci cos (2af?t + 1) + ++ -, (6) 


where C(t) is the time variation of the pumped diode. 
If the diode is to be resonated in shunt, it is convenient 
to use the short-circuit admittances!3 


Bu = InfCG, Bia IGi 
Bay Oat’ Cx, Beta f'Co. (7) 
In the discussion to follow, subscript zeros, such as 
(Bir)o; (Biz)o, (Bai)o, (B22)o, 


will be used to indicate parameters evaluated at the 
signal input midband frequency fo, or at the sideband 
output midband frequency fo’, whichever is appropriate 
for the given parameter. 

The equivalent shunt diode loss conductances are 


defined as 
XY 2rfoCo a (Bii)o 


a= (8) 
Qa Qa 
seen by input signal components and 
2afo'C Boo \2 
f= TE Et (2) 0) 
Qa Qa fo 


seen by the sideband components. Qu is the operating Q 
of the diode at frequency fo, and Qa’ is the operating Q 
at fo’. Since the diode loss resistance is effectively in 
series with the capacitance, 
fo 
Qa’ = Qa— 
fo’ 
Analogously, if the diode is to be resonated in series, 
it is convenient to use the open-circuit impedances* 


(10) 


1 yy iby 
a= dare 
nce S Daf’Ca* 
1 1 
X94 Xo (11) 


i InfC?  Wnf’Co® 
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where 
Co(1 — a? Co 
Co® = Co(1 — a), Cy = a = (12) 
a a 
and 
Cy 
a=—- (13) 
Co 


In this case, the effects of losses are represented by series 
resistors 


i 1 se (Xi)o _ (X22) 0 
Qa2tfoCo® Qa Qa’ 


which are the same at both the signal and sideband fre- 
quencies. 


Kesh, 


ees) 


III. DEFINITION oF FILTER PARAMETERS 


In order to reduce the number of degrees of freedom 
involved in the amplifier design, the design of the vari- 
ous band-pass filters required will be based on the use of 
lumped-element, low-pass prototypes. When these low- 
pass prototypes have been specified, the only parame- 
ters which remain to be specified in the corresponding 
band-pass filters are their center frequencies, their 
impedance levels, and their fractional bandwidths. 
Their center frequencies, of course, will be determined 
by the desired operating bands. The impedance level of 
the input and pump channel filters will be determined 
from considerations of their desired bandwidths and 
the impedance-matching properties which are required. 

Fig. 1 shows a typical low-pass prototype and a 
typical Tchebycheff response for such a filter. Tables of 
normalized element values for Tchebycheff and maxi- 
mally flat low-pass filters of this kind are presently 
availablet> Fig. 2 shows the corresponding band-pass 
filter obtained by the standard low-pass to band-pass 
transformation given in Fig. 2(a). 

In an actual microwave device it will usually be 
necessary to use resonators which are a combination 
of semi-lumped and transmission-line elements rather 
than lumped-element resonators such as are shown in 
Fig. 2. For this reason, it is convenient to use the equiv- 
alent representation shown in Fig. 3. In this case, the 
actual nature of the resonators is left unspecified, but it 
is assumed that in the vicinity of fo they exhibit reso- 
nance characteristics similar to those of the resonators in 
Fig. 2 [as is implied by Fig. 3(a) and (d)]. Then the 
properties of the individual resonators are specified in 
terms of the resonant frequency fo and a slope parame- 


4, Weinberg, “Additional tables for design of optimum ladder 
networks, Pts. 1 and 2,” J. Franklin Inst., vol. 264, pp. 7-23, 127— 
138; July, August, 1957. - 

6 W. J. Getsinger, et al., “Research on Design Criteria for Micro- 
wave Filters,” Stanford Res. Inst., Menlo Park, Calif., Tech. Rept. 
No. 4, SRI Project 2326, Contract DA 36-039 SC-74862; 1958. 


Matthaei: Design Theory of Up-Converters for Use as Electronically-Tunable Filters 


428 IRE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES seprempber 


Ae a [Ps 
E6® T? te, aA. 
Aug, — ohms RESISTANCE Ln 


|; ——henries INDUCTANCE Or eH 
€ —— farads CAPACITANCE --- 


(a) 
a 
7 WHERE 
AMETER) vee ae ) (b) 
2 z a ne — 
5 x, #(SLOPE PAR Fale eee 
S 
Zz AND = 
WwW R,. = 
c 
E bo oe (e)) 
FOR SHUNT BRANCHES, f NEAR f° 
0) 2 
Q, RADIAN PREGuENCY A-2550-5IRA B se b (d) 
, RADI k Ff i i 
Fig. 1—Low-pass prototype filter and a typical Tchebycheff response. WHERE 
dB | 
3, (SLOPE PARAMETER)=—>— =—— (A, 2¢,) Ce) 
kh 2 dwiw-w, wR, Cuz 
Fig. 3—General description of band-pass filters in terms of resonator 
slope parameters (%, J;, c, and rz are defined in Fig. 1). 
aes ( i ‘.) = ; 
ae an UNS oe a 
w I, fi 
WHERE SOUTPUT 
——= ene 
5 a Vee (b) a 7} 
| cc |_| Resa 
DIO 
2 DIODE : 
5 cig esis. (c) Ro GR 1: aR 
= = Tr} Hn eee AI pens 
s jt | alle C8 La a | 
2 i iS) 8, {E > E+} 8) if 
“ , ( cule (eee | [ Z 
i = 2 BANDWIDTH Sheers z. ape eis PUME 
0 fide mace INPUT FILTER 
; | fo) 2 


PUMP FILTER 
f, FREQUENCY 


Fig. 4—A possible circuit for an up-converter for 


Sine ees electronically-tunable filter applications. 


(d) 


COUPLING BOX TO 

REPRESENT COUPLING 

BETWEEN ENERGY 

AT FREQUENCIES f 

AND f’ DUE TO CF SIDEBAND 
OUTPUT 


oa 


R, IS CHOSEN TO GIVE DESIRED. INPUT FILTER 
IMPEDANCE LEVEL, AND R,=R, (4,/,) 


SIDEBAND 
RESONATOR 


FOR SERIES BRANCHES 


Petes ears ata ae ge 


L L 9) 
IGNA , 
ee 202 
FOR SHUNT BRANCHES . Z, Zo Zo 
-+——-—-—-— nn eee a) 

( \ wR, | (£) AT INPUT FREQUENCY f AT SIDEBAND 
o2(— ime a, y2( )( ) : 
k ae yall hk Wo 2,6, %, FREQUENCY f 


Fig. S—A circuit which is approximately equivalent to that in Fig. 4 
Fig. 2—Summary of relations for design of lumped-element for energy components at the input f etal 
band-pass filters from low-pass prototypes. equanes Hike cae eam tatabccinest a1 ook 


~ ¥ 


1961 Matthaei: 


ter. For a series, lumped-element resonator, the reactance 
slope parameter is simply 
1 WO dX; 
woC; 2 dw #9 ; 


0 


(15) 


where wo = 2mf and X; is the reactance of the resonator 
circuit. The derivative form in (15) and in Fig. 3(b) 
gives the corresponding reactance slope parameter for 
any series-resonant circuit, whether of lumped-element 
form or not. As indicated in Fig. 3(e), for circuits ex- 
hibiting a shunt resonance an analogous susceptance slope 
parameter b;, applies. Fig. 3(b) and (e) relates the 
band-pass filter resonator slope parameters to the ele- 
ments of the low-pass prototype. 

The parameter Ryo defined in Fig. 3(c) equals Rz for 
the circuit as shown. However, the definition given for 
Ryo is introduced because some microwave filters give 
impedance-level transformations so that R,) may not 


necessarily equal Rr. 


IV. AN Up-CoNVERTER MODEL FOR PURPOSES 
OF ANALYSIS AND D1IscuUSSION 


The up-converter model discussed in this section in- 
corporates some features which have been used at SRI 
quite successfully on a degenerate parametric amplifier, 
and it is believed that they could also be used to ad- 
vantage in electronically tunable up-converters in cer- 
tain frequency ranges. However, the design equations 
and charts in the following sections should be applicable 
for any of a variety of coaxial, stripline, or waveguide 
types of up-converter structures. In this model, the 
diode will be resonated in series, although the design 
data to follow can be applied also for designs where the 
diode is resonated in shunt. 

Fig. 4 shows the model up-converter circuit. This 
circuit uses a three-resonator wide-band input filter, a 
three-resonator moderately wide-band pump circuit 
filter, and a narrow-band, single-resonator output filter. 
The diode plus the length of high-impedance transmis- 
sion line introduces multiple resonances, one of which 
provides the resonance for the No. 1 resonator of the 
input filter, and another of which provides the reso- 
nance for the No. 1 resonator of the pump filter. Note 
that the diode representation includes Co’, Ci*, and R, 
in series along with La, the diode internal inductance. 
Since the number of resonators for the signal-input 
and pump filters may vary under different circum- 
stances, it is convenient for both the input and pump 
filters to define the resonator formed by the diode cir- 
cuit as the No. 1 resonator, even though the generator 
is actually at the other end. Under these conditions, Ry 
becomes the generator internal impedance instead of 
R,, but the power transmission properties are unaf- 
fected since the filter circuit obeys the reciprocity 
theorem. 

The performance of the circuit in Fig. 4 can best be 
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understood by considering its operation at the various 
frequencies of interest. Fig. 5 shows the equivalent cir- 
cuit of the input filter circuit at frequency f and the 
output circuit at the sideband frequency f’, with a box 
labeled X1.X.2 to represent the coupling effect of Cy be- 
tween energy components at frequencies f and f’. The 
operation of this box is different for lower-sideband up- 
converters than for upper-sideband up-converters,! and 
its operation will be summarized in Sections V and VII. 
Since the pump filter center frequency fo? is much 
higher than the center frequency fy of the input filter, 
the susceptance B,? of the No. 2 resonator of the pump 
filter in Fig. 4 will be very large, so that the high- 
impedance line is effectively short-circuited. At fre- 
quency fo, the high-impedance line is less than \/8 long 
and acts much like a lumped inductance which brings 
the diode to series resonance at that frequency. Since 
Co’, La, Rs, the impedance of the short-circuited line, and 
the impedance reflected through C\* (Z_ in Fig. 5) are 
effectively in series, the diode and high Zo line together 
may be represented as shown in Fig. 5 by R,, Z2, and the 
series resonator X}. 

At the pump-channel band-center frequency fo”, the 
impedance Z, in Fig. 4 will be reactive and nearly zero, 
and the diode plus the high Zp line will again exhibit 
resonance. The impedance seen looking into the high 
Z, line from the right in Fig. 4 is given by 


(Re + jXa) + 7Zo tan 6 
= ) 
j(Rs + jXa) tan 0 

+ 
Zo 


(16a) 


1 


where Rk,+7Xa is the total impedance developed across 
the diode along with the input filter reactance Z, in 
Fig. 4 at the pump frequency f?. Choosing Zo to be 
large as compared to R, (a condition which is easy to 
fulfill since R, might typically be around 4 ohms), to a 
good approximation (16a) becomes 


R, 


(1 - Fein) 
1 — — tan 6 
Zo 


(Zo tan 6 + Xa) 


Xa 
(1 — — tan ) 
Zo 


In most cases, the pump band center fo? will be above 
the self-resonant frequency of the diode; therefore, at 


Z=ReP+jxyv = 


(16b) 


foe the reactance X4q will usually be inductive, and 6 will 


be somewhat less than az. Under these conditions, the 
quantity in the denominators of (16b) will be some- 
what greater than one (since tan @ is negative), and will 
be slowly varying since Xa/Z» will be increasing with 
frequency, while tan @ will be decreasing in magnitude 
as the frequency increases. 
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Since 6 will be roughtly of the order of a, in some 
cases the denominators of the terms in (16b) will be very 
nearly one. Assuming for the moment that this is the 
case and that Z, for the input filter in Fig. 4 is negligible 
at fo”, then (16b) is equivalent to the impedance of 
R,, Co’, La, and the short-circuited, high Z» line all con- 
nected in series as shown in Fig. 6(a). If the denomina- 
tor in (16b) cannot be neglected, the magnitude of 
R,e+jXv will be affected as indicated by the equation, 
but the general nature of the performance will be 
qualitatively the same. 

It is thus seen that for frequencies in the vicinity of 
fo and fo?, Co*, La and the high Zp line operate at least 
qualitatively in a manner similar to the series-resonant 
circuit shown in Fig. 6(a), which has the reactance 
characteristics shown in Fig. 6(b). The reactance slope 
in the vicinity of fy will be important in the design of 
the signal-input filter, while the reactance slope in the 
vicinity of fo? will be important in the design of the 
pump input filter. At the sideband frequency fy’, the 
diode and high Zp» line in Fig. 4 will see reactive im- 
pedances Z, and Z,? of small magnitude, but the diode 
and line will not be resonant. However, the presence of 
the lightly-coupled sideband resonator in Fig. 4 will 
introduce a narrow resonance at either the lower- or 
upper-sideband frequencies indicated in Fig. 6(b). It is 
desirable that this sideband-output resonator be coupled 
as lightly as possible, consistent with low transmission 
loss, so that it will cause a minimum increase in the 
reactance slope in the vicinity of fo?. Since the up- 
converter would probably be followed by a fixed-tuned, 
superheterodyne receiver, an additional high-Q filter 
following the sideband-output filter might be desirable 
to ensure high attenuation at the image frequency of 
the superheterodyne receiver. 

Fig. 7 shows portions of a possible strip transmission 
line amplifier of the type under consideration. The diode 
used is in a computer diode package having a glass en- 
velope and wire leads. (The Hughes 1N896 diodes are 
examples of varactor diodes packaged in this manner.) 
The wire diode leads then serve as the high Zp line for 
resonating the diode. Resonators Nos. 2 and 3 for the 
signal-input filter are realized as two semi-lumped shunt 
resonators separated by a quarter-wavelength (at fre- 
quency fo) of line. The quarter-wavelength line elimi- 
nates the need to construct resonator No. 3 as a series 
resonator. The signal-input filter would be of wide band- 
width and it could be designed using the theory in 
Sections V or VII along with techniques for design of 
wide-band filters which have been previously treated.® 

The pump frequency must be able to shift an amount 
equal to the bandwidth of the signal-input filter, but 
since the pump frequency band is centered considerably 
higher than that of the signal-input band, the pump 


6G. L. Matthaei, “Design of wide-band (and narrow-band) band- 
pass microwave filters on the insertion loss basis,” IRE TRANS. ON 
MIcROWAVE THEORY AND TECHNIQUES, vol. MIT-8, pp. 580-593; 
November, 1960. 
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Fig. 6—(a) Approximate equivalent circuit for the diode resonator. 
(b) Reactance properties of the resonator in Fig. 4. 
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Fig. 7—A possible strip-transmission-line embodiment of the circuit 
in Fig. 4. The input, pump, and sideband filters are realized in 
practical microwave filter structures which will appear from the 
diode resonator circuit to be equivalent to those in Fig. 4 for the 
frequencies of interest. (a) Cutaway side view. (b) Top view with 
cover plate removed. 


filter fractional bandwidth will be relatively small. 
Pump filter resonators Nos. 2 and 3 are realized as 
quarter-wavelength direct-coupled, two-port resonators, 
since this construction is convenient for narrow or mod- 
erate bandwidths, and will not have any spurious pass 
bands until about three times the frequency of the first 
pass band.’ Each resonator bar is series-capacitance- 


7G. L. Matthaei, “Direct-coupled, band-pass filters with do/4 
Sy 1958 IRE NationaL ConveNTION REcorD, pt. 1, pp. 
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coupled at one end and shunt-inductance-coupled at 
the other. Looking to the right from the high Zp» wire in 
the diode circuit, the impedance of the structure ap- 
pears much like Z,? seen looking to the right from the 
/ high Z) line in Fig. 4 for frequencies in the vicinity 
of fo”. 

The sideband-output resonator shown in Fig. 7 is one- 
quarter-wavelength long; it is inductively coupled to the 
diode and high Zp wire circuit, but series-capacitance- 
coupled to the sideband-output line. 

The circuit shown in Fig. 7 appears to be a possible 
realization of an electronically-tunable up-converter 
having input frequencies centered around 1000 Me. 
However, other configurations are also possible. 


V. DETERMINATION OF SIGNAL-INPUT AND SIDEBAND- 
OuTPuT Circuit PARAMETERS FOR UppER-SIDEBAND 
Up-CONVERTERS 


When power is taken out at the upper sideband while 
the other components of the mixing process are sup- 
pressed, then the X12.X2: coupling box in Fig. 5 operates 
like an impedance inverter so that! 


X 2X01 
22 = ig 
Z! (17) 
and 
X2X01 
Zi = ) 18 
Z, (18) 


where, as indicated in Fig. 5, the unprimed impedances 
are evaluated at the signal-input frequency f, while the 
primed impedances are evaluated at the sideband fre- 
quency f’. By the Manley-Rowe equations,’ the net 
power P» entering the left side of X1.X2 is related to the 
net power P,’ entering the right side of X12.Xm by 


(19) 


where the minus sign indicates that if power flows into 
the left side, power will flow out of the right side. 

Let us define (Rs)nom as the nominal value of the 
signal-input filter impedance Z in Fig. 5. The im- 
pedance (Rs)nom is purely real and is equal to the re- 
sistive terminating impedance which would give best 
transmission through the fiter. In terms of Fig.3, 


r 
(Rb)nom = Ry = See (20) 
E 


r 
where 7, and 7; are from the lumped-element prototype. 
For purposes of computing the nominal gain of the 
amplifier, the signal-input filter will be replaced by a 
Thevenin equivalent generator having an internal 
impedance equal to (Ry)nom- At the sideband-output fre- 
quency fo’, the impedance Z;’ in Fig. 5 will be purely real 
and will be defined as 

Rees ee 


(21) 
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Using (17)—(19) and the above definitions for (Rs) nom 
and Ryo’, it can easily be shown that the nominal power 
gain for input signals in the vicinity of fy is 


I 
Je out 
Teast 


- fo’ 4A( Ro) nom Roo’ (X 12X21) 0 (22) 
fo Ue om tg Ry] (Roo SS) (X12X21)o} : 

where R,=R,’ is the diode loss resistance, Pavai is the 
available power at the input frequency fo, and Pout’ is 
the nominal output power at the sideband frequency 
fo’. The gain will vary somewhat across the tuning 
range due to variations in Z;, X1.Xm, and f’/f, but (22) 
gives what will herein be referred to as the nominal gain 
of an upper-sideband up-converter. 

In the design of an amplifier, the input and pump 
channel filters serve primarily as wide-band impedance- 
matching networks. Thus, the first step in the design 
process is to design the diode resonator circuit and then 
determine the slope parameters x; and x,” at frequencies 
fo and fo?, respectively, by use of the derivative form 
in Fig. 3(b). By Fig. 3(b), for a desired fractional band- 
width w of the input circuit, Ryo must be 


a= a(F) 
= 20) ——= 
b0 1 Ol, ) 


and then (R%5)nom is given by (20). The slope parameters 
x; and b; for the other resonators of the input filter may 
then be computed from the lumped-element prototype 
filter parameters, Ryo, and w, by use of Fig. 3(b) and (e). 

Knowing (Rs)nom and fo’/fo, and having values for 
(Xu)o, R,=R.’, and a=Ci/Cy for the diode, the maxi- 
mum gain will be obtained if Ryo’ matches Z;’ in Fig. 5, 
which requires that 


(23) 


' ; (X12Xa1)o 
Reo he eo Geen (24) 
be satisfied. This can be expressed as 
1 
Reo = Re E ae serch Os, 
where 
Ry nom Re) nom d 
and 
v= = = . (27) 
Inserting (25) in (22) yields, after some manipulation, 
Peet 
Payail . ip (28) 
(io AGREE Oise er il 
fo 
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If this analysis is carried out on the dual basis for the 
case of the diode resonated in shunt, we find the dual 


expressions 


(Ag ch ; 
Gio = bw ( — }, (29) 
oes Ge. 
‘gy 
(Ga) nem == G, a 2 Goo, (30) 
ASE; 
G E Peeeees | (31) 
Ca ! 
60 d mere 1) 
G nom G nom 
Fea p)nom _ (Go) Qa (32) 
Ga (Bii)o 


and (27) and (28) apply as before. The parameter Ga’ in 
(31) is again the equivalent shunt diode loss conduct- 
ance given by (9). In order to facilitate design calcu- 
lations, a plot of (Pout /Pavait)/(fo /fo) vs ¥ is shown in 
Fig. 8 for various values of 7. The reference form for the 
band-pass filters is the dual of that in Fig. 3; the slope 
parameters and terminations for this dual filter may be 
conveniently specified directly in terms of the low-pass 
prototype in Fig. 1. 


VI. Pump Circuit IMPEDANCE-MATCHING FILTER 


In order that the pump channel have a nearly uni- 
form reflection loss across the pump frequency band, it 
is necessary to design an impedance-matching filter as 
shown in Fig. 9. In this case, the resistance RK,” is given 
by (16b). The No. 1 resonator, having a reactance Xi? 
[given by (16b) for the case of Figs. 4, 7 and 9] is 
characterized by the slope parameter x.” computed at 
frequency fo? in Fig. 6(b) using the derivative form 
in Fig. 3(b). (Recall that in this case, the resonance 
of the No. 1 resonator is created by the second series 
resonance of the diode circuit, as discussed in Section 
IV.) The fractional bandwidth w, of the pump filter is 


fo 


Uy 3 Ys) 


fo? 


where w is again the desired fractional bandwidth of the 
input circuit. By Fig. 3(b), 


(=) 
1L/ » 
where the subscript p on (Qil/rz)» is introduced to 
indicate that these parameters refer to the low-pass 


prototype for design of the pump filter. Since, in ac- 
cord with Figs. 1 and 3, 


(33) 


X1?Wp 


Roo? 


; (34) 


Roo” 
Vege 


(35) 
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Fig. 8—Chart for determining the gain of upper-sideband up-con- 
verters. [T= (Rs) nom/R, for the case of a diode resonated in series, 
and T= (Gs)nom/Ga for the case of a diode resonated in shunt. 
This chart does not apply accurately for the shunt resonance case 
if the diode has significant series inductive reactance. | 
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AT PUMP FREQUENCY f° 


Fig. 9—A circuit approximately equivalent to that in 
Fig. 4 for energy at the pump frequency f?. 


we niay write (34) as 


a) 

Tg /p 
where, for the case of Figs. 4, 7 and 9, R,? is given by 
(16b). The problem now focuses on finding a low-pass 
prototype having (Qy//rg)p as given by (36), with a de- 
sired amplitude of Tchebycheff pass-band attenuation 
ripple, and with a minimum amount of reflection loss in 
the pass band. We are stuck with the prototype im- 
pedance jQh+7,, and the other elements of the proto- 
type filter are introduced to optimize the impedance 
across the band Q=0 to Q=Q,. 

Certain aspects of impedance-matching problems of 
this type were treated by Bode.® Fano treated the gen- 
eral limitations on lossless impedance matching, and 
also the design of certain Tchebycheff impedance- 


XW 


Rp (36) 


wes H. W. Bode, “Network Analysis and Feedback Amplifier De- 
sgn D. Van Nostrand Co., Inc., New York, N. Y., pp. 363-368; 
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matching networks,® and Carlin,!? Green," Barton,” 
Matthaei,® and others have made various contribu- 
tions. 

The problem at hand can be treated as follows. Sup- 
pose that it is desired to use m resonators in the pump- 
circuit filter and that a Tchebycheff ripple of 7, db is 
specified. Compute 


If 
A = antilogi) — » Sih 
£10 10 (37) 
Se ar 
sinh— ——— 
Ney 
d = sinh — (38) 
n 
and 
19 eG 
e=d-2(74)sin= (39) 
Qu, 2n 


Then, for given 7,/(Qi,), T,, and n, the reflection co- 
efficient value at the maximum reflection points of the 
Tchebycheff pass band is given by 


U,(e) 


U,(d) - 


| p Ee =: 
where U,(e) and U,(d) are polynomials tabulated by 
Matthaei and plotted in Fig. 10 of this paper. Under 
these conditions, the minimum pump power delivered 
to the diode in the pump filter pass band is 


Pransmitted / 
ee [oa 
CP axe)” 

times the available power (Payaii)” of the pump signal 
generator. It will be found that the larger the slope 
parameter x,” and the larger the fractional bandwidth 
W», the smaller A must be and the more available power 
will be required. 

Having values for 7,/(Qih), T;, and |p ea tOnetie 
pump filter, the actual element values of the low-pass 
prototype are in this case obtained most easily by use 
of equations given by Green."* It should be under- 
stood that (37)—(41) apply only to simple ladder cir- 


(41) 


9R. M. Fano, “Theoretical limitations on the broadband match- 
ing of arbitrary impedances,” J. Franklin Inst., vol. 249, pp. 57-83, 
139-154; January and February, 1950. ; ; 

10 H. J. Carlin, “Synthesis techniques for gain-bandwidth opti- 
mization in passive transducers,” Proc. IRE, vol. 48, pp. 1705-1714; 
October, 1960. an 4 

1 E, Green, “Amplitude-Frequency Characteristics of Ladder 
Networks,” Marconi’s Wireless Telegraph Co., Ltd., Chelmsford, 
Essex, England, pp. 62-78; 1954. : 

2B. F, Barton, “Design of Efficient Coupling Networks,” Elec- 
tronic Defense Group, University of Michigan, Ann Arbor, Tech. 
Rept. No. 44, Contract DA 36-039 SC-63203; March, 1955. 

18 G. L. Matthaei, “Synthesis of Tchebycheff impedance-matching 
~ networks, filters, and interstages,” IRE TRANS. ON CIRCUIT THEORY, 
vol. CT-3, pp. 162-172; September, 1956. ; 

4 E, Green, “Synthesis of ladder networks to give Butterworth 
or Chebyshev response in the pass band,” Proc. IEE, pt. 4, Mono- 
graph no. 88; 1954. 
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Fig. 10—Chart for use in determining the reflection loss required in 
the pump circuit when using an optimum Tchebycheff impedance- 
matching filter. 


cuits such as that in Fig. 1 or their band-pass equiva- 
lents (such as that in Fig. 9). 


VII. DETERMINATION OF SIGNAL-INPUT AND SIDEBAND- 
OuTPuT Circuit PARAMETERS FOR LOWER-SIDEBAND 
Up-CONVERTERS 


If power is taken out at the lower sideband instead 
of at the upper sideband, (17)—(19) become, respec- 
tively, 


X42X 91 
23. ) (42) 
(Za')* 
y X 2X4 
Le aioe ) (43) 
(Z.)* 
and 
P / iy 
ee (44) 
P» i 


where the asterisk indicates the complex conjugate. 
Because of the negative sign in (42), this type of opera- 
tion gives gain at the input frequency f, due to negative 
input resistance. As is implied by (44), the power P,’ 
(at the lower-sideband frequency f’) going out the 
right side of Xy2X2 in Fig. 5 is larger than the power P2 
at frequency f going out the left side by the factor 
f'/f. Thus, whatever gain is achieved at the input fre- 
quency is increased at the XX box by the factor 
f'/f for power taken out at the sideband frequency f’. 

In this case, (20) and (21) apply as before, but (22) 
becomes 


/ 
P out 
avail 


— (=) 4( Rp) nomlpo’ (X 12X 21) 0 
fo [(Ro) nom a R,| (Roo! me R,’) al (X12Xo1)o} ; 


» (45) 
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which gives the nominal gain of the amplifier. The value 
of Ryo for the input filter is determined by (23) as be- 
fore, but Ryo’ cannot be determined by impedance 
matching to maximize the gain as was done to obtain 
(24), since the maximum possible gain is now infinite 
(yielding oscillation). In this case, we define a stability 
parameter 


ss (Ris) mens =F ce x (Re) nem = ae ; (46) 


| Zs lian [ee 
Wey ae aes 


which is seen to be the ratio of the nominal positive re- 
sistance of Z, in Fig. 5 to the magnitude of the negative 
resistance Z) seen looking in the opposite direction. 
Parameter D fixes the nominal value of the negative 
resistance gain of the amplifier, and for stability, D must 
be greater than one. By (46) 

D(X 12X01) 0 


Fe eee hE 7 (47) 
eo OLR 


ID, 
-aigey 
v(7T + 1) 


where 7’ and v have the same meaning as in (26) and 
(27), respectively. Substituting (48) in (45) gives 


(48) 


jee 4 D ( 21 D 
Der | mat eel 
bs ero) 
fe (D — 1) 
fo 


The ratio (Pout’/Povait)/(fo /fo) is plotted vs »v in Fig. 
11(a)-(c) for D=2, 3, 4 and various T values. These 
charts may also be used for the dual case where 
the diode is resonated in shunt by defining 7 as in (32). 
For the dual shunt-diode case, (48) becomes the dual 


equation 
D 
G'| ——— - 1]. 
vX(T + 1) 


VIII. EstIMATED PERFORMANCE OF SOME 
DEsIGN EXAMPLES 


Gyo = (50) 


Let us estimate the performance of some up-con- 
verters having an input band center of fy =1 kMc. The 
construction in Figs. 4 and 7 will be assumed using a 
Hughes 1N896 diode which has the computer-type 
package shown in Fig. 7. The lead wires on the diode 
will serve as the high Zo line. The wires are 0.020 inch in 
diameter, and if a 0.500-inch ground-plane spacing is 
used, they will have an impedance of 207 ohms. The 
diode will be assumed to have a Coy value of 1.2 upyf 
(corresponding to a zero bias capacitance of about 2.3 
uf). The series inductance La of the diode itself will be 
taken as 4X10-° h, and the operating Q will be taken as 
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Fig. 11—Charts for determining the gain of lower-sideband up-con- 
verters. [T= (Rs) nom/Rs for the case of a diode resonated in series, 
and T= (Gs)rom/Ga for the case of a diode resonated in shunt. 
Parameter D is the ratio of the magnitude of positive resistance 
to negative resistance seen at the diode (or for the shunt case, 
positive conductance to negative conductance). These charts 
do not apply accurately for the shunt resonance case if the diode 
has significant series inductive reactance.] 


Qa=35 at fo=1kMc (from the manufacturer’s data, this 
appears to correspond to a cutoff frequency of around 
70 kMc for this type of abrupt junction diode). The 
a=C,/Cy parameter of the diode will be taken to be 
0.25: 

In terms of the approximate equivalent circuit in Fig. 
6(a), for the diode resonator, 


2mfCo® 


A) = 


+ 2nfLa + Zo tan @} (51) 


where f, is the frequency shown in Fig. 6(b) for which 
the high Zo line is. \/4 long. Setting f=f)=109, 
Cot =Co(1 —a?) =1.17X10-“fd, La=4xX10-%h, and 
Zo = 207, we find that f,=3.21 kMc is required, which 
in turn calls for the Z) = 207-ohm line to be 7=0.922 inch — 


pihesthed BW. 
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long. The slope parameter x; at frequency fy is obtained 
from 


W—2rfy=—W, 


+ 2nfoLla 


4 (=) Zo . 
2h, Tho 


cos” (=) 
Die 


which gives x,=145 ohms. Solution of (51) shows that 
fo? in Fig. 6(b) will occur at about fo? =5.45 kMc, and by 
calculations similar to (52), the slope parameter there 
is about x? =430 ohms. 

Let us first consider an upper-sideband up-converter 
having an input tuning band-edge ratio of one-half 
octave (.e., for the input filter, fo/fi=+/2 where f; and 
fe are defined in Fig. 2). This calls for w=0.348. The in- 
put filter will be assumed to have =3 resonators and 
1-db Tchebycheff ripple. From tables 4° of low-pass 
prototype filter elements, we find that (Qyl;/r,) =2.024 
and r,=r, for this case. Then by (20) and (23), we ob- 
tain (Ry)nom = Roo =25 Ohms. Since R, = (X11)0/Q) = 3.86 
ohms, T= (R3)nom/R;=6.5. In this upper-sideband case, 
Fo =fo? +fo=5.45+1=6.45 kMc which gives, by (27), 
v=0.291. By use of 7 and v in Fig. 8, along with the 
fact that fo'/fo=6.45, we obtain P’out/Pavaii=3.48 or a 
nominal gain of 5.3 db. By (25), it is found that the 
sideband resonator should be adjusted to couple a re- 
sistance of Ryo’ =9.95 ohms into the diode circuit at 
the upper-sideband frequency fo’ = 6.45 kMc. 

Using x”=430 ohms, w,=wfo/fo? =0.0638, and 
R#=R,=3.86 ohms in (36), we get (Qh/r,),=7.1 for 
the pump filter low-pass prototype. A Tchebycheff rip- 
ple of 0.5 db using an x=3 resonator pump filter will be 
specified which by (38) calls for d=0.6265. Then by 
(39), e=0.4855. By Fig. 10, U;(0.6265)=3.05 and 
U;(0.4855) =2.18. By (40), |p|max=0.715. By (41), a 
minimum power of A (Payait)?=0.49 (Pavaii)? will be de- 
livered to the diode where (Puyaii)? is the available 


l 
2 2mfoCo8 


(52) 


4 In computing R,? and x”, the denominator in (16b) was neg- 
lected and the representation in Fig. 6(a) was used. This should 
cause little error in the performance estimates about to be computed 
since, as has been mentioned in Section IV, the denominator of (16b) 
will be varying quite slowly. It can therefore be regarded as constant 
within the frequency range of interest. As a constant, it would de- 
crease the size of both R,? and x” proportionally (about 25 per cent 
in this case), but would cancel out in the computation given by (36). 
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pump power. Assuming that the diode is to be pumped 
with a peak-to-peak voltage of 6 v, it is estimated that 
the diode will absorb about 30 mw of pump power when 
f? =fo? =5.45 kMc. Thus, the required available pump 
power will be about (Pavai)? = 30/A =30/0.49 =61 mw, 
neglecting dissipation loss in the pump input filter. A 
1-db dissipation loss would raise this figure to about 
77 mw required available power. 

Let us now make analogous approximate calculations 
as described in Section VII using the same diode and 
diode circuit having fo=1kMc, fo? =5.45 kMc, «1=145 
ohms, x.” =430 ohms, etc., as before. Again, w=0.348 
will be used (half-octave input tuning range). In this 
case, the output will be taken at the lower-sideband 
frequency: fo’ =fo? —fo =5.45 —1=4.45 kMc. Since this 
gives negative resistance operation which increases the 
sensitivity to input impedance variations, the n=3 
resonator input filter will be assumed to have its 
Tchebycheff ripple reduced to 0.25 db which calls for 
(Quh,/7,) = 1.3034. 

The calculations as outlined in Section VII call for 
(Rs) nom = Roo = 38.7 ohms, Ryo’ = 8.0 ohms (when D=2). 
Assuming a stability parameter of D=2, the nominal 
power gain would be 13.4 db, 8.1 db greater than for 
corresponding upper-sideband design. The pump power 
required should be the same in both cases. 

Similar approximate calculations for a one-octave- 
tuning-range, lower-sideband up-converter were made, 
again assuming the same diode circuit and the same 
So, fo?, fo’ =4.45 kMc values. A value of 0.25-db Tcheby- 
cheff ripple was assumed for the 2=3 resonator input 
filter. Since stability might be more of a problem using 
an octave input bandwidth, D=3 was used. In this 
case, (Rs)non= 18.8, Roo =5.45, the nominal gain’ is 
8.8 db, A =0.273, and the required available pump power 
is about (Payai)® =30/0.273 =110 mw, neglecting pump 
filter dissipation loss and assuming that the diode needs 
30 mw of absorbed pump power. If the pump filter has 
1-db dissipation loss, the required available power 
would be about 138 mw. 


IX. CONCLUSIONS 


Properly designed devices of this type using voltage- 
tunable pump oscillators should have tuning ranges of 
the order of a half-octave to an octave, fast tuning 
capability, a useful amount of gain, no image response, 
and a low noise figure. Though some reflection loss in 
the pump circuit will be unavoidable, it appears that 
in many cases existing voltage-tunable oscillators will 
have adequate incident power. 


436 


IRE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 


September 


Anisotropic Properties of Strip-Type Artificial Dielectric* 


N. J. KOLETTIS} anp R. E. COLLIN{, SENIOR MEMBER, IRE 


Summary—Theoretical formulas for the propagation phase con- 
stant of a two-dimensional strip medium are presented for general 
directions of propagation. In addition a number of experimental re- 
sults are included that verify the validity of the theory. Some of the 
difficulties encountered in defining equivalent dielectric constants 
for this medium are also pointed out. 


duced by Kock,! have received extensive inves- 

tigation during the past decade. However, for 
the most part the detailed study of the anisotropic 
properties of typical structures has not been carried out. 
Most studies have been restricted to angles of propaga- 
tion through the medium and to polarizations, so that 
the artificial dielectric behaved essentially as an iso- 
tropic medium. Since the anisotropic properties can be 
useful in practice (for example, in wave polarizers), it is 
worthwhile to analyze these properties for typical 
structures. From another point of view, anisotropic 
properties are usually present and hence a knowledge of 
these properties is needed in order to determine the 
actual behavior of an artificial dielectric medium in a 
given application. 

The two-dimensional strip-type medium illustrated 
in Fig. 1 is a typical example of an artificial dielectric 
that exhibits pronounced anisotropic properties. For 
propagation through the medium in a direction normal 
to the broad side of the strips (along « axis) and with the 
electric field directed along the zg axis, the medium be- 
haves as an isotropic phase-delay dielectric. For the op- 
posite polarization the medium behaves as a_phase- 
advance dielectric. The properties of this medium with 
propagation restricted to be along the x-axis, has been 
analyzed by Brown,? and Cohn,? as well as others. 
If a wave with the electric field polarized in the xz plane 
is permitted to propagate through the medium at an 
effective angle 6,’ with respect to the x axis, then the 
medium exhibits electric anisotropic properties. The 
analysis for this more general case has been presented 
by Collin.* The purpose of the present paper is to give 


\ RTIFICIAL dielectric media, originally intro- 


* Received by the PGMTT, April 13, 1961; revised manuscript 
received, June 30, 1961. This work was sponsored by the Electronics 
Res. Directorate of the AF Cambridge Res. Lab. under Contract No. 
AF 19(604)-3887. 

ft RCA Lab., Princeton, N. J. Formerly with the Case Institute 
of Technology, Cleveland, Ohio. 
nat Elec. Engrg. Dept., Case Institute of Technology, Cleveland, 

io. 

1W. E. Kock, “Metallic delay lenses,” Bell Syst. Tech. J., vol. 27, 
pp. 58-82; January, 1958. 

2 J. Brown, “The design of metallic delay dielectrics,” Proc. IEE, 
vol. 97, pt. 3, pp. 45-48; January, 1960. 

3S. B. Cohn, “Analysis of the metal strip delay structure for 
microwave lenses,” J. Appl. Phys., vol. 20, pp. 257-262; March, 1949. 

4R. E. Collin, “Field Theory of Guided Waves,” McGraw-Hill 
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a summary of the theroetical results obtained to date, 
together with experimental results that show that the 
theoretical formulas specifying the anisotropic behavior 
of the medium give results of acceptable accuracy for 
design purposes. 


Fig. 1—Cross-section view of two-dimensional strip-type artificial 
dielectric, strips infinite in length along y-axis. 


MopEs OF PROPAGATION 


Before presenting the theoretical formulas that per- 
mit the equivalent dielectric constant of the medium to 
be evaluated, a qualitative discussion of the possible 
modes of propagation in the medium will be given. If 
propagation is restricted to be in the x direction with the 
E vector along the zg axis then the dominant mode is a 
TEM mode and all higher order modes (assumed 
evanescent) are TM modes with respect to the x axis. 
These may be labeled as (TM), modes. If propagation is 
permitted to be at an angle to the x axis but in the xz 
plane and furthermore with £ in the xz plane then all 
the modes are (TM), modes. For the opposite polariza- 
tion with the # vector in the xz plane all the modes are 
(TE), modes. For more general directions of propaga- 
tion through the medium, so that the propagation vec- 
tor has a component along all three axes, the possible 
modes are not of the above type. One set of possible 
modes of propagation is still characterized by the ab- 
sence of a y component of electric field, but it does have 
an x component of 7, and hence these may be labeled 
as (TE), modes, 7.e., transverse electric modes with 
respect to the y axis. These are the same types of modes 
that occur in connection with capacitive diaphragms in 
rectangular waveguides. These modes may be ex- 
pressed in terms of the field component H,. Between 


_ °H.M. Altschuler and L. O. Goldstone, “On network representa- 
tions of certain obstacles in waveguide regions,” IRE TRANS. ON 
MICROWAVE THEORY AND TECHNIQUES, vol. MTT-7, pp. 213-221; 
April, 1959. 
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any two strips the dominant mode in this class is of the 
form, with the z axis considered as the axis of propaga- 
tion, 


Hy = e—tho(nyytnez) ; 


where m,, m, are components of the unit wave normal n. 
In addition there are field components H, and E,. This 
mode has no low-frequency cutoff. 

The other class of possible modes is characterized by 
the absence of a y component of magnetic field and may 
be classed as (TM), modes. These modes may be ex- 
pressed in terms of E,. The dominant mode in this class 
that may exist between two strips, again considering the 
Z axis as the axis of propagation, is of the form 


TX 
Ei, = sin — e 2k0(nyytnzz) 
Ss 


where (Rom)? = ko?(1 —n,?) — (/s)?. For the usual values 
of s occurring in practice, (s <A 0/2), m, is imaginary and 
this mode is below cutoff. Thus for the (TM), modes the 
strip-medium has a low frequency stop band. Only for 
S>o/2 can a propagating (TM), mode exist, and when 
it does the medium functions as a phase-advance me- 
dium for these modes. 

In this paper we are concerned only with the case of 
the (TE), modes. For these, the component of the 
electric field along the z axis gives rise to a large z- 
directed electric dipole moment in each strip and hence 
an effective dielectric constant x, in the z direction. The 
strips are infinitely thin in the x direction; hence F, 
produces no dipole polarization in the x direction, so 
the effective dielectric constant in this direction is 
unity. The dielectric constant in the y direction may be 
arbitrarily chosen since the (TE), modes do not have 
an electric field component in this direction. For rea- 
sons to be explained later, it is convenient to chose the 
effective dielectric constant in the y direction as x. 

The specification of equivalent dielectric constants 
for an artificial dielectric medium is just a convenient 
artifice that will permit an analogy to be made between 
the artificial medium and a homogeneous anisotropic di- 
electric medium. Such an analogy is useful only if the 
propagation of a wave through the homogeneous di- 
electric follows the same law as that for the artificial 
dielectric. When this is true the refractive properties of 
the artificial medium are the same as for the homo- 
geneous dielectric and the design of a microwave lens 
may then be based on the properties of the equivalent 
homogeneous dielectric. It should be noted at this time 
that a homogeneous medium which is fully equivalent 
to the artificial dielectric medium for all possible modes 
of propagation can usually not be specified. Thus, what- 
ever analogy can be established is usually valid only 
under certain restrictive conditions on wave polariza- 
tion and direction of propagation. 

The homogeneous anisotropic dielectric medium 
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chosen as an analogy for the strip-medium, for a (TE), 
mode of propagation in the strip medium, is a uniaxial 
medium with dielectric constants 1, k., Kk. along the x, y 
and z axis respectively. The basis for this analogy is 
that for the strip medium, the low frequency solution 
for the propagation phase constant 6, along the wave 
normal n, satisfies the same eigenvalue equation as for 
the homogeneous medium. Furthermore, for propaga- 
tion in the xz plane, the extraordinary wave in the 
uniaxial medium has the same field components as the 
(TE), mode in the strip medium. The analogy is, how- 
ever, not a complete one, since for propagation through 
the medium with n having a y component, the (TE), 
mode has a zero y component of electric field while the 
extraordinary wave in the homogeneous uniaxial di- 
electric has a zero component of magnetic field along 
the optic axis (« axis in the present case).® Neverthe- 
less the eigenvalue equation and refractive properties 
of the two media are the same for the respective modes 
of propagation that they support. 


THEORETICAL RESULTS FOR EQUIVALENT 
DIELECTRIC CONSTANT 


The strip medium illustrated in Fig. 1 is a double 
periodic medium and may be analyzed in terms of 
modes propagating along the z axis or in terms of modes 
propagating along the x axis. In the first method, which 
will be referred to as the z-axis solution, each layer of 
strips presents a discontinuity at each interface. The 
reflection and transmission coefficients that charac- 
terize each interface are those associated with the inter- 
face between a set of semi-infinite parallel plates and 
free space. The latter problem has been solved by Carl- 
son and Heins.’ Using their results, it is found that the 
equivalent circuit of the strip medium for modes propa- 
gating in the z direction is that illustrated in Fig. 2.4 
The circuit in Fig. 2 is arrived at by assuming that the 
strip spacing s is small compared with a wavelength and 
that neither a nor } are so small that higher-order mode 
interaction between adjacent interfaces is important. 
For a TEM mode incident, in the xz plane, on the strip 
medium at an angle 6; relative to the z axis and with 
the electric vector in the xz plane, the propagation 
phase constant Bn, through the medium is given by 


(1 — p?) cos Bn,d 


5 
= cos Ex + bcos 6,;) — ko— (1 — cos ;) In a] 
Tv 
— p’ cos Ex — bcos 6) 


Fe ahs eer abs a], (1) 
au 


6 Collin, op. cit., Sec. 3.7. 

7 J. F. Carlson and A. E. Heins, “The reflection of an electromag- 
netic plane wave by an infinite set of plates,” Quart. Appl. Math., 
vol. 4, pp. 313-329; January, 1957, and vol. 5, pp. 82-88; April, 1947. 
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where 
p = (1 — cos 6;)/(1 + cos 4), 


and the other parameters are give in Fig. 1. This equa- 
tion results from an analysis of the periodic structure 
illustrated in Fig. 2. When the frequency is low so that 
kod approaches zero, the cosine terms in (1) may be 
approximated according to the relation cos x =1—x?/2. 
After simplification (1) reduces to 


b 
(81,)7 = Re E — (= + = In :) sin? a]. Q) 


T 


ko = 21/Xo 


b+2|In4 a-2|n4 


Ze = 2¢= Ko 
K, COS 8; 


Grae constant fe iaie 
= k, cos constant = Ky 


Fig. 2—Equivalent circuit for strip medium for 
modes propagating in the z direction. 


In the present case the propagation phase constant 
along the x axis is ky sin 6;=@n, and is known. If 
propagation takes place along the y axis also with a 
phase constant Bn,, then in (1) and (2) ko? must be re- 
placed by k? =k? —6?n,7 and Bn, becomes equal to F sin 
6;. The basis for this substitution is the same as that 
used to obtain the solution for a capacitive diaphragm 
in a rectangular guide from the solution for a diaphragm 
inan infinitely wide parallel-plate transmission line. Thus 
for general directions of propagation the low frequency 
limit of (1) is 


b Ss 
(iy Sha Gi) = ane (= += in 1). 3) 


The static solution for the dielectric constant of the 
strip medium in the gz direction, as obtained by con- 
formal mapping, is given by® 
d 
ae eee a (4) 
S 
6+ —In4 
z 
when both s/a and s/b are small. Substituting this into 
(3) and solving for 8? gives 


Keko- 
Na? = (1 i Nz?) Ke 


This, however, is the eigenvalue equation for the ex- 
traordinary wave in a homogeneous uniaxial aniso- 


8 J. Howes and E. A. N. Whitehead, “The Refractive Index of a 
Dielectric Loaded with Thin Metal Strips,” Elliot Brothters Ltd., 
London, England, Elliot Bros. Rept. 119; July, 1949. 

N. Kolettis, “Conformal Mapping Solution for Equivalent Rela- 
tive Permittivity of an Artificial Strip Dielectric Medium,” Case 
Inst. Tech., Cleveland, Ohio, Scientific Rept. No. 2, issued on Con- 
tract AF 19(604)-3887; January 15, 1959. 
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tropic dielectric with dielectric constants 1, ke, Ke, along 
the x, y, z axis, respectively. It is for this reason that 
the homogeneous uniaxial medium is chosen as an anal- 
ogy for the strip medium, and also the reason why 
the strip medium is assumed to have a dielectric con- 
stant x, in the y direction. Another similar analogy is 
also possible, and this is discussed in the Appendix. 

For high frequencies 8u, must be found from (1). By 
analogy with (3) the equivalent dielectric constant x, 
may then be defined by 


(Bn)? = ky? — (Bny)? = ea.) 
Or 


(Bnz)? 
ae 6 
ko? — (Bnz)? = (Bn,)* e 


In general, x. is a function of frequency except in the 
low frequency range where (4) is valid. 

The solution provided by (1) may be expected to 
give good accuracy for small values of s. For larger 
values of s higher mode interaction becomes significant 
and should be accounted for. On the other hand, by 
solving the same problem in terms of modes propagat- 
ing in the x direction, the condition that higher order 
mode interaction in this approach be negligible, is that 
the spacing s be large. Thus the x-axis solution is com- 
plementary to the z-axis solution in that it will yield 
good results for large s while the latter yields good ac- 
curacy for small s. A combination of the two methods 
will then enable accurate estimates of x. to be obtained 
for all values of s likely to be encountered. 

For the x-axis solution each row of strips in an x =con- 
stant plane constitute a capacitive susceptance 7B. Thus 
the equivalent circuit for modes propagating in the x 
direction is a transmission line loaded at intervals s by 
shunt susceptances 7B. The propagation phase constant 
Bn, in the x direction is a solution of 


cos Bnzs = cos hs — ms sin hs, (7) 


where 
h? = [Ro? — (Bmy)?] cos? 0,’ 


and 6,’ is the angle of propagation, relative to the x-axis, 
for the TEM mode between strips when ,=0. Alterna- 
tively h may be expressed as ky cos 6, where cos 0; is the 
direction cosine between the x axis and the wave normal 
of the TEM wave existing between adjacent gratings, 
an expression valid for all m,. In the analysis of the prob- 
lem Bn, and Bn, are usually specified and hence 


h? = ky? — (Bny)? a (Bn.)?, (8) 


since for the TEM mode, between gratings, the magni- 
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tude of the propagation phase constant must equal ko. 
An approximate value of B is given by? 
a 2kod cos ai bene a d (1 — a)?[(1 — a?/4)(Ay + A_) + 402A, A_] 
9 2 2/ 9) 
tg 2d 2 (1 — a?/4) + a2 (1 + 02/2 — at/8)(Ay + A_) + ie ; 
where mounted on polyfoam, then x,p!/?f/, is the resonant fre- 
Aoeein 40/20 quency of the same mode for a cavity filled with strips 
é surrounded by free space. k, is the dielectric constant of 
an polyfoam, about 1.08 in the experiments carried out. 
ie ie kod cos 0,/\ 27-12 This frequency scaling was used to convert all results 
Az =|1+—sin@,’ — )| — 1. to that for strips surrounded by free space, so as to 
Tv 27 


In this equation, ky cos @;’=h and ky sin 6,’ = (ky? — hh)! 
and h is given by (8) in general. When 1,=0, 6,’ is the 
angle of incidence relative to the x axis in Fig. 1, but 
for n,~0, k,” is to be replaced by k?=ky?—(@n,)? and 

h=ky cos 6,’ becomes (ko? —6?n,?)! cos 6,’ while ko sin 
0,’ becomes (Ro? — B?,?)"/? sin 8,’ = (Ro? — Bn,2 — h?) 2 = Bn,. 
Once $n, has been found from (7) the equivalent dielec- 
tric constant kx, may be obtained from (6) since Bn, and 
Gn, are determined by the known incident wave on the 
medium. 


EXPERIMENTAL RESULTS 


Measurement of the propagation constant Bn. was 
carried out for six different samples three of which were 
eight periods long, while the other three were ten pe- 
riods long. The samples were constructed by mounting 
metallic foil strips between sheets of polyfoam as 
illustrated in Fig. 3. Each sample was made to fit inside 
a rectangular guide of internal dimensions 1 7/8” by 
7/8” and having a removable top plate with a centered 
slot. The foil strips were bent over along the sides of 
each sample in order to make a reasonably good contact 
with the side walls of the guide. At each end of the 
sample the guide was closed off by a short-circuiting 
block so as to form a resonant cavity. The cavity was 
excited through a small circular aperture in one end 
plate and with an Ey; mode incident in the main guide. 
When the cavity resonates Bu,=n7/Nd where n is the 
nth resonant mode and J is the total number of sec- 
tions in the sample. The available pass band of the main 
guide for the E11 mode permitted from 4 to 7 resonant 
frequencies to be found. Resonance was detected by 
means of a probe inserted into the cavity through the 
centered slot along the top surface of the cavity. 

The measured values of 6, for strips mounted on 
polyfoam may be converted to equivalent results for 
strips in free space by frequency scaling. If fi is a 
resonant frequency for the cavity filled with strips 


9N. Marcuvitz, “Waveguide Handbook,” McGraw-Hill Book 
Co., Inc., New York, N. Y., vol. 10 of M.I.T. Rad. Lab. Series, Sec. 
5.18, Eq. (1a); 1951. 


avoid the need to consider the properties of the support- 
ing medium in the comparison of the theoretical and 
experimental results. 

In the z-axis solution for the loaded guide Bu,=7/d; 
and Bn,=7/d2 where d; is the guide height (7/8”) and 
d, is the guide width (1 7/8”). In the x-axis solution 
Bn, and Bn, are also fixed by the guide dimensions and 
hence (7) and (9) must be solved simultaneously for the 
parameter h. Eq. (8) is then readily solved for Bn,. The 
solution for h was obtained by plotting B vs h as de- 
termined by (7) and (9). The point of intersection of the 
two curves determines the value of h# that satisfied both 
equations. 


polyfoam spacers 


Fig. 3—Typical strip medium sample with 3 layers and 8 sections. 


In Figs. 4 and 5 the theoretical and experimental 
values of Bn, are plotted as a function of frequency. The 
frequency has been scaled so that the results are for 
strips surrounded by free space. For Fig. 4 the parame- 
ters of the three samples measured are a=1.143 cm, 
b=0:381 cm, s=0.554 cm, 0-739 cm, 1:11 cm, respec: 
tively and N=8, 7.e., the samples are 8 sections long. 
The samples for Fig. 5 are 10 sections long and a=0.762 
cm,-6=0.254 cm, and s=0:554 cm, 0.739'cm, and" 1711 
cm, respectively. An examination of the experimental 
results shows very clearly that for small values of s the 
experimental values lie very close to the theoretical 
curve based on the z-axis solution, while for the largest 
value of s used the experimental results lie close to the 
theoretical curve based on the x-axis solution. Fig. 6 is 


a plot of the computed value of x. [based on (0) | asa 


function of frequency for the two samples for which 
s=0.739 cm. At the higher frequencies k, rises rapidly 
in value since the edge of the first pass band is being 
approached. 


440 IRE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES September 


+ 


Bnz rad/cm. 
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x-axis solution 
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Z-axis solution “t_z-axis solution 
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2.6 


x-axis solution 


Xaxis solution 


z-axis solution “2_ z-axis solution 
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Fig. 4—Computed and measured values of Bn, for a=1.143 i — S 
pee a el cans Ftexcerimientallipointa, (a) «= 0.554 Gui ec < ccUsee d Oh ets or ae eae 
0.739 cm, (c) s=1.11 cm. s=0.739 cm, (c) s= 111 cm. points. (a) s=0.554 cm, (b) 


Ke 
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Fig. 6—Computed value of x. for two sample i 
] s for which s=0. 
cm, based on experimental values of Bn; oon Figs. 4(b) se se 
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The theoretical solution for Bm., in the case n,=0, 
yields Bn, as an implicit function of the free space wave 
number ko. When 2,0 the solution for Bn, at any given 
frequency f; may be found by simply replacing 
ko =2fi(uoeo)? by k= (ko? —B2n,2)'2. Thus the curves 
presented in Figs. 4, 5, and 6 as a function of frequency 
with Bn,=7/d2, may be interpreted as results for the 
case 1, =0, 2.e., propagation in xz plane only, by replac- 
ing the frequency scale f by a new frequency scale f, 


where 
1 1/2 
Hs sf E Be | , 


CONCLUSIONS 


Experimental results have been presented which show 
that for small strip spacing s, the z-axis solution gives 
good results for the propagation phase constant of the 
strip-type artificial dielectric medium. For larger values 


_of the spacing s, a solution based on modes propagating 


normal to the broad face of the strips gives good results. 
A combination of the two methods thus provides a 
method of estimating the value of the propagation con- 
stant without having to resort to the laborious proce- 
dure of taking higher-order mode interaction into ac- 
count. In order to estimate the best value of the prop- 
agation phase constant it is necessary to plot Bm, as a 
function of the strip spacing s with other parameters 
fixed. When this is done it is found that for a certain 
range of s, the two methods give results that are in 
close agreement. For smaller values of s the best es- 
timate curve of Bn. should be drawn so as to approach 
the z-axis solution, while for large values of s the best es- 
timate curve for 82, should be made to approach the 
x-axis solution. Typical results that are obtained by this 
procedure are given in Collin.‘ 

From a knowledge of the propagation constant 6n,, 
an equivalent dielectric constant for the medium, along 
the z and y axis, may be defined and evaluated by con- 
sidering an analogy between the strip-type artificial di- 
electric and a homogeneous uniaxial anisotropic di- 
electric medium. Some of the difficulties encountered in 
formulating such an analogy were pointed out. In par- 
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ticular, it was shown that a complete analogy with a 
homogeneous anistropic dielectric medium was not 
possible in general. 


APPENDIX 


The strip-medium could be assumed to have a dielec- 
tric constant of unity in both the x and y directions. In 
this case one would attempt to make an analogy with a 
uniaxial medium with dielectric constants 1, 1, x, along 
the x, y and gz axis. The z-axis now becomes the optic 
axis. In this case the eigenvalue equation for the homo- 
geneous dielectric medium is 


2 
Keke 
By elv0 


ae N2ke + (1 — n,”) ; 


Only if x,y=0°so n?2=1-—n, is this equation the same 
as (5). Since (5) arises naturally from the low-frequency 
limit of (1) it seems preferable to use the analogy dis- 
cussed in the paper. 

A further point of interest is that if the analogy dis- 
cussed in this Appendix is used then in place of (6), the 
dielectric constant is in general given by 


_ (na)? + (Bry)? 
ko? ni (Bn)? 


Ke 


which is not the same as that given by (6). The analogy 
discussed in this Appendix is the same as that used by 
Kolettis.'° 

Although the two analogies lead to different results 
for the equivalent dielectric constant x,., they are both 
of equal usefulness in practice since, invariably, the 
parameter xk, is used to compute 6 and, provided the ap- 
propriate eigenvalue equation is used, the same value 
of 6 obviously results. Clearly, since the strip medium 
is not fully analogous to a uniaxial homogeneous aniso- 
tropic dielectric medium, not too much physical signif- 
icance can be attached to the defined equivalent dielec- 
tric constant in either case. 


10 N. J. Kolettis, “Electric Anisotropic Properties of the Metallic- 
Strip-Type Periodic Medium,” Case Inst. Tech., Cleveland, Ohio, 
Scientific Rept. No. 17 issued on Contract No. AF 19-(604)-3887; 
October, 1960. 
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Optimization of Waveguide Tapers Capable 


of Multimode Propagation* 
OG, 8b WANG: 


Summary—By converting Maxwell’s equations, the general case 
of mode conversion in tapered waveguides is treated by matrix 
formulation in terms of an infinite set of coupled differential equa- 
tions with nonuniform coupling coefficients and varying phase con- 
stants. An “orthogonalization” or “diagonalization” process is intro- 
duced through a nonlinear matrix transformation which is a function 
of taper length. The general matrix solution of the problem is obtained 
through a perturbation method in the form of an integral equation 
of the Volterra type, and the integral equation is solved by an itera- 
tion method. In view of the difficulties in finding eigenvalues, the 
problem is then reduced to the two-mode case, and the mode con- 
version is obtained in an explicit form revealing certain information 
which characterizes the choice of “mode-conversion distribution 
function.” Optimization is first obtained through proper choice of 
the mode-converison distribution function. In an attempt to approxi- 
mate a Tchebycheff mode-conversion response, further optimization 
is realized by creating “new zeros” and thereby changing the density 
of the distribution of zeros in the vicinity of the origin of the mode- 
conversion curve and the nature of the optimization procedure es- 
sentially becomes that of synthesis. Through using the optimized dis- 
tribution function, a total reduction of about 50 per cent in taper 
length is realized (when compared with the cosine-squared dis- 
tribution) for the case of 50-db prescribed-mode discrimination 
in a taper connecting a j-in ID waveguide to a 2-in ID waveguide 
operating in the circular electric mode up to 75 kmc. 


I. INTRODUCTION 


PERUSAL of recent literature shows that the 
AN tapered waveguide capable of multimode propa- 

gation has been a subject of interest for the past 
few years, since in long-distance transmission by use of 
the low-loss circular electric mode (TEo1) in circular 
waveguide, tapers are necessary in several important 
applications. Tanaka,! expanding the field of conical 
guide in terms of the eigenfunctions of the uniform cir- 
cular guide and then matching the fields at the junc- 
tion, obtained a general expression for a mode conver- 
sion through the taper. Solymar,? using a similar tech- 
nique with some approximations, gave the design proce- 
dure of a one-section conical taper and a multisection 
conical taper. Savvirykh,? combining perturbation 
theory with the method of W.K.B., treated the field in 
the tapered waveguide in terms of an eigenfunction 
series expansion in an “artificial” orthogonal coordinate 
system. Unger‘ obtained an improved design of the cir- 


* Received by the PGMTT, April 17, 1961; revised manuscript 
received, June 23, 1961. 
t Bell Telephone Labs., Murray Hill, N. J. 

_ 1K. Tanaka, “Mode Conversion Through the Tapered Section of 
Circular Waveguide System,” presented at the Congrés Internat. 
Circuits et Antennas Hyperfrequences, Paris, France; October, 1957, 

2 L. Solymar, “Monotonic multisection tapers for over-moded 
circular waveguides,” Proc. IEE, vol. 106, pt. B, suppl. no. 13, pp. 
121, 1959. 

8S. K. Savvirykh, “On the theory of tapered circular waveguides,” 
Radioteckh. i Elektron, vol. 4, pp. 972-1002; April, 1959. 

“H. Unger, “Circular waveguide taper of improved design,” 
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cular waveguide taper with varying cone angle by using 
two coupled differential equations. The present paper 
treats the general case of an infinite set of coupled 
differential equations in matrix formulation, and the 
general matrix solution of the problem is obtained 
through solving an integral equation of the Volterra 
type by an iteration method. A procedure to create 
“new zeros” in mode-conversion characteristics is in- 
troduced in order to “optimize” the taper for “mini- 

mum” length. 

The presence of a taper in a waveguide inevitably 
introduces spurious modes. On the assumption that the 
taper possesses perfect symmetry and its axis is per- 
fectly straight, only TE», modes will be excited in the 
tapered region. If the taper axis has a slight curvature, 
additional TM, modes will be excited due to the de- 
generacy between the TE), and TM:., modes. Our 
goal in designing the taper is to reduce the spurious 
modes to a prescribed level in the operating frequency 
range with a taper length as short as possible. 

To describe electromagnetic fields in a perfectly con- 
ducting curved waveguide with nonuniform and arbi- 
trary cross section, we have to solve Maxwell’s equa- 
tions with appropriate boundary conditions. Schel- 
kunoff® has shown that certain field problems with com- 
plicated boundary conditions can be handled more 
readily by converting Maxwell’s equations into gen- 
eralized telegraphist’s equations. To convert Maxwell's 
equations into telegraphist’s equations, we introduce a 
complete set of orthogonal functions most appropriate 
to the particular geometry of the problem in question. 
Fields within the guide are then expanded in a series 
in terms of these orthogonal functions of a complete 
set. Subsituting these fields into Maxwell’s equations 
with due care given to the convergence property of the 
series on the boundary and taking advantage of the 
particular orthogonality relations of the functions, we 
obtain the following generalized telegraphist’s equa- 
tions: 


dVm - — 

= | = Lens E mnrn 
dw > ~ ae 
dln 


> 25 VimnaVn ate De ST anls 


dw R= n= 


(Mess REP hice nice cu (1) 


i Se peemeett “Conversion of Maxwell’s equations into gen- 
eralized telegraphist’s equations,” Bell Sys. Tech. J. : 
995-1403; September, 1955. TG ee ee 
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where w is the coordinate along the assumed direction 
of propagation, Zn, is the self or mutual series im- 
pedances, Y,,, is the self or mutual shunt admittances, 
*Tmn is the “voltage-transfer coefficients,” and !Ty,» is 
‘the “current-transfer coefficients.” Vm and Im are, .re- 
spectively, the voltages and currents related to the 
amplitudes of electric and magnetic intensities asso- 
ciated with each particular function. For the circular- 
waveguide taper, we choose as the complete set the 
orthogonal modes of a straight-circular waveguide with 
perfectly conducting walls and filled with a homo- 
geneous dielectric as in Fig. 1. It is obvious that the 


? a 
l 4 Sal 


Fig. 1—Circular waveguide taper. 


choice is arbitrary, although most appropriate in this 
case. The infinite set of first-order differential equa- 
tions (1) with coupling among all possible modes re- 
duces, in this case, to an infinite set with coupling 
among TE», modes only and takes the following form‘ 
in cylindrical cordinates: 


fae out, 1 8S, Bhnbe_ 
lg De 
is Len pect Rees ae 
cy 1 da 2, 2kmkn 
fo ETS SCG ee ieetay a 


Wek ae (2) 


where the summations are extended over all 1 ex- 
cept 2=M. Ym is the propagation constant of the mth 
mode, and k,, is the mth root of the Bessel function J,. 
If we let A, and R,, be the amplitudes of the forward 
and backward waves of a typical mode TEom, the fol- 
lowing equations are always true: 


Vee —= Rie a == Ke.) 5 


1 
Thess SS = Am ihe: ; 
Va. (3) 


where Z,, is the wave impedance of TEom mode 


_ jew 
Ym 


Substitution of (3) into (2) results in a new set of equa- 
tions in A,, and R,,. If the taper is very gradual, we can 
assume backward waves and multiple reflections neg- 
ligibly small and obtain the following infinite set of 
coupled differential equations in forward-wave ampli- 
tudes: 
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dA, ES 
et Ss amin or : cans, 
dz a (4) 


where £,,,s are the coupling coefficients defined by 


AEE oe Zn Zh 2 
a Rel 


II. GENERAL FORMULATION IN MaTRIX 
REPRESENTATION 


€ 
Smn 


The infinite set of equations (4) for a section of a 
tapered waveguide can be conveniently cast in matrix 
form 


D{ Ay =[M]{4}, {4@} = {c}, (6) 


where D, { } and [ | represent, respectively, an op- 
erator, a column matrix and a square matrix as follows: 


A,(z) 
ieee = ae) 
dz A3(z) 
—71(z) — Ero(z) — Exs(z) > - 
[] ee fo1(z) —o(z) foa(z) ++ - 
— r3(2) Sipeuks 


£a1(2) Es2(2) 


[M(z)] is assumed to be bounded continuous square- 
matrix function in a certain interval of the z axis. 
{C} is a column matrix of constants. It is seen that the 
presence of tapering introduces the coupling terms 
Emn(z) in [M], and the 4,, is no longer orthogonal. In a 
uniform circular waveguide, the coupling coefficients 
Emn Vanish, the propagation constants Ym are independ- 
ent of g, and the amplitudes A, form a complete 
orthogonal set. 

For the lossless case, we can show with the aid of the 
law of conservation of energy that the matrix [M | is, in 
general, complex and &nn=—&nm*.6 In our particular 
case, the elements £m, of the matrix are all real, and 
Emn = —£nm since the impedances Z,, and Z,, are all real 
for all modes far away from the cut-off as seen from 
(5). The propagation constant y,, is purely imaginary 
for perfectly conducting tapers and has the form 
Y=jBm, where Bm is real. The square matrix [IM] then 
takes the form 


—jB.(2) C12(2) Cis(z) - + - 
—Cix(2) —jB2(z)  Cas(z) + - - 


Le] it — Cis(2) — Cos(2) —JB3(2) Aig! 


where Emn = Cmn, real functions of z. 


6 The star * implies the complex conjugate. 
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We now introduce a nonlinear, nonsingular matrix 
transformation of the following form: 


{A} = [P®@]{B@}, (8) 


in an attempt to obtain a new orthogonal set of ampli- 
tudes B,,(z) along the tapered waveguide in at least a 
localwise sense. Substituting (8) into (6) and following 
the rules of matrix calculus, we obtain the following: 


[P]D{ B} = [m][P]{B} — D[P]{B}. (9) 
Multiplication of (9) by [P| yields 
D{B} = [El{B}, {BO} =[POPNC}, 


where 


(10) 


and 


Clie a reals ac fai Bed onze bah 


If the transformation matrix [P] is independent of z, 
the second term in (11) would vanish and the matrix 
[Q] is truly diagonal and the B,,’s would be orthogonal 
as seen from (10). In general, the matrix [Q] is not 
diagonal because of the presence of the term [P|“D[P], 
which represents the effect of tapering. However, if 
[P]“D|P] can be considered as a small perturbing 
term to the uniform guide, the matrix differential equa- 
tion (10) can be solved approximately by perturbation 
methods. It is worthwhile to note that the transformed- 
matrix differential equation (10) is of the same form as 
the original differential equation (6); however, the sig- 
nificance lies in the fact that the weights of the non- 
diagonal terms of the matrix |] have been shifted into 
the diagonal terms of the matrix [Q] by the transforma- 
tion of (8). Eq. (10), therefore, can yield a more accu- 
rate solution by the perturbation method than (6). 

To find the matrix [Q], we must first obtain the 
eigenvectors comprising the diagonalizing matrix [P]. 
This involves the determination of eigenvalues from 
the secular equation of the matrix 


| [4] — arlz]| =, (12) 


where [J] is the unit matrix. The values of \ for which 
the equation is satisfied are the desired eigenvalues. 
Since [M] is, in general, complex and nonsymmetrical, 
the eigenvalues and the associated eigenvectors will 
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also, in general, be complex. From (11) and (8) we see 
that the required transformation matrix [P(z)| also 
performs the role of a similarity transformation to the 
matrix [M(z)]. If the eigenvalues of the matrix [M(z) | 
are distinct, there exists a matrix [P(z)] that diagonal- 
izes the matrix [M(z)]. The eigenvectors belonging to 
the corresponding distinct eigenvalues then form an 
orthonormal basis, 2.e., 
Dy Pun mea= One, 


m 


NS = 1,0),-3) ae een 
where * denotes the complex conjugate and 6,, the 
Kronecker delta symbol. When [P] is found, the evalu- 
ation of [P]~D|P| is straightforward. Our task now is 
to solve the infinite set of coupled differential equations 
(10) with the appropriate boundary conditions. 

The matrix integral of the system of (10) is asquare 
matrix [B(z)|, the columns of which are n-linearly- 
independent solutions of the system. Since each column 
of the matrix [B(z)] satisfies (10), the matrix integral 
[B(z) | also satisfies the equation 


D[B] = [Q][B], [BO] = [1]; 


where [J] is the identity matrix. The formal solution of 
(10) can now be written in the form 


{B= [B\ (Poy fe}. (15) 


We now seek the solution of system (14)? which in 
view of (11) can be written as 


D[B] = [q][B] — [e][B], [BO] = UI. 


Under the condition that |e] may be considered as a 
very small perturbing matrix relative to the diagonal 
matrix [gq], we can take [e][B] as the perturbing non- 
homogeneous term of the matrix differential equation 
(16). The approximate solution of (16) in matrix repre- 
sentation is 


[B] = exp 1 J tales — {ex | f tales 
exp q- f “[a(e)]a't [e(e")][B(e) Jae’ 


Eq. (17) is a Volterra integral equation of the second 
kind and can be solved by an iteration method. We ob- 
tain, therefore, the final solution of (16) as an infinite 
series in powers of [e| 


(14) 


(16) 


(17) 


is] = ex { f Tales} $2) — [exp {- f ‘Tate [dedlexp { f é fae) |as'h ag 
+ [ [ex {- if lac |ae [«@)] exp | {i “La(eD]a't | 
Ap is {- if “[a(e"y]as"t [se")] exp | J taerntaert da") dg 


+ terms of higher powers in (a ; 


(18) 


7 See Appendix I. 
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The uniform convergence of the series in (18) in the 
closed interval can be shown according to standard 
methods.® The final matrix solution {A(z)} is obtained 
through (18), (15), and (8). 

The above general formulation can, in principle, be 
applied to any number of coupled modes; however, the 
explicit solutions of the eigenvalues of (12) are already 
too clumsy to handle even in the case of three coupled 
modes. Accordingly, only the solution for two coupled 
modes (TEm and TE), which have the strongest 
coupling, will be carried out in the following equation 
for a gradual taper. Inspection of the coupling coeffi- 
cient in (5) justifies the preceding statement in addi- 
tion to the fact that the phase-constant difference is 
much larger for all other higher-order modes. For TEo: 
and TE 2 coupling, (12) yields the two eigenvalues 


Xi = 25[—(61 + 6.) + V/(Gr— Bo)? + 4012? 


eae] (Se) Ie) (19) 
and 
ho = 4j[—(81 + B2) — V(Bi — Bx)? + 4C12?| 
sie 1); 
where 
T= VAs? + Ci? (20) 
and 


Ap = 3(B1 — Bo); B= 3(B1 + B:). 


2 2 ay 
~jexp jsf @+ roast f 
0 0 dz 
| -exp {2 f r(e)as') de 
0 


The transformation matrix [P] then takes the form 


) Wie | (21) 
| i S| 


We see that the unitary condition of (13) is satisfied 
by the [P] matrix, and its determinant is —1 in agree- 
ment with the fact that the determinant of an orthonor- 
mal matrix can only be +1 (corresponding to rotation) 


| 
| 
| 


eh 


8F. G. Tricomi, “Integral Equations,” Interscience Publishers, 
Inc., New York, N. Y., p. 10; 1957. 
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or —1 (corresponding to reflection). If use is made of 
(20) by letting the angle between the sides A8 and T' be 
(28), so that cos 20=A8/T, then the matrix [P] of (21) 
takes the simple form 


[P] = ee coe 

sin@ 7cosé 
For small angle (26), the variation of (20) is directly pro- 
portional to that of the coupling coefficient Cj, and in- 
versely proportional to that of the phase-constant dif- 
ference AG along the taper. Accordingly, this variation 
of (26) can appropriately be interpreted as “mode-con- 


version distribution” along the taper. The matrix [e] of 
(11) in terms of 0 becomes a symmetric matrix 


do 
dz 


(22) 


[e] = [P]}*D[P] = -j (23) 


dé 

dz 
It is clear now that the coupling to the TEo: mode be- 
comes much smaller for higher-order modes. For a 


gradual taper, we see that [e| is, indeed, a small per- 
turbing term to the diagonal matrix 


l= (Pbael=-[" 5° 4, | 28 


Note that the trace of the diagonalized matrix remains 
equal to ($:+ 2). For a solution to the first order of 
approximation, we obtain from (18) via (23) and (24) 


~jexp if e- rash f = 


“exp \-4 f rast dz’ 
0 


exp \—j [e+ nya 


(25) 


The normalized boundary condition for the two-mode 
case requires that 


|{4@}] = |{Bs@}| = |[POr{c}| 


[} 


where the vertical bars denote the norm of a vector of 
complex elements. For a gentle, smooth taper with a 
mode-conversion distribution function vanishing at 
both ends of the taper, we obtain the explicit solution 
for A,(1) from (8), (22), (15), (25), and (26) 


I I dé 
A) = exp {jf @+ nak fo 
“exp {24 f rcyae'k dz. 


(26) 


(27) 
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Note that | 42(2)| is essentially in the form of a Fourier 
transform of d6/dz. Letting 


p@) = [ Peas (28) 
0 
and integrating (27) by parts, we have 
pl ; 
| Ax(on| = | f (20d. (29) 
0 


If we integrate (27) by parts in another way, the fol- 
lowing expression is obtained: 


| 42(0)| = Salen Ge “i Be 
sf oe 
tages Gh Ea a 
ea ae 


a= = [ Teds. 


(30) 
where 


(31) 


Eq. (29) is in a form suitable for the computation of 
“mode discrimination” when the “conversion-distribu- 
tion function” (26) is given in terms of the parameter p. 
Eq. (30), on the other hand, gives us a clue that in 
order to obtain a higher-mode discrimination, it is ad- 
visable to choose a conversion-distribution function with 
vanishing first and higher derivatives at both taper 
ends. A detailed discussion will be given in Section V 
in this respect. When the distribution function (26) is 
chosen, the waveguide radius a(p) can be obtained from 
(5) and (20) and the taper length of the guide z(p) from 
(28). 


III. THe CHoIcE of CONVERSION DISTRIBUTION 
FUNCTION 


It is evident that the choice of distribution function 
is not unique. Under the stipulation that the function 
itself vanishes at both ends as expressed in (8), (22) 
and (26), a simple choice of such a function in the form 
of an infinite series in (p/p) is 


TT 
20 = K, sin” (“), 
Pi 


where sin (7p/p:) can be considered as a “generating” 
function, and m may or may not be an integer. Substi- 
tuting (32) into (29) and using (47) in Appendix II, we 
can show that the mode discrimination is given by 


(32) 
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Fig. 2—(a)—(c) Mode conversion in waveguide tapers. 
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where 


. T 2Rmen (0) 
a = 2), bi and C= LOS 
P1 Ri coca ae ay 


(34) 


‘and k,, and k, are defined in (2). It is evident from (33) 
that C corresponds to the mode conversion of a step 
discontinuity in the diameter of the waveguide. For 
cases where ” is noninteger, general solutions for the 
mode discrimination in closed form are impossible, and 
we have to resort to numerical integration. Eq. (33) for 
integral values of 1 are plotted in Fig. 2(a)—2(c). 

It is seen from these figures that there is always an 
“optimum” integer which minimizes the length of a 
taper for a prescribed discrimination. Alternatively, if 
we prescribe the taper length, there is always an opti- 
mum integer 2 that provides the highest-mode discrim- 
“ination. From Fig. 2(a) we see that, for a fixed taper 
length of pi=18.5, optimum discrimination occurs for 
-n=7. On the other hand, if we prescribe a —50-db 
mode discrimination within a frequency range up to 75 
kMc, the required minimum length of the taper for the 
case (a(0)=7"/16 and a(/)=1”") is about 3 ft* using 


OPTIMIZED n=3 


0.9 


0.8 


RADIUS a IN INCH 
te} 
N 


0.5 


ieee eee all 2S Sees Se 
O40 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 
ABSOLUTE LENGTH z IN INCHES 
Fig. 3—Comparison of profile of tapers of same 
mode conversion (50 db). 
3430 
z 


lAal 
Gu 


as 


p 
lAal «fe Gn(42)"(-5-Y'dp 


MINATION 
® 
fe) 


SCRI 
8 


Ol 
4 


MOD! 


lAal 
SCRIMINA TON, 


MODE Di 


3 
Fig. 4—Mode conversion in waveguide tapers. 


Tang: Optimization of Waveguide Tapers Capable of Multimode Propagation 


447 


n=2, but only about 2 ft using m=3 as shown in Fig. 3. 
The actual computation of radius a and length zg is 
shown in Appendix II. 

Another simple choice for the distribution function is 
a polynomial in (p/p1) 


4 n n 
09 (0-2). 
Pi Pl 


where the “generating” function is (49/p1)(1—p/p1), 
and m may or may not be an integer. Evaluating the 
integral of (29) for (35) with integer m, it can be shown 
that the mode discrimination is equal to 


SA 2 2 Pi 
| 42(p1) | = 5 ( ) (2: (e7! + 1) 
RPL ON Pi 


2 
tees (e771 — »| tomas 
Ai 


(35) 


| ta(o,) | (=) [(F +5) 1) 
As = eel — 
(p1 ! 23- py pi a a? 
12 | 
— Een + »| TOL Ze 
a‘ | 
Ser (DN Oph /aes 00) 
| Ao(p1) | a ( (o <r Sand CA + 1) 
Deo Non, at a? 


(36) 


forn = 3, 


72 10 | 
ae hear (r+ (ear! es 1) 
a? a | 


where C is defined in (34). Eq. (36) is plotted in Fig. 4. 
It can be shown that both (33) and (36) have the value 
unity times C at p1=0, as they should. Eqs. (33) have all 
their zeros at multiples of 7/2; on the other hand, (36) 
has its initial zeros shifted closer toward the origin pi=0, 
as shown in Fig. 4. Inspection of Fig. 4 shows that the 
zeros at large values of p; also gradually shift to positions 
at multiples of r/2 as (33) do. Comparison of Fig. 2 and 
Fig. 4 for curves of corresponding values of 7 shows that 
the sine distribution has better over-all discrimination 
than the polynomial distribution, except in the region 
between p; =0 and the first zero of (33). 


IV. OPTIMIZATION 


In an attempt to further generalize and optimize the 
mode-conversion distribution function, we expand it in 
a symmetrical Fourier series (with the origin of the co- 
ordinate system shifted to the center of the taper) 


(37) 


co) Tp 
26-= >> Dm cosm(™), 


m=0 Pl 


With suitably-chosen coefficients D,,’s, it is obvious that 
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(37) can include all cases described by (32) and (39). 
Substituting (37) into (29), we have 


ype 2 Tp 
| A 2(pi) | = f De D, cosm (=) cindy 
Pilg m=0 Pl 


TT 
2D» cos m — (e%! — 1) 


| Aa(pr) | 7 »y ‘ , 


9 9L9 
m=0 a” ae mb” 


meven (38a) 


Tv 
mb D, sinm 5 (e271 + 1) 


, modd (38b) 


| 42(o1) | 7 NE 


m=1 
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where a and 6 are defined in (34). It is important to 
note that only odd or even values of m are required to 
represent the mode-conversion distribution function 
given in (33), depending on whether is odd or even. 
Equating (38) to (33) for the corresponding case, we 
can obtain the D,, in (38) as shown in Table I. It is in- 
teresting to note that the ratio of D,,’s for a particular 
n correspond to the coefficients of the binomial expan- 
sion. We are now in a position to further optimize the 
mode-conversion curve defined by (33), using the values 
of Table I as a guide. The aim is to reshape the mode- 
conversion curve in such a way that the first few maxima 
of the “side lobes” of (33) will be leveled and at the 
same time lowered optimally. Before carrying out the 
optimizing procedure, it is appropriate at this point to 
discuss the incentive for this further optimization in 
more detail in the following paragraph. 


TABLE I 
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The class of distribution functions that we have been 
considering have the property that the distribution 
function and all its derivatives are single-valued, uni- 
formly bounded, and continuous in the interval of in- 
terest. Any distribution function of this class can be 
transformed to a general form in terms of the zeros of 
the function. With the coordinate origin at the center 
of the taper, the function of interest has the form 


26 = f(p) E = (2) | ‘ 


where f(p) is an even function due to the symmetry of 
the function (26), and f(p) does not vanish at the taper 
ends (9 = +p;/2) or any other value of p in the interval. 
Taylor? has shown that the Fourier transform of (39) 
has the following asymptotic form as p; approaches 


infinity: 
( n+ 1 ) 
cos { pi — re 
ree 


(pi) 


(39) 


(40) 


F(p:) ~ (#) T(x + 1) 


It is, therefore, seen that the mode discrimination of 
very long tapers is only trivially different, no matter 
what form the distribution function has. On the other 
hand, as we have seen earlier, the initial slope of the 
function and accordingly the value of are of consid- 
erable importance in “optimizing” the taper. The initial 
slope of the function (2@) determines the positions of the 
zeros near the origin of the mode discrimination curve 
A,(p). The smaller the initial slope, the further will the 
first zero be from the origin. Likewise, the value of 
alters the zero positions, since ” changes the slope of 
the function. Eq. (40) indicates clearly that spurious 
mode conversion or discrimination decays as 1/| p:| ar 
and zeros appear at p:=(7/2) at large values of py. 
Inspection of (33), (36) and (38) confirms this decay 
rate and the position of zeros at large values of p;. At 
this point, it is particularly appropriate to compare 
this decay rate with the nondecaying characteristics of 
the Tchebycheff polynomial of infinite degree. If we are 
to simulate the Tchebycheff polynomial of infinite de- 
gree by the function of (40), we see that the only choice 
to make (40) nondecaying is to make = —1. However, 
for this choice of m, the function of (39) will have poles 
at the taper ends and will no longer be uniformly 
bounded. This violates our basic requirement, and, 


therefore, it is clear that a smooth transition taper with 


its mode-conversion characteristics described by a 
Tchebycheff polynomial of infinite degree is unrealiza- 
ble. In fact, this unrealizability is simply a consequence 
of the law of conservation of energy. This is why even 


_°? T. T. Taylor, “Design of line-source antenna for narrow band- 
width and low side lobes,” IRE TRaNs, ON ANTENNAS AND PROPAGA- 
TION, vol. AP-3, pp. 16-28; January, 1955. 
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in the transmission line case either steps!” have to be 
introduced at taper ends or a modified!! Tchebycheff 
polynomial of infinite degree has to be used in the syn- 
thesis procedure. Although the frequency range of the 


design using infinite degree Tchebycheff polynomial ex- 


tends to infinity, for band-pass applications this is un- 
necessary. Accordingly, the Tchebycheff design is opti- 
mum only in the sense of infinite bandwidth, even for a 
transmission-line taper. 

From the above exposition, we see that in order to 
further optimize, we need to “flatten” or “level” the 
decay rate of the first few “minor lobes” of the mode- 
conversion curve as much as permissible, after choos- 
ing the “optimum” value of m for a prescribed discrim- 
ination level from Fig. 2(c). It is clear that the nature of 
this procedure essentially becomes that of synthesis. 
At first sight, this might look rather aimless if we do it 
in a heuristic way by adjusting the coefficients D,, of 
(38). Knowing, however, that the shape of the mode- 
conversion curve of (38) depends very much on the 
density of the distribution of zeros near the origin, we 
are thus led to create an extra “zero” at such a position 
that the “side lobes” near the origin will be leveled and 
at the same time optimally lowered. It is important to 
note that the new zeros will be introduced by properly 
choosing the coefficients D,, in (38), while the original 
zeros were determined only by the term (e771+1). This 
will be done first for 7 =3. With reference to (38b), we 
have 


m[e-(o2)] 


N2 
as 3D3 fe Fal (o, =) | zs 0, (41) 
T 
where 
bc aC F 
De, as Da ——i Dy, . 
2p 


(42) 


We can solve for the required coefficients D,’ and Dx 
when the position of the new zero is intelligently se- 
lected. Inspection of Fig. 2(a) for »=3 shows that it is 
advisable to locate the “new zero” around p,=6(7/2) 
in order to achieve the desired results. With this value 
of pi, we get Dy’ =1.09375 and D;’ =0.28125 from (41) 


10 R. W. Klopfenstein, “A transmission line taper of improved 
design,” Proc. IRE, vol. 44, pp. 31-35; January, 1956. | : 

1 R. E. Collin, “The optimum tapered transmission line matching 
section,” Proc. IRE, vol. 44, pp. 539-548; April, 1956. 
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and (42). A plot of (38) with these values of D,’ and 
D;' as shown in Fig. 5 indeed gives the desired results. 
Fig. 6, showing the relation between the position of the 
hew zero and the maximum of the maxima of the side 
lobes, also confirms the fact that the optimum location 
of the new zero should be in the vicinity of 6(/2). 
Mode-conversion curves for 


3 
cos? () and E (1 — *)] 
Pl Pi P1 


are also plotted in Fig. 5 for comparison. It is seen that 
the optimized conversion curve does have its first few 
side lobes maxima leveled, and the maxima of the re- 
maining side lobes decay according to a rate asymptoti- 
cally proportional to 1/p;4. The improvement is about 
another 30 per cent reduction in taper length for 50-db 
discrimination. The actual optimized length of the 
taper is plotted in Fig. 3 for comparison. Through us- 
ing the optimized cos* (7p/pi) distribution function, a 
total reduction of about 50 per cent in taper length is 
realized for the case of 50-db prescribed mode discrim- 
ination in a taper connecting a 4-in ID waveguide to 
a 2-in ID waveguide. 
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Fig. 5—Comparison of mode conversion in waveguide tapers, 
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The same procedure can be applied, for example, to 
the n=5 case and we have 


2\2 TNE 
Tv Tv 
and 
IDE Ds 
| De meee =| = (44) 
3 5 
where 
aC 
JD, = ID! 
2p. 


Inspection of Fig. 2(a) shows that the logical choice for 
the location of the new zero is at p1=8(7/2). In this 
case, however, we have only two equations for three 
unknowns. It is necessary to assume an appropriate 
value of D,’ so that D;’ and D;’ can be determined. The 
first suitable choice of D,’ might be D,’=75/64 ob- 
tained from Table I for the m=5 case. Calculation shows 
again that the new zero should indeed be around 
p:=8(2/2), and the optimum choice of D, is indeed 
(75/64). This case is shown in Fig. 7 together with the 
cos® (1p/p1) case for comparison. 

Further investigation of the 7=5 case reveals a bet- 
ter value for Di’ because we can now create two new 
zeros in view of the extra undetermined coefficient. 
The assignment of the second zero will give us an extra 
equation and thus determine the three coefficients 
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Fig. 7—Comparison of mode conversion in waveguide tapers. 
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uniquely. A judicial choice of two zeros at pi = 8(7/2) 
and p:=10(m/2), respectively, should yield even better 
mode discrimination at a fixed taper length than the 
above case. This result is also plotted in Fig. 7 for com- 
parison and shows that the prediction is valid. A reloca- 
tion of the first “created zero” at pi=7.72(4/2) gives 
the best discrimination for this case, as shown in Fig. 7. 
Accordingly, it is seen that we can further optimize a 
taper by creating new zeros near the first few zeros of 
the mode-conversion curve, and the number of new 
zeros allowed to be created increases directly as the 
number of undetermined coefficients. On the other hand, 
the freedom to have more created zeros is only available 
for higher values of ~ which dictates higher mode dis- 
crimination and longer taper lengths. 


V. CONCLUSIONS 


In an attempt to approximate a Tchebycheff mode- 
conversion response in a wideband waveguide taper, we 
optimize the taper by creating new zeros in the mode- 
conversion response. This response had been initially se- 
lected to yield the “shortest” taper length in the “sine 
or cosine distribution function family” at a prescribed 
level of mode discrimination. 

We note that the first and second derivative of a dis- 
tribution function with ratio D3;/D,=4 [this is the cos* 
(xp/pi) case] vanish at the taper ends, but those with 
ratios other than } do not vanish. Thus, change of end 
slopes evidently is a consequence of the addition of new 
zeros in the mode-conversion curve, since the end slopes 
are closely related to the position and density of the zero 
distribution near the origin. Because the optimized dis- 
tribution functions do not have vanishing derivatives 
at the taper ends, it is seen that the distribution func- 
tions with vanishing derivatives at the taper ends may 
not be most desirable. We now return to (30) for further 
information in this respect. It can be easily shown that 
for the symmetrical distribution functions we used, the 
nth derivatives at the taper ends decreases as 1/p,"+, 
where pi is the equivalent taper length. This implies 
that longer tapers have smaller derivatives and that the 
(n+-1)th derivative is smaller than the mth derivative. 
Investigation of (30) also shows that if a symmetrical 
distribution function has its first and second derivative 
vanishing, the third derivative will represent a signifi- 
cant part of the mode conversion. On the other hand, if 
the first two derivatives are nonvanishing but very 
small, the total mode conversion may still be smaller 
than the vanishing derivative case due to alternate 
+ signs and — signs in the real and imaginary parts of 
(30). 

From what has been shown here, it is obvious that we 
cannot claim to have synthesized the absolute optimum 
taper, but we can claim the taper to be very close to 
the absolute optimum. The procedure discussed above 
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may lead to a further optimization, but it is quite clear 
that the small reduction of taper length due to such a 
procedure might not justify the amount of computa- 
tional labor involved. The starting point is still (32). For 
instance, for a certain prescribed spurious mode level, 
there is in (32) an “optimum” n which may be non- 
integer. The evaluation of (29) with a distribution func- 
tion of noninteger power will require numerical integra- 
tion, and the evaluation of the radius and length of the 
taper as shown in Appendix II will again require numer- 
ical integration. 


APPENDIX I| 


If the matrix [6] can be considered as a perturbing 
matrix to a matrix [N], the matrix exponential e!¥)+0) 
can be expressed in several ways. If [N] and [6] com- 
mute, then 


e(N1+18] = pINIe[8). 


Considering the general case where [NV] and [6] do not 
commute, we can write 


e+] — 


By expanding the series, we can collect the terms in 
powers of the perturbing term [6] in the form of 


e(NI1+(8] = e(N] + ss [6]*f.([V’]). 


k=1 


This, however, cannot be accomplished in a neat fash- 
ion, and it is necessary to solve equations of the form 
' of (16) by perturbation techniques. 


APPENDIX II 


The coefficient Cy is obtained from (5) under the stip- 
ulation of the TEo; and TEo2 modes that are far away 
from cut-off, and it takes the following simple form: 


1 da  2kik» k da 
AS > = are oe a 
a dz ky — ky? a dz 
where 
2Rik2 
ay he 


From the preceding equation and (20) and (28), we have 


le da ip oe 
ay a 0 k 


for a gentle taper. 


(45) 


Tang: Optimization of Waveguide Tapers Capable of Multimode Propagation 
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The actual length of the taper in terms of the parameters 
pi is obtained from (28). For a gentle taper it is 


; fee are 
ia 46 fa tei 1) 


1 P dp 
Ro? ky4 f =|. 49 
18 (Rot ki") eee (49) 


where fo is the phase constant in free space. The last 
term in (49) is negligibly small, in general. Substitution 
of (48) into (49) gives 
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A Broad-Band Coaxial Ferrite 
Switch* 


A broad-band strip-line reflective ferrite 
switch has been described by Johnson and 
Wiltse,! who also referred to the possibility 
of a similar switch in coaxial line. This note 
describes a coaxial on-off switch which will 
operate over the band 2500-4100 Mc; two 
such units can be combined to make a two- 
way switch. An isolation of 40 db was 
achieved, with a very low loss for the trans- 
mitting path, which, in both these devices, 
was obtained by magnetizing the ferrite well 


i pecs by the PGMTT, June 5, 1961. 

M. Johnson and J. C. Wiltse, “A broad-band 
ene: reflective switch,” IRE TRANs. on MIcrRO- 
WAVE THEORY AND TECHNIQUES, vol. MTT-8, pp. 
466-467; July, 1960. 


beyond the value for isolation.” 

Each switching element employed a small 
slug of a developmental ferrite (Bsar= 2280 
gauss) which completely filled a half-inch 
section of air spaced coaxial line of 9/32 inch 
outer diameter and ¢ inch inner diameter, 
the ends of which oe directly coupled to 
Type “C” coaxial connectors. Fig. 1 shows 
the attenuation obtained with a field of 
about 400 oe compared with pads of approxi- 
mately 30 db and 40 db, while Fig. 2 
shows the attenuation obtained with a field 
of about 2500 oe compared with a 3-db pad. 
The attenuation of about 40 db was ob- 
tained at 100°C in a convection cooled sole- 


noid, but greater attenuation was achieved 


at lower temperatures. The VSWR under 
reflecting conditions in Fig. 3 is about 0.15, 


2The General Electric Company, Ltd., Brit. 
Patent Application No, 19948/59; June 10, 1959. 


but a lower value may be obtained over a 
smaller bandwidth. The solenoid power re- 
quired was rather high, but was reduced by 
constructing the coaxial line of iron with a 
brass section in the position of the ferrite 
slug. Permanent magnet bias can also be 
used without unduly slowing the switching 
speed. 

Fig. 4 shows the variation in the position 
of the effective short-circuit planes in front 
of a slug of a similar ferrite material. A two- 
way switch may be constructed therefore by 
arranging for a high impedance to appear at 
the T junction. Two ferrite slugs were used 
in the top of the “T,” one of which was 
magnetized for isolation and the other for 
transmission at any instant. The transmis- 
sion loss for the two-way switch is shown in 
Big. 9: 

The development of this switch was part 
of the work done under a contract for the 
Admiralty. 
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Fig. 2—Transmission loss of ferrite slug. 
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Fig. 4—Short-circuit planes produced under 
constant reflecting conditions. 


Fig. 5—Transmission loss of T junction. 
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Conditions for Maximum Power 
Transfer* 


It is sometimes of interest to ask for the 
conditions of maximum power transfer from 
a fixed source into a load constrained to vary 
over an arbitrary contour in the impedance 
plane. There exists a simple graphical solu- 
tion to this question as shown below. Con- 
sider the circuit shown in Fig. 1. 

If P is the power delivered to the load 
and Po=(E#2/4R,) is the available power 
from the source, then 

le 4r 
Se a ee (1) 
Py (r+1)?4+ (x + 2,)? 


giving 


Po 2 
E _ (2 P — 1) | + (x + x,)2 
ee 1) il, @ 
= (02> 2919) 12) 
iP 
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R Xx 2€, 
— ) r= —>) i 
IK R; Re 
Z ; Ls 
SRM BER 


Eq. (2) represents a family of circles in the 
z plane of radius »/(2Po/P—1)?—1, whose 
centers lie along the line x=—x, through 
the point z,* as shown in Fig. 2. 
If a is the distance from unity to the 
center of a circle along the line x=—wx,, the 
equation for the family of circles becomes 


[yr — (a+ 1))? + (@ + 4,)? = (a2 + 1)? — 1, 


where a=2(Po/P—1). 

To find the impedance for maximum 
power transfer to Z, on the plane in which z 
is drawn, strike a line parallel to the 7 axis 
through the point z.*=1—jx,. Move an arbi- 
trary distance a from unity along this line, 
and from this point draw a circle of radius 


+jX 


Fig. 


Zz plane 
k 


Fig. 3. 


* Received by the PGMTT, June 15, 1961. 
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V(a+1)—1. Repeat until the circle is just 
tangent to the given z curve. The point of 
tangency gives the impedance Z,, for maxt- 
mum transfer, while 

12 2 


Po a+ 2 
as shown in Fig. 3. 

Note that maximum power transfer does 
not occur at the point of closest approach 
OA 

CarL SHULMAN 
School of Technology 
Dept. of Elec. Engrg. 

The City College 
New York, N. Y. 


A Comment on the Scattering 
Matrix of Cascaded 2n-Ports* 


Epprecht' calculated the scattering ma- 
trix of two cascaded two-ports. Redheffer? 
does the same for the 2”-port using non- 
standard notation. This note will comment 
on the physical interpretation of the con- 
stituents of the resultant-scattering matrix. 
To use the notation of Fig. 1, the scattering 
matrix constituents are 


Su = Su! + S1'Su"(1 = Soo/S’”) Sar 

Sto = Syo! (1 — Si1'S20') 4S 19” 

Sor = Sar! (1 = S20/S11) S91" 

So = Soo!” + So1/So9/(1 — Sy1/"S 99’) 3S 19". (1) 


r 
| 
| 
l 
IE 
] 
| 
| 
= 


Fig. 1—S’ and S”’ are n Xn scattering matrices of the 
respective networks. S is the scattering matrix of 
the resultant network: 

a 
Saat 2 


= Sur’ =] iS Su” 
Sor’ Soo! Ser” 

The interpretation given to these formu- 
las is that Sy is the bilinear transformation 
of Sy through the single primed network, 
and S» is the bilinear transformation 
of S2'’ through the double primed network. 
Both of these results also follow from the 


St 


* Received by the PGMTT, June 21, 1961. This 
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5, pp. 169-173; 1959. 

2R. M. Redheffer, “Inequalities for a matrix Ric- 
cati equation:” J. Math. and Mech., vol. 3, pp. 349- 
367; May, 1959. 
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definition of scattering matrix terms on the 
basis of matched termination (7.e., if the out- 
put has a matched load, S=0, the input co- 
efficient of the double primed network is Su”, 
which is the output coefficient of the primed 
network). Sj and Sx are similarly inter- 
pretable, with the special case of bilaterally 
matched networks being the “star” multi- 
plication of Altschuler and Kahn.* 

It should also be noted that formulas (1) 
are valid when an n-port and an m-port are 
cascaded? (or interconnected). 

Dla RerStoci 

L. J. Kaplan 

Dept. Elec. Engrg. 
New York University 
New York, N. Y. 
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TECHNIQUES, vol. MTT-4, pp. 228-233; October, 
1956. 

4L. J. Kaplan. and D. J R. Stock, “A generaliza- 
tion of the matrix Riccati equation and the ‘Star’ 
multiplication of Redheffer,” J. Math. and Mech., 
vol._6; November, 1962. 


Use of Flow Graphs to Evaluate 
Mistermination Errors in Loss and 
Phase Measurements* 


The purpose of this note is to show how 
the signal flow graph technique illustrated 
by Hunton! leads quite naturally to an ex- 
pression for error due to mistermination 
when measuring insertion loss and phase. 

We start with the flow graph used by 
Hunton to represent the tandem connection 
of generator, network and load: 


b 

ees ay Sor 2 
Ty, Su Soo re 

by Sie a2 


With the aid of the nontouching loop rule to 
solve the graph, Hunton easily obtained the 
result 

CS Se 

E  1-YTySu—VrS22+T lr (Si1S22—S12S21) 


Very little extra work is needed to compute 
insertion loss and phase measurement errors 
due to mistermination, once the above equa- 
tion is availabie. 

E is the wave amplitude at the output 
port of the generator when terminated in a 
matched load Zo. If V, znd Z, represent the 
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NIOUES vol. MTT-8, pp. 206-212; March, 1960. 
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Thevenin generator voltage and impedance, 
then 


— Zo 
Zee a 
Since 
1+T, 
Lh , 
Peet. 
i= 1; 
R= Wr 
2 g 


The above equations together give 
Vy Soi(1—Ty) 
2 1-TpSu—PrS2+T Pz (SiS22—Si2S21) 


bo= 


From the flow graph we see that a2=bI'z- 
The total wave amplitude across the load is 
therefore 


Vz 
Vom as tbs 


| sat 1pa Ftp | 
1—T,Su—TrS2+T 1 (SuS22—Si2S21) 


Now the measured insertion ratio Rm is ob- 
tained by dividing the load voltage with net- 
work removed by the load voltage with net- 
work inserted. To remove the network, we 
set Sn, Sx equal to zero and Siz, Sx to unity. 
The result is 


“i = SoeTp—SuTg tT Pr (S1$22—S12S 21) : 


Rn 
(1—T,Tr)S21 


If the source and load were reflectionless 
(T,=lrzx=0), the corresponding insertion 
ratio Ry would be just 


Ro = : 
Soy 
Hence, the quotient 
Rm 
= z 
1—SooTz—Suly+T lz (S 1S22—S 12501) 
x i—Tjiz - 


provides the measurement error due to net- 
work mistermination. In the common case 
where [, and I'y, are <1, Q simplifies to 


Q~1-+A 


where A, the fractional error in nepers and 
radians, is given by 


A= — Sul, — Soot 
+ PyPr(1 + SirS22 — S12S21). 


For reciprocal structures, Sj is equal to 
So; these in turn are equal to the reciprocal 
of the design insertion ratio Ro. 

As an example of the application of the 
expression for A, consider the measurement 
of a network having | Su| =|Sx2| =0.3 (cor- 
responding to a VSWR of 1.85) and |Si| 
=|S»|=1. Then, if source and load were 
such that |I,|=|I'z|=0.02 (VSWR of 
1.04), we could expect maximum errors of 
0.11 db or 0.73 degrees, depending on the 
phases of the S’s and I’s. 

DANIEL LEED 
Bell Telephone Labs. 
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Francaise des Electroniciens et des Radio- 
électriciens. 
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A. E. Karbowiak was born in Warsaw, 
Poland, on March 1, 1923. He received the 
B.S. degree in 1950 and the Ph.D. degree in 
electrical engineer- 
ing in 1954 from the 
University of Lon- 
don, England. 

From 1951 to 
1954 he was at Uni- 
versity College, Lon- 
don, carrying out re- 
search on_ surface 
waves. In 1954 
he joined Standard 
Telecommunication 
Labs., where he was 
responsible for  re- 
search on various aspects of long-distance 
communication by waveguide. Since 1958 
he has been head of the Microwave Depart- 
ment. 

Dr. Karbowiak is an associate member 


of the IEE. 
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Nicholas J. Kolettis was born in Calim- 
nos, Greece, on February 15, 1933. He re- 
ceived the B.S.E.E. and M.S.E.E. degrees 
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in 1958 and 1961, respectively, from the Case 
Institute of Technology, Cleveland, Ohio. 
While attending graduate school, he 
was a Research As- 
Same Win were 
Microwave Group of 
the Electrical Engi- 
neering Department 
at that Institute. At 
present, he isa Mem- 
ber of the Technical 
Staff of the RCA 
Laboratories, Prince- 
ton, N. J., where he 
is engaged in research 
work on microwave 
applications of solid 


N. J. KoLettis 


state devices. 

Mr. Kolettis is a member of Eta Kappa 
Nuand Tau Beta Pi, and an associate member 
of Sigma Xi. 


George L. Matthaei (S’49—A’52—M’57), 
for a photograph and biography, please see 
page 675 of November, 1960, issue of these 
TRANSACTIONS. 
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Charles C. H. Tang was born in Shang- 
hai, China, on September 27, 1924. He re- 
ceived the B.S. degree in physics in 1946, 
from the University 
of Shanghai. He re- 
ceived the M.S. de- 
gree in electrical en- 
gineering in 1952, 
from Oklahoma State 
University, Stillwa- 
ter, Okla. He entered 
Harvard University, 
Cambridge, Mass., 
on a Gordon McKay 
Scholarship, and re- 
ceived the Ph.D. in 
applied physics in 
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1956. 

From 1946 to 1949, he taught in the 
physics department at the University of 
Shanghai. From 1950 to 1951, he was Chief 
Accountant at the Universal Textile Com- 
pany, Ltd., Hongkong, China. From 1956 to 
1957, he was at Harvard University as a 
Post-Doctoral Fellow. From 1957 to 1958, 
he was an Assistant Professor at the Uni- 
versity of California, Berkeley, Calif. From 
1958 to 1959, he was an Associate Professor 
in Physics at Tunghai University, Taiwan. 
Since 1959, he has been with the Bell Tele- 
phone Labs., Murray Hill, N. J. | 

Dr. Tang is a member of Sigma Xi and 
Phi-Kappa Phi. 
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AIRBORNE INSTRUMENTS LABORATORY FXR, INC. 
A Division of Cutler-Hammer, Inc. 25-26 50th Street, Woodside 77, N.Y. 


Deer Park, Long Island, New York as F 
Precision Microwave Test Equipment, High Power Microwave 


Research and Development Electronics, Microwave Components & Instrumentation 


AIRTRON, INC. ITT FEDERAL LABORATORIES 


A Division of Litton Industries : 
200 East Hanover Ave., Morris Plains, N.J. ee eae oe GvecpNatien 40, Ha: 


Line-of-Sight and Over-the-Horizon Microwave 


Designers and Producers of Complete Line of Systems; Test Equipment and Components 


Microwave Electronic and Aircraft Components 


LITTON INDUSTRIES 
Electron Tube Division 
960 Industrial Rd., San Carlos, Calif. 


Magnetrons, Klystron, TWT's, Noise Sources, BWO's, 
Display Devices, CFA's, Switch Tubes, MM Wave Tubes 


ALFORD MANUFACTURING COMPANY 
299 Atlantic Ave., Boston 10, Mass. 


RF Instruments, Coaxial Components, 
Antennas and Air Navigation Aids 


CASCADE RESEARCH MELABS 
Div. of Lewis & Kaufman Electronics Corp. illvi : 
5245 San Fernando Rd., Los Angeles, Calif. EMU RIES ed ACS SONOS 


R A Dove! <aeaa eeay Ferrite Devi Microwave Components and Systems for Communications, 
esearch, Deve opment, Production: Microwave Ferrite Devices, Countermeasures, Reconnaissance, Radar, Etc. 
Microwave Components & Subsystems 


EIMAC TUBES, EITEL-McCULLOUGH, INC. MICROLAB 
301 Industrial Way, San Carlos, Calif. 570 West Mount Pleasant Ave., Livingston, NJ. 


Microwave Tubes, TWT-VTM-Reflex Klystrons, Designers and Manufacturers of a Complete Line of 
Power Grid Tubes, Amplifier Klystrons Coaxial Microwave Components and Cavity Filters 


ELECTRONIC SPECIALTY CO. MICROWAVE ASSOCIATES, INC. 
5121 San Fernando Rd., Los Angeles 39, Calif. South Avenue, Burlington, Mass. 


Airborne and Ground Antennas, Towers, Waveguides, Waveguide Devices, Microwave Diodes, Magnetrons, 
Microwave Components, Complete Radiating Systems Duplexers, TWT's, Ferrite & Solid-State Devices 
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THE MICROWAVE JOURNAL 
1330 Beacon St., Brookline 46, Mass. 


A Magazine Devoted to the Interests of Engineers 
Working at Microwave Frequencies 


NARDA MICROWAVE CORP. 
Plainview, L.I., New York 


Microwave Test Equipment & Components, Modulators, 
Bolometer & Thermistors, Ferrite Devices 


PHELPS DODGE COPPER PRODUCTS CORP. 
300 Park Ave., New York 22, N.Y. 


Styroflex, Foamflex Air Dielectric Aluminum Sheathed_ 
Coaxial Cables, Connectors and Waveguides 


STEWART ENGINEERING COMPANY 
Santa Cruz, Galif. 


Backward Wave Oscillators & Other TW Type Tubes; 
Controlled Atmosphere Furnaces & Precision Spot Welder 
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SYLVANIA MICROWAVE DEVICE OPERATIONS 
Sylvania Electric Products Inc. 
500 Evelyn Ave., Mountain View, Calif. 


Magnetrons, Klystron, TWT's, BWO's, Ferrite Devices, 
Waveguide Windows, Microwave Diodes 


TECHNICRAFT DIVISION, Electronic Specialty Co. 


Thomaston, Conn. 


Designers, Developers and Producers of Microwave 
Components, Transmission Lines, and Assemblies 


PRD ELECTRONICS, INC. 
202 Tillary St., Brooklyn 1, N.Y. 


Complete Line of Microwave and Electronic Test Equipment 
Waveguide and Coaxial Components 


SAGE LABORATORIES, INC. 
3 Huron Dr., Natick, Mass. 


Microwave Attenuators, Couplers, Crystal Holders, 
Filters, Hybrids, Mixers, Rotary Joints 


SCIENTIFIC-ATLANTA, INC. 
2162 Piedmont Rd., N.E., Atlanta 19, Ga. 


Complete Antenna Pattern Range Instrumentation, 
Special RF and Antenna Systems Development 


WATKINS-JOHNSON COMPANY 
3333 Hillview Ave., Palo Alto, Calif. 


Res., Dev., Mfg. Microwave Electron Devices, TWT's, 
BWO's, Parametric Amplifiers, Microwave Systems 


WEINSCHEL ENGINEERING COMPANY, INC. 
10503 Metropolitan Ave., Kensington, Md. 


Attenuation Standards, Coax Attenuators, Insertion 
Loss Test Sets, Voltage & Power Calibrators 


WHEELER LABORATORIES, INC. 
Great Neck, N.Y. 
Antenna Lab., Smithtowa, N.Y. 


Consulting Services, Research & Development, 
Microwave Antennas & Waveguide Components 
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